This CellML model represents the synaptic coupling between the four different neuronal populations which make up the brain stem respiratory network computational model described by Smith et al. 2007, and modified by Dr Carlo Laing.
This is a CellML1.1 model, the units have been checked and they are consistent. In order to fully reproduce the results of the original published model we need to duplicate each neuron 50 times to represent a population of cells. We also need to introduce heterogeneity between cells to make the model physiologically realistic. This heterogeneity can be introduced by randomly sampling certain parameter values from a population. In addition, initial conditions for each variable should be randomly chosen from a uniform distribution. Finally, the published model descibes three oscillatory mechanisms; namely the three-phase, two-phase, and one-phase rhythms. The CellML model described here represents the two-phase rhythm. Equations for gsynE and gsynI will have to be adjusted to produce the one-, and three-phase rhythms.
Abstract: Mammalian central pattern generators (CPGs) producing rhythmic movements exhibit extremely robust and flexible behavior. Network architectures that enable these features are not well understood. Here we studied organization of the brain stem respiratory CPG. By sequential rostral to caudal transections through the pontine-medullary respiratory network within an in situ perfused rat brain stem-spinal cord preparation, we showed that network dynamics reorganized and new rhythmogenic mechanisms emerged. The normal three-phase respiratory rhythm transformed to a two-phase and then to a one-phase rhythm as the network was reduced. Expression of the three-phase rhythm required the presence of the pons, generation of the two-phase rhythm depended on the integrity of Botzinger and pre-Botzinger complexes and interactions between them, and the one-phase rhythm was generated within the pre-Botzinger complex. Transformation from the three-phase to a two-phase pattern also occurred in intact preparations when chloride-mediated synaptic inhibition was reduced. In contrast to the three-phase and two-phase rhythms, the one-phase rhythm was abolished by blockade of persistent sodium current (I(NaP)). A model of the respiratory network was developed to reproduce and explain these observations. The model incorporated interacting populations of respiratory neurons within spatially organized brain stem compartments. Our simulations reproduced the respiratory patterns recorded from intact and sequentially reduced preparations. Our results suggest that the three-phase and two-phase rhythms involve inhibitory network interactions, whereas the one-phase rhythm depends on I(NaP). We conclude that the respiratory network has rhythmogenic capabilities at multiple levels of network organization, allowing expression of motor patterns specific for various physiological and pathophysiological respiratory behaviors.
|Schematic of the computational model of the brain stem respiratory network. Model includes interacting neuronal populations within major brain stem respiratory compartments (Pons, BotC, pre-BotC, and rVRG). Spheres represent neuronal populations (excitatory, red; inhibitory, blue; motoneuronal, brown); green triangles represent sources of tonic excitatory drives (in pons, RTN/BotC, and pre-BotC compartments) to different neural populations. Excitatory and inhibitory synaptic connections are indicated by arrows and small circles, respectively. Simulated 'transections' (dashed lines) mimic those performed in situ.|
The original paper reference is cited below:
Spatial and functional architecture of the mammalian brain stem respiratory network: a hierarchy of three oscillatory mechanisms, Smith JC, Abdala AP, Koizumi H, Rybak IA, Paton JF, 2007, Journal of Neurophysiology, 98, 3370-3387. PubMed ID: 17913982.