- Auckland Bioengineering Institute, University of Auckland, Auckland, New Zealand
The stomach is extensively innervated by the vagus nerve and the enteric nervous system. The mechanisms through which this innervation affects gastric motility are being unraveled, motivating the first concerted steps towards the incorporation autonomic regulation into computational models of gastric motility. Computational modeling has been valuable in advancing clinical treatment of other organs, such as the heart. However, to date, computational models of gastric motility have made simplifying assumptions about the link between gastric electrophysiology and motility. Advances in experimental neuroscience mean that these assumptions can be reviewed, and detailed models of autonomic regulation can be incorporated into computational models. This review covers these advances, as well as a vision for the utility of computational models of gastric motility. Diseases of the nervous system, such as Parkinson’s disease, can originate from the brain-gut axis and result in pathological gastric motility. Computational models are a valuable tool for understanding the mechanisms of disease and how treatment may affect gastric motility. This review also covers recent advances in experimental neuroscience that are fundamental to the development of physiology-driven computational models. A vision for the future of computational modeling of gastric motility is proposed and modeling approaches employed for existing mathematical models of autonomic regulation of other gastrointestinal organs and other organ systems are discussed.
1. Background physiology
Gastrointestinal (GI) motility is driven by complex interactions between nerves and myogenic electrophysiological mechanisms. Investigations have uncovered the physiological mechanisms of interaction between enteric nerves and interstitial cells of Cajal (ICC) as gastrointestinal pacemaker cells (Furness, 2012, 2022). In silico investigations of gastric motility have been used to bridge the understanding of the interactions between the neural and myogenic components of gastric motility and their implications in pathological developments. In particular, using simulations based on sound biophysical principles and specific experimental data, it is possible to simulate normal and pathological gastric electrophysiology and estimate biomarkers which match these scenarios. However, in these models the role of regulatory mechanisms, like the autonomic nervous system, has been omitted. Simplifying assumptions about the link between gastric motility and electrophysiology have been made in these models, however, with a burgeoning corpus of physiological research it is now possible to revise these assumptions and explore the impacts of autonomic regulation within a mathematical modeling framework (Du et al., 2013b,2018; Cheng et al., 2021).
The core of the ICC network exists in the myenteric (ICC-MY) and submucosal plexus (ICC-SM), with further intramuscular ICCs (ICC-IM) situated within the circular (ICC-CM) and longitudinal muscle layers (ICC-LM). ICC-IM are aligned with the muscle direction in their respective layers. ICC also exist in other GI organs, but their distribution varies. In general, one of the key physiological functions of ICC is to generate the major component of the gastric slow wave which controls rhythmic contractions of GI smooth muscle cells (SMCs) (Komuro, 1999; Kito, 2011). As a result, smooth muscle in the circular and longitudinal muscle exhibit coordinated depolarizations in the presence of ICC-MY, however in the absence of ICC-MY the depolarization of SMCs in the two layers is not coordinated (Huizinga, 2018). Gap junctions between cells of the gastric wall allow the flow of ionic currents between SMC, ICC, and PDGFRα+ cells in the stomach wall, forming a tissue structure that has been termed the SMC, ICC, PDGFRα+ syncytium (SIP syncytium) (Yeoh et al., 2016; Sanders, 2019).
Bioelectrical depolarization of SMC causes an influx of calcium ions and initiates a cascade of actions that lead to the contraction of muscle filaments (Hirst and Edwards, 2004). It is generally understood that gastric slow waves with inputs from the enteric nervous system (ENS) set the basic rhythmic contractions of the stomach. In the intact organ, the origin of the normal entrained slow waves is the pacemaker region along the greater curvature of the stomach. The waves propagate in coherent wavefronts in the aboral direction towards the pylorus. At the cellular level, the ICC network is believed to be responsible for setting the basic rhythmicity of slow waves and ensures that contractions are organized in an annular fashion (Mah et al., 2021). While ICC are known to exhibit decreasing intrinsic frequencies towards the pylorus of the stomach (Kelly et al., 1969), entrainment of slow waves occurs when the depolarization of a group of ICC with higher intrinsic frequency “phase locks” adjoining ICC with lower intrinsic frequencies, which results in a single frequency in the stomach, as well as over extended segments in the intestine (van Helden et al., 2010; Parsons and Huizinga, 2018). Persistent deviations of activities from the natural pacemaker region have been shown as a biomarker of diseases in gastroparesis and chronic nausea and vomiting (O’Grady et al., 2012; Angeli et al., 2015).
The brain-gut axis is used to label the plethora of interactions between GI organs and the brain. These interactions are mediated through the autonomic nervous system, of which the ENS is a part, and through the endocrine system (Lyte and Cryan, 2014). The effects of the brain-gut axis on the transport function of the stomach are covered in this work, but other functions of the stomach, such as secretion, are also influenced by the brain-gut axis. The vagus nerve is the primary, bi-directional communication pathway for the nervous system component of the brain-gut axis. Efferent vagal neurons originate at the brain stem and target cells in peripheral organs, and afferent vagal neurons originate at peripheral organs and target neurons in the brain stem. Observations made by microscopic imaging have shown that efferent vagal neurons release neurotransmitters at ENS ganglia, specifically targeting only a portion of a single ganglion (Powley et al., 2019). The ENS integrates extrinsic vagal input and intrinsic input from ENS interneurons, sensory neurons, and neuronal circuits thereby regulating gastric motility by neurotransmission to effector cells in the stomach wall (Furness et al., 2020). On the other hand, afferent vagal neurons are sensory neurons, chiefly sensing the stretch and strain of the stomach wall (Powley et al., 2019). Sensory neurons also exist within the ENS (Furness, 2012). Sensory neurons form feedback pathways to control stomach function, in conjunction with efferent nerves. Sensory neuronal activity is also a significant aspect of gastric neuro-circuitry but is not covered in detail in this review.
Extrinsic innervation of the stomach is chiefly provided by the vagus nerve though some sympathetic neurons originating in the T6-T9 level (thoracic) of the spinal cord innervate densely innervate gastrointestinal sphincters adjacent to the stomach (Browning and Travagli, 2014). Detailed quantification of the number and size of vagus nerve bundles was conducted by Prechtl and Powley (1985, 1990), with the measured diameter of neurons being consistent with the relatively slow conduction velocities observed in other parasympathetic nerves. Using cell tracer labeling with DiL and in vivo stimulation, Berthoud et al. (1991) deduced the branches of the vagus nerve from which innervation of gastrointestinal organs arises. Thoracic organs, such as the heart and lungs, are also innervated by the vagus nerve but branches to these organs diverge from the main vagus nerve at the cervical level.
The central nervous system influence on the upper gastrointestinal tract, including the stomach and esophagus, is complemented by reflexive control mediated at the level of the brainstem, termed vago-vagal reflexes. Afferent vagus nerve signaling to the nucleus tractus solitarius is integrated by neurocircuitry in the brainstem, and reflexes are relayed via efferent vagus neurons which originate at the dorsal motor nucleus of the vagus nerve. The gastric accommodation reflex results in a relaxation in gastric smooth muscle to accommodate food intake upon the distension of the proximal stomach. Another reflex, the esophagogastric reflex, results in gastric distension in response to stretch of the distal esophagus (Travagli and Anselmi, 2016). The nature of esophagogastric reflex disruption was noted as being similar between vagotomised patients and those diagnosed with functional dyspepsia (Troncon et al., 1995). The authors of the study suggested that vagus nerve defects may be responsible for some effect of functional dyspepsia.
The interface of neural and slow wave electrophysiology, with its complex interactions between multiple cooperating mechanisms, presents a complex challenge to the advance of in silico simulations. A previous review has covered aspects related to neuromodulation modeling of the gut (Barth and Shen, 2018), but there is significant knowledge gap and challenges to modeling the mechanistic components of automatic regulation. Specifically, the detailed interactions between sub-categories of ICC and the ENS. ICC have the ability to generate their own intrinsic slow waves and maintain an entrained wave propagation over an extended area of the GI tract, yet the influence of neural innervation of ICC over the same spatiotemporal scales remains incompletely understood (Forrest et al., 2006; Furness, 2012; Sanders and Ward, 2019). Furthermore, the ENS also has a level of direct control over GI smooth muscles, and it is unclear exactly how innervation interacts with slow waves under different physiological and patho-physiological conditions.
2. Mathematical models
A number of mathematical models have proposed to model various aspects of GI electrophysiology and motility (Du et al., 2010; Mah et al., 2021). The purpose of this section is to review approaches for continuum modeling of the stomach and their application in coupled-autonomic models; other reviews have provided more in-depth coverage of the different gastric function models (Du et al., 2010; Mah et al., 2021).
2.1. Cell electrophysiology models
The gastric slow wave is composed of multiple currents which are thought to originate in different cells (Sanders et al., 2006). Some models, (Edwards and Hirst, 2006; Faville et al., 2008; Youm et al., 2014) include detailed components that describe unitary potentials. A model of intracellular calcium by Means and Sneyd (2010) used stochastic terms to model IP3R activation. Models incorporating a variety of cellular pacemaking mechanisms have also been published. Youm et al. (2006) and Corrias and Buist (2008) modeled pacemaker potential initiation as a result of non-selective cation channel (NSCC) currents. Lees-Green et al. (2014) investigated the role of calcium-dependent chloride channels, encoded by the ANO1 gene, in initiating the pacemaker potential. These models have been used to explore the dynamics of cellular electrophysiology, however, many models omit this level of detail since the characteristics of the slow wave in aggregate rather than its components are relevant for modeling functional outcomes at the organ level. Many of the computed currents and gating variables in biophysical models can be simplified or lumped depending on the application of the model to reduce the computational load when simulating organ level functional behavior. Therefore, simplified models have often been used to model slow waves in multiscale models. Since ICC contribute the major component of the slow wave, models of ICC are often viewed as analogs of slow wave models, particularly when used as part of multiscale simulations. For reference, simulations of selected ICC models are shown in Figure 1.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1146097/fnins-17-1146097-HTML/image_m/fnins-17-1146097-g001.jpg)
Figure 1. Examples of membrane potential simulated by cell models of gastrointestinal ICC. Cell model implementations were retrieved from the Physiome Model Repository. Panel (A) shows the membrane potential in a single gastric ICC-MY cell using the Corrias and Buist (2008) model. Panel (B) shows the membrane potential of a single small intestine ICC-MY using the Youm et al. (2006) model. Panel (C) shows the membrane potential of a single jejunal ICC using the Lees-Green et al. (2014) model.
Du et al. (2013a) developed a simplified model of gastric ICC by incorporating ionic currents from the Corrias & Buist model into a theoretical model originally developed by Imtiaz et al. (2002). The resulting model achieved a trade-off between biophysical relevance and computational load. This was particularly suited to use in multiscale organ models and has been used to compare in vivo recordings with theoretical understandings of slow waves (Wang et al., 2018). Other multiscale organ models use simpler phenomenological models which capture organ-level slow wave dynamics but cannot model cellular ion concentrations at all. For example, an early multi-scale slow wave model used the Aliev et al. (2000) phenomenological cell model to simulate the entrainment of excitable gastric cells (Pullan et al., 2004).
2.2. Whole-organ electrophysiology models
Multiscale modeling extends single-cell models by applying governing equations that spatially couple cellular potential. Continuum (multi-domain) modeling and network modeling are the two main approaches for simulating multi-cellular level events. The bidomain model has widely been used to model the propagation of electrical activity across the heart, where investigations have focused on the detailed relationship between tissue structure and function (Austin et al., 2006). Recent advances in microscopic imaging techniques have also allowed the incorporation of 3D imaging into models from the cellular level to the whole-organ level (Sands et al., 2022).
An early example a whole stomach model (Pullan et al., 2004) used an approach where the organ was represented as a continuum-averaged mesh with mesh elements coupled by governing equations. This was the first model of gastric bioelectrical propagation which used the bidomain model. The cell model, in this case the Aliev et al. (2000) cell model, was solved for the continuum element. The bidomain formulation was used as the governing equation, where a reaction-diffusion equation relates the change in cell potential to the flow of current in space and another equation enforces the conservation of charge (Buist and Poh, 2010).
The bidomain model models the flow of current between an intracellular and extracellular domain. The bidomain model was initially extensively used to model the propagation of bioelectrical events across the heart (Geselowitz and Miller, 1983). As noted in Figure 2, the intracellular domain (subscript i) and extracellular domain (subscript e) have conductivity tensors M and potential u. Equations 1 and 3 describe the bidomain model. The physical basis of the bidomain model is that charge in the two domains is conserved. This means that the change in current density must be equal between the intracellular and extracellular domains. Additionally, the changes in current density must travel through the cell membrane, therefore they are equal to the total membrane current (It). Mathematically this is expressed by
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1146097/fnins-17-1146097-HTML/image_m/fnins-17-1146097-g002.jpg)
Figure 2. Domains and variables in the bidomain model. The intracellular and extracellular domains are linked by a series of ionic currents, which give rise to changes in membrane potentials of ICC and SMC.
where ∇⋅ is the divergence operator and ∇ is the vector gradient operator. Since total membrane current is composed of ionic current (Iion) and current due to membrane capacitance, we can form two equations for the bidomain model,
where χ is the ratio of the surface area to the volume of the integration unit and Cm is the capacitance of the membrane.
Two variations on the bidomain model have been used for modeling of gastric electrophysiology. The monodomain model simplifies the bidomain model by assuming that conductivity tensors of the two domains are linearly proportional, this reduces the number of computations required to solve the model at each time step.
Another approach is to model cells as nodes in a network and couple ionic currents as a transfer of charge between nodes or some extracellular space. The coupling is achieved through gap junction models which are typically Ohmic resistor in a network. A recent paper by Ahmed et al. (2021) demonstrated this for gastric cells but the approach has been used in other biological modeling. For example, models of uterine electrophysiology used a grid lattice network with gap junction coupling to explore organ level bioelectrical propagation (Xu et al., 2015). A previous study has also proposed an extended bidomain or “tridomain” model to simulate the sharing of extracellular space by ICC and SMC (Buist and Poh, 2010; Sathar et al., 2015), which would be relevant given the close association between ENS and ICC in the SIP syncytium. A more detailed derivation of both the bidomain model and extended bidomain model for the gastrointestinal context are presented in Buist and Poh (2010).
2.3. Electromechanical cell models
Mechanical contraction has been modeled (Figure 3) by coupling with electrical potential (Du et al., 2011; Brandstaeter et al., 2018; Klemm et al., 2020, 2023) where the relationship between electrical potential and mechanical contraction is modeled as a continuous one. For successful biomechanical modeling the constitutive parameters that define tissue properties need to be determined. Bauer et al. characterized the biomechanical properties of different regions of the porcine stomach wall (Bauer et al., 2020). Their work showed that the stomach wall was spatially heterogeneous and mechanical properties were anisotropic. Reconciliation of this data with anatomical map of muscle thickness and fibre orientation (Avci et al., 2022; Di Natale et al., 2022), will progress work on biomechanical modeling of the stomach through finite element approaches.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1146097/fnins-17-1146097-HTML/image_m/fnins-17-1146097-g003.jpg)
Figure 3. An electromechanical model of gastric slow waves and contraction published by Klemm et al. (2023) reproduced from Klemm et al. (2023) licenced under CC-BY-NC-ND 4.0. The simulations show the emergence of an ectopic pacemaker of slow waves in the gastric antrum due to stretching and the subsequent change in slow waves and mechanical contractions at t1: 60 s, t2: 62 s, t3: 64 s, t4: 65 s and t5: 69 s.
Another approach is a meshless smooth particle hydrodynamics approach as demonstrated by Alexiadis et al. (2021) where intestinal contents were modeled. The particles representing intestinal contents interacted with a lattice of distinct particles that represented the intestinal wall. The intestinal wall particles were constrained by linear springs. This model simplified the mechanics of the intestinal wall because the problem was primarily concerned with modeling the velocity of intestinal content flow in normal and diseased intestinal ENS health.
2.4. Coupled autonomic-organ models
While no models of stomach motility with vagal neural control have been developed, there are examples in the literature of neural control systems being modeled in combination with motility modeling for other organs of the GI tract. Many of these examples are from the past decade, underlining the novelty of this research domain. Good examples of earlier models of enteric nerve physiology and their impact of motility include the series of models developed by Chambers et al. (Chambers et al., 2008, 2011, 2014a,b), which have also been covered in a separate review article (Chambers et al., 2014b).
The aforementioned model of intestinal contraction and fluid transport by Alexiadis et al. (2021) whose biomechanics approach was described earlier also incorporated the effect of the ENS. Intestinal wall particles were coupled to the output layer of an artificial neural network which enforced the contraction or relaxation of the intestinal wall particle. The computational model implemented a theorized mechanism of peristalsis where the presence of a bolus caused stretching distal to the peristaltic wave and the resulting in a feedback response via the ENS. The ENS was modeled as an artificial neural network, which was trained by reinforcement learning to favor a maximal velocity of intestinal content transport.
Another example is a model of colonic motility by Barth et al. (2017). This model uses multiscale approach, modeling various cell populations separately and connecting them in a network with interactions defined by fitting to experimental measurements. The neural control was implemented as a network of artificial neurons positioned in multiple sub-populations longitudinally along the colon. Sub-populations were positioned such that they resembled the arrangement of neurons observed by imaging of rodent tissue. ICC were simulated using the model proposed by Edwards and Hirst (2006). However, the innervation of ICC by the ENS was not modeled. Instead, only the generation of junction potentials at circular smooth muscle cells was modeled. Neural stimulation in this model resulted in the muscles reaching or failing to reach a contraction threshold when ICC mediated slow wave depolarization occurred.
While there are few GI specific models using coupled neural electrophysiology, smooth muscle interactions with neural stimulation have been modeled for other smooth muscle organs such as the bladder and uterus. The bladder receives central nervous system innervation from the spinal cord. Innervation of the bladder is more direct than the stomach due to the lack of an intrinsic bladder nervous system similar to the ENS. Sympathetic neurons directly synapse with smooth muscle of the bladder to cause excitatory and inhibitory junction potentials (Fowler et al., 2008). The existence of interstitial cells in the bladder is well established, but their role in neurotransmission is an area of ongoing research (Koh et al., 2018).
Computational models of the propagation of electrophysiological activity have been published for both organs. In these examples, published by Appukuttan et al. (2015) and Xu et al. (2015) for the bladder and uterus, respectively, cell-cell coupling by gap junctions was modeled as an Ohmic resistor. Appukuttan et al. (2015) modeled bladder SMCs as passive single compartment cells which responded to external electrical stimulation, including stimulation occurring as a result of synaptic input from neurons. Results showed the time course of the spread of junction potentials induced by synaptic transmission.
In contrast, the model of uterus electrophysiology by Xu et al. (2015) uses a biophysical model of SMCs to demonstrate entrainment of SMC electrophysiological activity in the uterus. This model does not incorporate any neurotransmission, but it does model passive interstitial cells. The authors found that their model was able to simulate entrainment without a defined pacemaker region only if there was sufficient connectivity between adjacent cells.
Finally, Dokos et al. (1996) published a model of neurotransmitter release at the sinoatrial node. The mathematical model published by Dokos et al. (1996) was very detailed, and suitable for theoretical investigation at a cellular scale. The results of their simulations suggested that sustained negative chronotropic effects of cholinergic stimulation were not the result of the time course of neurotransmitter. Instead, the model results suggested that potassium currents controlled the longevity of the chronotropic response. The incorporation of a detailed biophysical model such as this one directly into a whole organ, multi-scale model would be computationally difficult. However, analysis of complex biophysical models can show how they can be simplified to yield new mathematical models suitable for multi-scale computational modeling. The simplified biophysically-based model of ICC gastric slow wave activity by Du et al. (2013a) is an example of this approach.
3. Future directions
3.1. Models of autonomic control of the stomach
There is clear experimental evidence that autonomic control of gastric motility occurs through the ENS, ICC, and SMC. The key gap in knowledge exists in identifying the details of these how various components, from the network connectivity of cells, to the expression of particular cells, and regional differences within the stomach, interact to bring about a complex and coordinated control of motility under various physiological stimuli. A comprehensive evaluation of these interactions is critical, as every component is a potential target of disease biomarkers or therapies, which may have cascading effects on the tightly co-regulatory system of control. Mirroring the availability of experimental data, a number of biophysically based cell models of ICC and SMC have been developed in recent decades, with further attempts to link them to biomechanical contractions. These investigations provide the critical building blocks for incorporation of neural regulation of GI motility, which has remained relatively crude by assuming direct innervation of SMC without interactions with ICC. The main remaining challenge is a lack of consistency in studies specifically investigating interactions between the ENS, ICC, and SMC and an understanding of the variation in these interactions between functional regions of the stomach and between animal models.
3.2. Applications of mathematical models
One potential application of mathematical models is to understand the integrated mechanisms of diseases related to the brain-gut axis. For example, it has been highlighted that in Parkinson’s disease, propagation of pathological α-synuclein aggregates throughout the nervous system, and α-synuclein propagation along the brain-gut axis is a potential disease initiation mechanism (i.e., Braak’s hypothesis) (Braak et al., 2003). There is significant evidence that the vagus and enteric nerves facilitate this propagation in the central nervous system and from the gut, respectively (Fasano et al., 2015). Such changes in the enteric nerves have been proposed as a potential biomarker and treatment target of Parkinson’s disease. However, the details of how neuromodulation influences and normalizes gastric slow wave dysrhythmias remain unclear. An aspect that the modeling can contribute is to develop a detailed structural-functional (neuro-electromechanical) model of stomach, with a vagal neurocircuit model complementing the existing whole-organ gastric slow wave model. The permutations of stimulation protocols and their effects of changes of gastric slow waves and motility can then be tested. This will enable future in silico hypothesis testing for optimal neuromodulation parameters targeted at specific symptoms caused by Parkinson’s disease.
Secondary clinical applications can include generating hypothetical gastric signatures for far-field body-surface mapping investigations (Calder et al., 2022). The whole-organ models representing various disease states and responses to neuromodulation parameters can be placed inside a volume conductor to simulate resultant potentials on the body-surface, which can be detected clinically using far-field devices such as body-surface gastric mapping (BSGM) or a superconducting quantum interference device (SQUID) (Kim et al., 2010; Bradshaw et al., 2016; Calder et al., 2022). Given the relatively low signal-to-noise ratio of gastric slow waves compared to cardiac activity, synthetic signals have the advantage of being noise-free so that the true response of the stomach can be explored, and a refined target can be generated for analysis of real signals. Another potential application of the models is to understand the impact of vagotomy on gastric slow waves and motility, which has been shown to significantly alter neuroendocrine peptide levels in the GI tract (El-Salhy et al., 2000), and potentially impacts the response of the stomach to neuromodulation.
One clinical application of computational physiology models is in closed loop control of neurostimulation devices (Payne et al., 2019). Closed-loop control requires the stimulation system to measure a functional quantity and predict the stimulation parameters that will bring the functional quantity towards a desirable set point. The predictive model required for this needs to reach a suitable tradeoff for accuracy and computation speed. The approach taken by Branen et al. (2022) to developing such a model for the cardiovascular system was to use outputs from a mathematical model as training data for several neural networks. The neural networks were used to predict heart rate and mean arterial blood pressure outcomes over hundreds of cardiac cycles. The trained hyperparameters in the neural network model are well suited for the fast prediction of cardiovascular function given inputs for stimulation parameters (pulse width and frequency at three stimulation locations), and the current measured mean arterial pressure and heart rate. The feasibility of this approach for gastric clinical applications is contingent on the availability of a high-quality, predictive, and accurate mathematical model with which to generate valid training data. At present, neuromodulation of the brain-gut axis is receiving active interest, with studies that have reported recovery of gastric functions in functional dyspepsia patients (Zhu et al., 2021) as well as new metrics from non-invasive body-surface gastric mapping studies for definitive of normative values and classifications of diseases (Gharibans et al., 2022; Varghese et al., 2022). Together, these emerging metrics can be used as objective functions for brain-gut axis models to generate targets to device closed-loop protocols for controlling gastric functions via neuromodulation. Closed-loop control of gastric pacing, where electrical stimulation is delivered to the smooth muscle syncytium rather than nerves, has been demonstrated using a model of gastric slow wave propagation (Wang et al., 2022).
4. Conclusion
In conclusion, gastric functions are highly dependent on activity from the vagus nerve and enteric nervous system, yet the exact interactions between the gut and brain remain under-investigated. Advances in experimental and computational techniques will lead the development of predictive multi-scale models that can be used to explore the impacts of various pathological conditions related to the brain-gut axis and generate potential treatment targets for the next generation of neuromodulation devices.
Author contributions
OA, RA, LC, and PD were involved in all aspects of manuscript preparation. All authors contributed to the article and approved the submitted version.
Funding
This work was supported by the Marsden Council Fund managed by the Royal Society Te Apārangi, NIH SPARC award (1OT2OD030538-01) and the Ministry of Business, Innovation and Employment’s Catalyst: Strategic fund. OA was supported by a University of Auckland Doctoral Scholarship.
Acknowledgments
The authors thank Drs. John Furness and Zhongming Liu for their expert guidance and inputs on physiology and experimental recording related to this review.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
References
Ahmed, M., Venugopal, S., and Jung, R. (2021). Engaging biological oscillators through second messenger pathways permits emergence of a robust gastric slow-wave during peristalsis. PLoS Comput. Biol. 17:e1009644. doi: 10.1371/journal.pcbi.1009644
Alexiadis, A., Simmons, M., Stamatopoulos, K., Batchelor, H., and Moulitsas, I. (2021). The virtual physiological human gets nerves! How to account for the action of the nervous system in multiphysics simulations of human organs. J. R. Soc. Interface 18:20201024. doi: 10.1098/rsif.2020.1024
Aliev, R., Richards, W., and Wikswo, J. P. (2000). A simple nonlinear model of electrical activity in the intestine. J. Theor. Biol. 204, 21–28. doi: 10.1006/jtbi.2000.1069
Angeli, T., Cheng, L., Du, P., Wang, T., Bernard, C., Vannucchi, M., et al. (2015). Loss of interstitial cells of Cajal and patterns of gastric dysrhythmia in patients with chronic unexplained nausea and vomiting. Gastroenterology 149, 56–66.e5. doi: 10.1053/j.gastro.2015.04.003
Appukuttan, S., Brain, K., and Manchanda, R. (2015). A computational model of urinary bladder smooth muscle syncytium : Validation and investigation of electrical properties. J. Comput. Neurosci. 38, 167–187. doi: 10.1007/s10827-014-0532-6
Austin, T., Hooks, D., Hunter, P., Nickerson, D., Pullan, A., Sands, G., et al. (2006). Modeling cardiac electrical activity at the cell and tissue levels. Ann. N. Y. Acad. Sci. 1080, 334–347. doi: 10.1196/annals.1380.025
Avci, R., Wickens, J., Sangi, M., Athavale, O., Di Natale, M., Furness, J., et al. (2022). A computational model of biophysical properties of the rat stomach informed by comprehensive analysis of muscle anatomy. Annu. Int. Conf. IEEE Eng. Med. Biol. Soc. 2022, 4954–4957. doi: 10.1109/EMBC48229.2022.9871314
Barth, B. B., Henriquez, C., Grill, W., and Shen, X. (2017). Electrical stimulation of gut motility guided by an in silico model. J. Neural. Eng. 14:066010. doi: 10.1088/1741-2552/aa86c8
Barth, B. B., and Shen, X. (2018). Computational motility models of neurogastroenterology and neuromodulation. Brain Res. 1693(Pt B), 174–179. doi: 10.1016/j.brainres.2018.02.038
Bauer, M., Morales-Orcajo, E., Klemm, L., Seydewitz, R., Fiebach, V., Siebert, T., et al. (2020). Biomechanical and microstructural characterisation of the porcine stomach wall: Location- and layer-dependent investigations. Acta Biomater. 102, 83–99. doi: 10.1016/j.actbio.2019.11.038
Berthoud, H., Carlson, N., and Powley, T. (1991). Topography of efferent vagal innervation of the rat gastrointestinal tract. Am. J. Physiol. 260, R200–R207. doi: 10.1152/ajpregu.1991.260.1.R200
Braak, H., Rüb, U., Gai, W., and Del Tredici, K. (2003). Idiopathic Parkinson’s disease: Possible routes by which vulnerable neuronal types may be subject to neuroinvasion by an unknown pathogen. J. Neural. Transm. 110, 517–536. doi: 10.1007/s00702-002-0808-2
Bradshaw, L., Cheng, L., Chung, E., Obioha, C., Erickson, J., Gorman, B., et al. (2016). Diabetic gastroparesis alters the biomagnetic signature of the gastric slow wave. Neurogastroenterol. Motil. 28, 837–848. doi: 10.1111/nmo.12780
Brandstaeter, S., Gizzi, A., Fuchs, S. L., Gabauer, A., Aydin, R. C., and Cyron, C. J. (2018). Computational model of gastric motility with active-strain electromechanics. Z. Angew. Math. Mech. 98, 2177–2197. doi: 10.1002/zamm.201800166
Branen, A., Yao, Y., Kothare, M., Mahmoudi, B., and Kumar, G. (2022). Data driven control of vagus nerve stimulation for the cardiovascular system: An in silico computational study. Front. Physiol. 13:798157. doi: 10.3389/fphys.2022.798157
Browning, K., and Travagli, R. (2014). Central nervous system control of gastrointestinal motility and secretion and modulation of gastrointestinal functions. Compr. Physiol. 4, 1339–1368. doi: 10.1002/cphy.c130055
Buist, M., and Poh, Y. (2010). An extended bidomain framework incorporating multiple cell types. Biophys. J. 99, 13–18. doi: 10.1016/j.bpj.2010.03.054
Calder, S., Cheng, L., Andrews, C., Paskaranandavadivel, N., Waite, S., Alighaleh, S., et al. (2022). Validation of noninvasive body-surface gastric mapping for detecting gastric slow-wave spatiotemporal features by simultaneous serosal mapping in porcine. Am. J. Physiol. Gastrointest. Liver Physiol. 323, G295–G305. doi: 10.1152/ajpgi.00049.2022
Chambers, J., Bornstein, J., and Thomas, E. (2008). Insights into mechanisms of intestinal segmentation in guinea pigs: A combined computational modeling and in vitro study. Am. J. Physiol. Gastrointest. Liver Physiol. 295, G534–G541. doi: 10.1152/ajpgi.90303.2008
Chambers, J., Bornstein, J., and Thomas, E. (2011). Multiple neural oscillators and muscle feedback are required for the intestinal fed state motor program. PLoS One 6:e19597. doi: 10.1371/journal.pone.0019597
Chambers, J., Bornstein, J., Gwynne, R., Koussoulas, K., and Thomas, E. A. (2014a). A detailed, conductance-based computer model of intrinsic sensory neurons of the gastrointestinal tract. Am. J. Physiol. Gastrointest. Liver Physiol. 307, G517–G532. doi: 10.1152/ajpgi.00228.2013
Chambers, J., Thomas, E., and Bornstein, J. (2014b). Mathematical modelling of enteric neural motor patterns. Clin. Exp. Pharmacol. Physiol. 41, 155–164. doi: 10.1111/1440-1681.12209
Cheng, L., Nagahawatte, N., Avci, R., Du, P., Liu, Z., and Paskaranandavadivel, N. (2021). Strategies to refine gastric stimulation and pacing protocols: Experimental and modeling approaches. Front. Neurosci. 15:645472. doi: 10.3389/fnins.2021.645472
Corrias, A., and Buist, M. (2008). Quantitative cellular description of gastric slow wave activity. Am. J. Physiol. Gastrointest. Liver Physiol. 294, G989–G995. doi: 10.1152/ajpgi.00528.2007
Di Natale, M., Patten, L., Molero, J., Stebbing, M., Hunne, B., Wang, X., et al. (2022). Organisation of the musculature of the rat stomach. J. Anat. 240, 711–723. doi: 10.1111/joa.13587
Dokos, S., Celler, B., and Lovell, N. (1996). Vagal control of sinoatrial rhythm: A mathematical model. J. Theor. Biol. 182, 21–44. doi: 10.1006/jtbi.1996.0141
Du, P., Calder, S., Angeli, T., Sathar, S., Paskaranandavadivel, N., O’Grady, G., et al. (2018). Progress in mathematical modeling of gastrointestinal slow wave abnormalities. Front. Physiol. 8:1136. doi: 10.3389/fphys.2017.01136
Du, P., O’Grady, G., Davidson, J., Cheng, L., and Pullan, A. (2010). Multiscale modeling of gastrointestinal electrophysiology and experimental validation. Crit. Rev. Biomed. Eng. 38, 225–254. doi: 10.1615/critrevbiomedeng.v38.i3.10
Du, P., O’Grady, G., Gao, J., Sathar, S., and Cheng, L. (2013b). Toward the virtual stomach: Progress in multiscale modeling of gastric electrophysiology and motility. Wiley Interdiscip. Rev. Syst. Biol. Med. 5, 481–493. doi: 10.1002/wsbm.1218
Du, P., Gao, J., O’Grady, G., and Cheng, L. K. (2013a). A simplified biophysical cell model for gastric slow wave entrainment simulation. Annu. Int. Conf. IEEE Eng. Med. Biol. Soc. 2013, 6547–6550. doi: 10.1109/EMBC.2013.6611055
Du, P., Poh, Y., Lim, J., Gajendiran, V., O’Grady, G., Buist, M., et al. (2011). A preliminary model of gastrointestinal electromechanical coupling. IEEE Trans. Biomed. Eng. 58, 3491–3495. doi: 10.1109/TBME.2011.2166155
Edwards, F., and Hirst, G. (2006). An electrical analysis of slow wave propagation in the guinea-pig gastric antrum. J. Physiol. 571, 179–189. doi: 10.1113/jphysiol.2005.100743
El-Salhy, M., Danielsson, A., Axelsson, H., and Qian, B. (2000). Neuroendocrine peptide levels in the gastrointestinal tract of mice after unilateral cervical vagotomy. Regul. Pept. 88, 15–20. doi: 10.1016/s0167-0115(99)00107-x
Fasano, A., Visanji, N., Liu, L., Lang, A., and Pfeiffer, R. (2015). Gastrointestinal dysfunction in Parkinson’s disease. Lancet Neurol. 14, 625–639. doi: 10.1016/S1474-4422(15)00007-1
Faville, R., Pullan, A., Sanders, K., and Smith, N. P. (2008). A biophysically based mathematical model of unitary potential activity in interstitial cells of Cajal. Biophys. J. 95, 88–104. doi: 10.1529/biophysj.107.122507
Forrest, A., Ordög, T., and Sanders, K. (2006). Neural regulation of slow-wave frequency in the murine gastric antrum. Am. J. Physiol. Gastrointest. Liver Physiol. 290, G486–G495. doi: 10.1152/ajpgi.00349.2005
Fowler, C., Griffiths, D., and de Groat, W. (2008). The neural control of micturition. Nat. Rev. Neurosci. 9, 453–466. doi: 10.1038/nrn2401
Furness, J. (2012). The enteric nervous system and neurogastroenterology. Nat. Rev. Gastroenterol. Hepatol. 9, 286–294. doi: 10.1038/nrgastro.2012.32
Furness, J. (2022). Comparative and evolutionary aspects of the digestive system and its enteric nervous system control. Adv. Exp. Med. Biol. 1383, 165–177. doi: 10.1007/978-3-031-05843-1_16
Furness, J., Di Natale, M., Hunne, B., Oparija-Rogenmozere, L., Ward, S., Sasse, K., et al. (2020). The identification of neuronal control pathways supplying effector tissues in the stomach. Cell Tissue Res. 382, 433–445. doi: 10.1007/s00441-020-03294-7
Geselowitz, D., and Miller, W. T. (1983). A bidomain model for anisotropic cardiac muscle. Ann. Biomed. Eng. 11, 191–206. doi: 10.1007/BF02363286
Gharibans, A., Calder, S., Varghese, C., Waite, S., Schamberg, G., Daker, C., et al. (2022). Gastric dysfunction in patients with chronic nausea and vomiting syndromes defined by a noninvasive gastric mapping device. Sci. Transl. Med. 14:eabq3544. doi: 10.1126/scitranslmed.abq3544
Hirst, G., and Edwards, F. (2004). Role of interstitial cells of Cajal in the control of gastric motility. J. Pharmacol. Sci. 96, 1–10. doi: 10.1254/jphs.crj04002x
Huizinga, J. D. (2018). “The physiology and pathophysiology of interstitial cells of cajal: Pacemaking, innervation, and stretch sensation,” in Physiology of the gastrointestinal tract: Sixth edition, eds H. Said, J. K. Kaunitz, F. Ghishan, J. Merchant, and J. Wood (Amsterdam: Elsevier), 305–335. doi: 10.1016/B978-0-12-809954-4.00013-X
Imtiaz, M., Smith, D., and van Helden, D. F. (2002). A theoretical model of slow wave regulation using voltage-dependent synthesis of inositol 1,4,5-trisphosphate. Biophys. J. 83, 1877–1890. doi: 10.1016/S0006-3495(02)73952-0
Kelly, K., Code, C., and Elveback, L. (1969). Patterns of canine gastric electrical activity. Am. J. Physiol. 217, 461–470. doi: 10.1152/ajplegacy.1969.217.2.461
Kim, J., Bradshaw, L., Pullan, A., and Cheng, L. (2010). Characterization of gastric electrical activity using magnetic field measurements: A simulation study. Ann. Biomed. Eng. 38, 177–186. doi: 10.1007/s10439-009-9804-0
Kito, Y. (2011). The functional role of intramuscular interstitial cells of Cajal in the stomach. J. Smooth Muscle Res. 47, 47–53. doi: 10.1540/jsmr.47.47
Klemm, L., Seydewitz, R., Borsdorf, M., Siebert, T., and Böl, M. (2020). On a coupled electro-chemomechanical model of gastric smooth muscle contraction. Acta Biomater. 109, 163–181. doi: 10.1016/j.actbio.2020.04.007
Klemm, L., Seydewitz, R., Siebert, T., and Böl, M. (2023). Three-dimensional multi-field modelling of gastric arrhythmias and their effects on antral contractions. Comput. Biol. Med. 153:106488. doi: 10.1016/j.compbiomed.2022.106488
Koh, S., Lee, H., Ward, S., and Sanders, K. (2018). The mystery of the interstitial cells in the urinary bladder. Annu. Rev. Pharmacol. Toxicol. 58, 603–623. doi: 10.1146/annurev-pharmtox-010617-052615
Komuro, T. (1999). Comparative morphology of interstitial cells of Cajal: Ultrastructural characterization. Microsc. Res. Tech. 47, 267–285.
Lees-Green, R., Gibbons, S., Farrugia, G., Sneyd, J., and Cheng, L. (2014). Computational modeling of anoctamin 1 calcium-activated chloride channels as pacemaker channels in interstitial cells of Cajal. Am. J. Physiol. Gastrointest. Liver Physiol. 306, G711–G727. doi: 10.1152/ajpgi.00449.2013
Lyte, M., and Cryan, J. F. (2014). “Microbial endocrinology: The microbiota-gut-brain axis in health and disease,” in Advances in experimental medicine and biology, eds M. Lyte and J. F. Cryan (New York, NY: Springer), doi: 10.1007/978-1-4939-0897-4
Mah, S., Avci, R., Cheng, L., and Du, P. (2021). Current applications of mathematical models of the interstitial cells of Cajal in the gastrointestinal tract. Wires Mech. Dis. 13:e1507. doi: 10.1002/wsbm.1507
Means, S., and Sneyd, J. (2010). Spatio-temporal calcium dynamics in pacemaking units of the interstitial cells of Cajal. J. Theor. Biol. 267, 137–152. doi: 10.1016/j.jtbi.2010.08.008
O’Grady, G., Angeli, T., Du, P., Lahr, C., Lammers, W., Windsor, J., et al. (2012). Abnormal initiation and conduction of slow-wave activity in gastroparesis, defined by high-resolution electrical mapping. Gastroenterology 143, 589–598.e3. doi: 10.1053/j.gastro.2012.05.036
Parsons, S., and Huizinga, J. (2018). Slow wave contraction frequency plateaux in the small intestine are composed of discrete waves of interval increase associated with dislocations. Exp. Physiol. 103, 1087–1100. doi: 10.1113/EP086871
Payne, S., Furness, J., and Stebbing, M. (2019). Bioelectric neuromodulation for gastrointestinal disorders: Effectiveness and mechanisms. Nat. Rev. Gastroenterol. Hepatol. 16, 89–105. doi: 10.1038/s41575-018-0078-6
Powley, T., Jaffey, D., McAdams, J., Baronowsky, E., Black, D., Chesney, L., et al. (2019). Vagal innervation of the stomach reassessed: Brain-gut connectome uses smart terminals. Ann. N. Y. Acad. Sci. 1454, 14–30. doi: 10.1111/nyas.14138
Prechtl, J., and Powley, T. (1985). Organization and distribution of the rat subdiaphragmatic vagus and associated paraganglia. J. Comp. Neurol. 235, 182–195. doi: 10.1002/cne.902350204
Prechtl, J., and Powley, T. (1990). The fiber composition of the abdominal vagus of the rat. Anat. Embryol. 181, 101–115. doi: 10.1007/BF00198950
Pullan, A., Cheng, L., Yassi, R., and Buist, M. (2004). Modelling gastrointestinal bioelectric activity. Prog. Biophys. Mol. Biol. 85, 523–550. doi: 10.1016/j.pbiomolbio.2004.02.003
Sanders, K. (2019). Spontaneous electrical activity and rhythmicity in gastrointestinal smooth muscles. Adv. Exp. Med. Biol. 1124, 3–46. doi: 10.1007/978-981-13-5895-1_1
Sanders, K., and Ward, S. (2019). Nitric oxide and its role as a non-adrenergic, non-cholinergic inhibitory neurotransmitter in the gastrointestinal tract. Br. J. Pharmacol. 176, 212–227. doi: 10.1111/bph.14459
Sanders, K., Koh, S., and Ward, S. (2006). Interstitial cells of cajal as pacemakers in the gastrointestinal tract. Annu. Rev. Physiol. 68, 307–343. doi: 10.1146/annurev.physiol.68.040504.094718
Sands, G., Ashton, J., Trew, M., Baddeley, D., Walton, R., Benoist, D., et al. (2022). It’s clearly the heart! Optical transparency, cardiac tissue imaging, and computer modelling. Prog. Biophys. Mol. Biol. 168, 18–32. doi: 10.1016/j.pbiomolbio.2021.06.005
Sathar, S., Trew, M., OGrady, G., and Cheng, L. K. (2015). A multiscale tridomain model for simulating bioelectric gastric pacing. IEEE Trans. Biomed. Eng. 62, 2685–2692. doi: 10.1109/TBME.2015.2444384
Travagli, R., and Anselmi, L. (2016). Vagal neurocircuitry and its influence on gastric motility. Nat. Rev. Gastroenterol. Hepatol. 13, 389–401. doi: 10.1038/nrgastro.2016.76
Troncon, L., Thompson, D., Ahluwalia, N., Barlow, J., and Heggie, L. (1995). Relations between upper abdominal symptoms and gastric distension abnormalities in dysmotility like functional dyspepsia and after vagotomy. Gut 37, 17–22. doi: 10.1136/gut.37.1.17
van Helden, D., Laver, D., Holdsworth, J., and Imtiaz, M. (2010). Generation and propagation of gastric slow waves. Clin. Exp. Pharmacol. Physiol. 37, 516–524. doi: 10.1111/j.1440-1681.2009.05331.x
Varghese, C., Schamberg, G., Calder, S., Waite, S., Carson, D., Foong, D., et al. (2022). Normative values for body surface gastric mapping evaluations of gastric motility using gastric alimetry: Spectral analysis. Am. J. Gastroenterol. [Epub ahead of print]. doi: 10.14309/ajg.0000000000002077
Wang, L., Malik, A., Roop, P., Cheng, L., and Paskaranandavadivel, N. (2022). A framework for the design of a closed-loop gastric pacemaker for treating conduction block. Comput. Methods Programs Biomed. 216:106652. doi: 10.1016/j.cmpb.2022.106652
Wang, T., Du, P., Angeli, T., Paskaranandavadivel, N., Erickson, J., Abell, T., et al. (2018). Relationships between gastric slow wave frequency, velocity, and extracellular amplitude studied by a joint experimental-theoretical approach. Neurogastroenterol. Motil. 30:e13152. doi: 10.1111/nmo.13152
Xu, J., Menon, S., Singh, R., Garnier, N., Sinha, S., and Pumir, A. (2015). The role of cellular coupling in the spontaneous generation of electrical activity in uterine tissue. PLoS One 10:e0118443. doi: 10.1371/journal.pone.0118443
Yeoh, J., Corrias, A., and Buist, M. L. (2016). A mechanistic model of a PDGFRα(+) cell. J. Theor. Biol. 408, 127–136. doi: 10.1016/j.jtbi.2016.08.004
Youm, J., Kim, N., Han, J., Kim, E., Joo, H., Leem, C., et al. (2006). A mathematical model of pacemaker activity recorded from mouse small intestine. Philos. Trans. A Math. Phys. Eng. Sci. 364, 1135–1154. doi: 10.1098/rsta.2006.1759
Youm, J., Leem, C., Lee, S., Song, I., Kim, H., Heo, H., et al. (2014). Modeling of stochastic behavior of pacemaker potential in interstitial cells of Cajal. Prog. Biophys. Mol. Biol. 116, 56–69. doi: 10.1016/j.pbiomolbio.2014.09.002
Keywords: gastroenterology, brain-gut axis, multi-scale modeling, electromechanical modeling, enteric nerves, vagus nerve, electrophysiology
Citation: Athavale ON, Avci R, Cheng LK and Du P (2023) Computational models of autonomic regulation in gastric motility: Progress, challenges, and future directions. Front. Neurosci. 17:1146097. doi: 10.3389/fnins.2023.1146097
Received: 17 January 2023; Accepted: 27 February 2023;
Published: 15 March 2023.
Edited by:
Joel C. Bornstein, The University of Melbourne, AustraliaReviewed by:
Lakshmi Kuttippurathu, Lexicon Pharmaceuticals, United StatesBin Xu, Nanjing University of Chinese Medicine, China
Copyright © 2023 Athavale, Avci, Cheng and Du. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Peng Du, peng.du@auckland.ac.nz