Operation of spinal sensorimotor circuits controlling phase durations during tied-belt and split-belt locomotion after a lateral thoracic hemisection.

Locomotion is controlled by spinal circuits that interact with supraspinal drives and sensory feedback from the limbs. These sensorimotor interactions are disrupted following spinal cord injury. The thoracic lateral hemisection represents an experimental model of an incomplete spinal cord injury, where connections between the brain and spinal cord are abolished on one side of the cord. To investigate the effects of such an injury on the operation of the spinal locomotor network, we used our computational model of cat locomotion recently published in eLife (Rybak et al., 2024) to investigate and predict changes in cycle and phase durations following a thoracic lateral hemisection during treadmill locomotion in tied-belt (equal left-right speeds) and split-belt (unequal left-right speeds) conditions. In our simulations, the "hemisection" was always applied to the right side. Based on our model, we hypothesized that following hemisection, the contralesional ("intact", left) side of the spinal network is mostly controlled by supraspinal drives, whereas the ipsilesional ("hemisected", right) side is mostly controlled by somatosensory feedback. We then compared the simulated results with those obtained during experiments in adult cats before and after a mid-thoracic lateral hemisection on the right side in the same locomotor conditions. Our experimental results confirmed many effects of hemisection on cat locomotion predicted by our simulations. We show that having the ipsilesional hindlimb step on the slow belt, but not the fast belt, during split-belt locomotion substantially reduces the effects of lateral hemisection. The model provides explanations for changes in temporal characteristics of hindlimb locomotion following hemisection based on altered interactions between spinal circuits, supraspinal drives, and somatosensory feedback.

Operation regimes of spinal circuits controlling locomotion and role of supraspinal drives and sensory feedback.

Locomotion in mammals is directly controlled by the spinal neuronal network, operating under the control of supraspinal signals and somatosensory feedback that interact with each other. However, the functional architecture of the spinal locomotor network, its operation regimes, and the role of supraspinal and sensory feedback in different locomotor behaviors, including at different speeds, remain unclear. We developed a computational model of spinal locomotor circuits receiving supraspinal drives and limb sensory feedback that could reproduce multiple experimental data obtained in intact and spinal-transected cats during tied-belt and split-belt treadmill locomotion. We provide evidence that the spinal locomotor network operates in different regimes depending on locomotor speed. In an intact system, at slow speeds (< 0.4 m/s), the spinal network operates in a non-oscillating state-machine regime and requires sensory feedback or external inputs for phase transitions. Removing sensory feedback related to limb extension prevents locomotor oscillations at slow speeds. With increasing speed and supraspinal drives, the spinal network switches to a flexor-driven oscillatory regime and then to a classical half-center regime. Following spinal transection, the model predicts that the spinal network can only operate in the state-machine regime. Our results suggest that the spinal network operates in different regimes for slow exploratory and fast escape locomotor behaviors, making use of different control mechanisms.

Spinal control of locomotion before and after spinal cord injury.

Thoracic spinal cord injury affects long propriospinal neurons that interconnect the cervical and lumbar enlargements. These neurons are crucial for coordinating forelimb and hindlimb locomotor movements in a speed-dependent manner. However, recovery from spinal cord injury is usually studied over a very limited range of speeds that may not fully expose circuitry dysfunction. To overcome this limitation, we investigated overground locomotion in rats trained to move over an extended distance with a wide range of speeds both pre-injury and after recovery from thoracic hemisection or contusion injuries. In this experimental context, intact rats expressed a speed-dependent continuum of alternating (walk and trot) and non-alternating (canter, gallop, half-bound gallop, and bound) gaits. After a lateral hemisection injury, rats recovered the ability to locomote over a wide range of speeds but lost the ability to use the highest-speed gaits (half-bound gallop and bound) and predominantly used the limb contralateral to the injury as lead during canter and gallop. A moderate contusion injury caused a greater reduction in maximal speed, loss of all non-alternating gaits, and emergence of novel alternating gaits. These changes resulted from weak fore-hind coupling together with appropriate control of left-right alternation. After hemisection, animals expressed a subset of intact gaits with appropriate interlimb coordination even on the side of the injury, where the long propriospinal connections were severed. These observations highlight how investigating locomotion over the full range of speeds can reveal otherwise hidden aspects of spinal locomotor control and post-injury recovery.

Spinal control of locomotion before and after spinal cord injury.

Thoracic spinal cord injury affects long propriospinal neurons that interconnect the cervical and lumbar enlargements. These neurons are crucial for coordinating forelimb and hindlimb locomotor movements in a speed-dependent manner. However, recovery from spinal cord injury is usually studied over a very limited range of speeds that may not fully expose circuitry dysfunction. To overcome this limitation, we investigated overground locomotion in rats trained to move over an extended distance with a wide range of speeds both pre-injury and after recovery from thoracic hemisection or contusion injuries. In this experimental context, intact rats expressed a speed-dependent continuum of alternating (walk and trot) and non-alternating (canter, gallop, half-bound gallop, and bound) gaits. After a lateral hemisection injury, rats recovered the ability to locomote over a wide range of speeds but lost the ability to use the highest-speed gaits (half-bound gallop and bound) and predominantly used the limb contralateral to the injury as lead during canter and gallop. A moderate contusion injury caused a greater reduction in maximal speed, loss of all non-alternating gaits, and emergence of novel alternating gaits. These changes resulted from weak fore-hind coupling together with appropriate control of left-right alternation. After hemisection, animals expressed a subset of intact gaits with appropriate interlimb coordination even on the side of the injury, where the long propriospinal connections were severed. These observations highlight how investigating locomotion over the full range of speeds can reveal otherwise hidden aspects of spinal locomotor control and post-injury recovery.

Ipsilateral and Contralateral Interactions in Spinal Locomotor Circuits Mediated by V1 Neurons: Insights from Computational Modeling.

We describe and analyze a computational model of neural circuits in the mammalian spinal cord responsible for generating and shaping locomotor-like oscillations. The model represents interacting populations of spinal neurons, including the neurons that were genetically identified and characterized in a series of previous experimental studies. Here, we specifically focus on the ipsilaterally projecting V1 interneurons, their possible role in the spinal locomotor circuitry, and their involvement in the generation of locomotor oscillations. The proposed connections of these neurons and their involvement in different neuronal pathways in the spinal cord allow the model to reproduce the results of optogenetic manipulations of these neurons under different experimental conditions. We suggest the existence of two distinct populations of V1 interneurons mediating different ipsilateral and contralateral interactions within the spinal cord. The model proposes explanations for multiple experimental data concerning the effects of optogenetic silencing and activation of V1 interneurons on the frequency of locomotor oscillations in the intact cord and hemicord under different experimental conditions. Our simulations provide an important insight into the organization of locomotor circuitry in the mammalian spinal cord.

The role of V3 neurons in speed-dependent interlimb coordination during locomotion in mice.

Speed-dependent interlimb coordination allows animals to maintain stable locomotion under different circumstances. The V3 neurons are known to be involved in interlimb coordination. We previously modeled the locomotor spinal circuitry controlling interlimb coordination (Danner et al., 2017). This model included the local V3 neurons that mediate mutual excitation between left and right rhythm generators (RGs). Here, our focus was on V3 neurons involved in ascending long propriospinal interactions (aLPNs). Using retrograde tracing, we revealed a subpopulation of lumbar V3 aLPNs with contralateral cervical projections. V3OFF mice, in which all V3 neurons were silenced, had a significantly reduced maximal locomotor speed, were unable to move using stable trot, gallop, or bound, and predominantly used a lateral-sequence walk. To reproduce this data and understand the functional roles of V3 aLPNs, we extended our previous model by incorporating diagonal V3 aLPNs mediating inputs from each lumbar RG to the contralateral cervical RG. The extended model reproduces our experimental results and suggests that locally projecting V3 neurons, mediating left-right interactions within lumbar and cervical cords, promote left-right synchronization necessary for gallop and bound, whereas the V3 aLPNs promote synchronization between diagonal fore and hind RGs necessary for trot. The model proposes the organization of spinal circuits available for future experimental testing.

Computational Modeling of Spinal Locomotor Circuitry in the Age of Molecular Genetics.

Neuronal circuits in the spinal cord are essential for the control of locomotion. They integrate supraspinal commands and afferent feedback signals to produce coordinated rhythmic muscle activations necessary for stable locomotion. For several decades, computational modeling has complemented experimental studies by providing a mechanistic rationale for experimental observations and by deriving experimentally testable predictions. This symbiotic relationship between experimental and computational approaches has resulted in numerous fundamental insights. With recent advances in molecular and genetic methods, it has become possible to manipulate specific constituent elements of the spinal circuitry and relate them to locomotor behavior. This has led to computational modeling studies investigating mechanisms at the level of genetically defined neuronal populations and their interactions. We review literature on the spinal locomotor circuitry from a computational perspective. By reviewing examples leading up to and in the age of molecular genetics, we demonstrate the importance of computational modeling and its interactions with experiments. Moving forward, neuromechanical models with neuronal circuitry modeled at the level of genetically defined neuronal populations will be required to further unravel the mechanisms by which neuronal interactions lead to locomotor behavior.

Neural Interactions in Developing Rhythmogenic Spinal Networks: Insights From Computational Modeling.

The mechanisms involved in generation of rhythmic locomotor activity in the mammalian spinal cord remain poorly understood. These mechanisms supposedly rely on both intrinsic properties of constituting neurons and interactions between them. A subset of Shox2 neurons was suggested to contribute to generation of spinal locomotor activity, but the possible cellular basis for rhythmic bursting in these neurons remains unknown. Ha and Dougherty (2018) recently revealed the presence of bidirectional electrical coupling between Shox2 neurons in neonatal spinal cords, which can be critically involved in neuronal synchronization and generation of populational bursting. Gap junctional connections found between functionally-related Shox2 interneurons decrease with age, possibly being replaced by increasing interactions through chemical synapses. Here, we developed a computational model of a heterogeneous population of neurons sparsely connected by electrical or/and chemical synapses and investigated the dependence of frequency of populational bursting on the type and strength of neuronal interconnections. The model proposes a mechanistic explanation that can account for the emergence of a synchronized rhythmic activity in the neuronal population and provides insights into the possible role of gap junctional coupling between Shox2 neurons in the spinal mechanisms for locomotor rhythm generation.

Spinal V3 Interneurons and Left-Right Coordination in Mammalian Locomotion.

Commissural interneurons (CINs) mediate interactions between rhythm-generating locomotor circuits located on each side of the spinal cord and are necessary for left-right limb coordination during locomotion. While glutamatergic V3 CINs have been implicated in left-right coordination, their functional connectivity remains elusive. Here, we addressed this issue by combining experimental and modeling approaches. We employed Sim1Cre/+; Ai32 mice, in which light-activated Channelrhodopsin-2 was selectively expressed in V3 interneurons. Fictive locomotor activity was evoked by NMDA and 5-HT in the isolated neonatal lumbar spinal cord. Flexor and extensor activities were recorded from left and right L2 and L5 ventral roots, respectively. Bilateral photoactivation of V3 interneurons increased the duration of extensor bursts resulting in a slowed down on-going rhythm. At high light intensities, extensor activity could become sustained. When light stimulation was shifted toward one side of the cord, the duration of extensor bursts still increased on both sides, but these changes were more pronounced on the contralateral side than on the ipsilateral side. Additional bursts appeared on the ipsilateral side not seen on the contralateral side. Further increase of the stimulation could suppress the contralateral oscillations by switching to a sustained extensor activity, while the ipsilateral rhythmic activity remained. To delineate the function of V3 interneurons and their connectivity, we developed a computational model of the spinal circuits consisting of two (left and right) rhythm generators (RGs) interacting via V0V, V0D, and V3 CINs. Both types of V0 CINs provided mutual inhibition between the left and right flexor RG centers and promoted left-right alternation. V3 CINs mediated mutual excitation between the left and right extensor RG centers. These interactions allowed the model to reproduce our current experimental data, while being consistent with previous data concerning the role of V0V and V0D CINs in securing left-right alternation and the changes in left-right coordination following their selective removal. We suggest that V3 CINs provide mutual excitation between the spinal neurons involved in the control of left and right extensor activity, which may promote left-right synchronization during locomotion.

Modulation of Respiratory System by Limb Muscle Afferents in Intact and Injured Spinal Cord.

Breathing constantly adapts to environmental, metabolic or behavioral changes by responding to different sensory information, including afferent feedback from muscles. Importantly, not just respiratory muscle feedback influences respiratory activity. Afferent sensory information from rhythmically moving limbs has also been shown to play an essential role in the breathing. The present review will discuss the neuronal mechanisms of respiratory modulation by activation of peripheral muscles that usually occurs during locomotion or exercise. An understanding of these mechanisms and finding the most effective approaches to regulate respiratory motor output by stimulation of limb muscles could be extremely beneficial for people with respiratory dysfunctions. Specific attention in the present review is given to the muscle stimulation to treat respiratory deficits following cervical spinal cord injury.

Computational modeling of brainstem circuits controlling locomotor frequency and gait.

A series of recent studies identified key structures in the mesencephalic locomotor region and the caudal brainstem of mice involved in the initiation and control of slow (exploratory) and fast (escape-type) locomotion and gait. However, the interactions of these brainstem centers with each other and with the spinal locomotor circuits are poorly understood. Previously we suggested that commissural and long propriospinal interneurons are the main targets for brainstem inputs adjusting gait (Danner et al., 2017). Here, by extending our previous model, we propose a connectome of the brainstem-spinal circuitry and suggest a mechanistic explanation of the operation of brainstem structures and their roles in controlling speed and gait. We suggest that brainstem control of locomotion is mediated by two pathways, one controlling locomotor speed via connections to rhythm generating circuits in the spinal cord and the other providing gait control by targeting commissural and long propriospinal interneurons.

State-dependent rhythmogenesis and frequency control in a half-center locomotor CPG.

The spinal locomotor central pattern generator (CPG) generates rhythmic activity with alternating flexion and extension phases. This rhythmic pattern is likely to result from inhibitory interactions between neural populations representing flexor and extensor half-centers. However, it is unclear whether the flexor-extensor CPG has a quasi-symmetric organization with both half-centers critically involved in rhythm generation, features an asymmetric organization with flexor-driven rhythmogenesis, or comprises a pair of intrinsically rhythmic half-centers. There are experimental data that support each of the above concepts but appear to be inconsistent with the others. In this theoretical/modeling study, we present and analyze a CPG model architecture that can operate in different regimes consistent with the above three concepts depending on conditions, which are defined by external excitatory drives to CPG half-centers. We show that control of frequency and phase durations within each regime depends on network dynamics, defined by the regime-dependent expression of the half-centers' intrinsic rhythmic capabilities and the operating phase transition mechanisms (escape vs. release). Our study suggests state dependency in locomotor CPG operation and proposes explanations for seemingly contradictory experimental data. NEW & NOTEWORTHY Our theoretical/modeling study focuses on the analysis of locomotor central pattern generators (CPGs) composed of conditionally bursting half-centers coupled with reciprocal inhibition and receiving independent external drives. We show that this CPG framework can operate in several regimes consistent with seemingly contradictory experimental data. In each regime, we study how intrinsic dynamics and phase-switching mechanisms control oscillation frequency and phase durations. Our results provide insights into the organization of spinal circuits controlling locomotion.

Computational modeling of spinal circuits controlling limb coordination and gaits in quadrupeds.

Interactions between cervical and lumbar spinal circuits are mediated by long propriospinal neurons (LPNs). Ablation of descending LPNs in mice disturbs left-right coordination at high speeds without affecting fore-hind alternation. We developed a computational model of spinal circuits consisting of four rhythm generators coupled by commissural interneurons (CINs), providing left-right interactions, and LPNs, mediating homolateral and diagonal interactions. The proposed CIN and diagonal LPN connections contribute to speed-dependent gait transition from walk, to trot, and then to gallop and bound; the homolateral LPN connections ensure fore-hind alternation in all gaits. The model reproduces speed-dependent gait expression in intact and genetically transformed mice and the disruption of hindlimb coordination following ablation of descending LPNs. Inputs to CINs and LPNs can affect interlimb coordination and change gait independent of speed. We suggest that these interneurons represent the main targets for supraspinal and sensory afferent signals adjusting gait.

Central control of interlimb coordination and speed-dependent gait expression in quadrupeds.

KEY POINTS: Quadrupeds express different gaits depending on speed of locomotion. Central pattern generators (one per limb) within the spinal cord generate locomotor oscillations and control limb movements. Neural interactions between these generators define interlimb coordination and gait. We present a computational model of spinal circuits representing four rhythm generators with left-right excitatory and inhibitory commissural and fore-hind inhibitory interactions within the cord. Increasing brainstem drive to all rhythm generators and excitatory commissural interneurons induces an increasing frequency of locomotor oscillations accompanied by speed-dependent gait changes from walk to trot and to gallop and bound. The model closely reproduces and suggests explanations for multiple experimental data, including speed-dependent gait transitions in intact mice and changes in gait expression in mutants lacking certain types of commissural interneurons. The model suggests the possible circuit organization in the spinal cord and proposes predictions that can be tested experimentally. ABSTRACT: As speed of locomotion is increasing, most quadrupeds, including mice, demonstrate sequential gait transitions from walk to trot and to gallop and bound. The neural mechanisms underlying these transitions are poorly understood. We propose that the speed-dependent expression of different gaits results from speed-dependent changes in the interactions between spinal circuits controlling different limbs and interlimb coordination. As a result, the expression of each gait depends on (1) left-right interactions within the spinal cord mediated by different commissural interneurons (CINs), (2) fore-hind interactions on each side of the spinal cord and (3) brainstem drives to rhythm-generating circuits and CIN pathways. We developed a computational model of spinal circuits consisting of four rhythm generators (RGs) with bilateral left-right interactions mediated by V0 CINs (V0D and V0V sub-types) providing left-right alternation, and conditional V3 CINs promoting left-right synchronization. Fore and hind RGs mutually inhibited each other. We demonstrate that linearly increasing excitatory drives to the RGs and V3 CINs can produce a progressive increase in the locomotor speed accompanied by sequential changes of gaits from walk to trot and to gallop and bound. The model closely reproduces and suggests explanations for the speed-dependent gait expression observed in vivo in intact mice and in mutants lacking V0V or all V0 CINs. Specifically, trot is not expressed after removal of V0V CINs, and only bound is expressed after removal of all V0 CINs. The model provides important insights into the organization of spinal circuits and neural control of locomotion.

Organization of flexor-extensor interactions in the mammalian spinal cord: insights from computational modelling.

KEY POINTS: Alternation of flexor and extensor activity in the mammalian spinal cord is mediated by two classes of genetically identified inhibitory interneurons, V1 and V2b. The V1 interneurons are essential for high-speed locomotor activity. They secure flexor-extensor alternations in the intact cord but lose this function after hemisection, suggesting that they are activated by inputs from the contralateral side of the cord. The V2b interneurons are involved in flexor-extensor alternations in both intact cord and hemicords. We used a computational model of the spinal network, simulating the left and right rhythm-generating circuits interacting via several commissural pathways, and extended this model by incorporating V1 and V2b neuron populations involved in flexor-extensor interactions on each cord side. The model reproduces multiple experimental data on selective silencing and activation of V1 and/or V2b neurons and proposes the organization of their connectivity providing flexor-extensor alternation in the spinal cord. ABSTRACT: Alternating flexor and extensor activity represents the fundamental property underlying many motor behaviours including locomotion. During locomotion this alternation appears to arise in rhythm-generating circuits and transpires at all levels of the spinal cord including motoneurons. Recent studies in vitro and in vivo have shown that flexor-extensor alternation during locomotion involves two classes of genetically identified, inhibitory interneurons: V1 and V2b. Particularly, in the isolated mouse spinal cord, abrogation of neurotransmission derived by both V1 and V2b interneurons resulted in flexor-extensor synchronization, whereas selective inactivation of only one of these neuron types did not abolish flexor-extensor alternation. After hemisection, inactivation of only V2b interneurons led to the flexor-extensor synchronization, while inactivation of V1 interneurons did not affect flexor-extensor alternation. Moreover, optogenetic activation of V2b interneurons suppressed extensor-related activity, while similar activation of V1 interneurons suppressed both flexor and extensor oscillations. Here, we address these issues using the previously published computational model of spinal circuitry simulating bilateral interactions between left and right rhythm-generating circuits. In the present study, we incorporate V1 and V2b neuron populations on both sides of the cord to make them critically involved in flexor-extensor interactions. The model reproduces multiple experimental data on the effects of hemisection and selective silencing or activation of V1 and V2b neurons and suggests connectivity profiles of these neurons and their specific roles in left-right (V1) and flexor-extensor (both V2b and V1) interactions in the spinal cord that can be tested experimentally.

Organization of the Mammalian Locomotor CPG: Review of Computational Model and Circuit Architectures Based on Genetically Identified Spinal Interneurons(1,2,3).

The organization of neural circuits that form the locomotor central pattern generator (CPG) and provide flexor-extensor and left-right coordination of neuronal activity remains largely unknown. However, significant progress has been made in the molecular/genetic identification of several types of spinal interneurons, including V0 (V0D and V0V subtypes), V1, V2a, V2b, V3, and Shox2, among others. The possible functional roles of these interneurons can be suggested from changes in the locomotor pattern generated in mutant mice lacking particular neuron types. Computational modeling of spinal circuits may complement these studies by bringing together data from different experimental studies and proposing the possible connectivity of these interneurons that may define rhythm generation, flexor-extensor interactions on each side of the cord, and commissural interactions between left and right circuits. This review focuses on the analysis of potential architectures of spinal circuits that can reproduce recent results and suggest common explanations for a series of experimental data on genetically identified spinal interneurons, including the consequences of their genetic ablation, and provides important insights into the organization of the spinal CPG and neural control of locomotion.

Organization of left-right coordination of neuronal activity in the mammalian spinal cord: Insights from computational modelling.

KEY POINTS: Coordination of neuronal activity between left and right sides of the mammalian spinal cord is provided by several sets of commissural interneurons (CINs) whose axons cross the midline. Genetically identified inhibitory V0D and excitatory V0V CINs and ipsilaterally projecting excitatory V2a interneurons were shown to secure left-right alternation at different locomotor speeds. We have developed computational models of neuronal circuits in the spinal cord that include left and right rhythm-generating centres interacting bilaterally via three parallel pathways mediated by V0D , V2a-V0V and V3 neuron populations. The models reproduce the experimentally observed speed-dependent left-right coordination in normal mice and the changes in coordination seen in mutants lacking specific neuron classes. The models propose an explanation for several experimental results and provide insights into the organization of the spinal locomotor network and parallel CIN pathways involved in gait control at different locomotor speeds. ABSTRACT: Different locomotor gaits in mammals, such as walking or galloping, are produced by coordinated activity in neuronal circuits in the spinal cord. Coordination of neuronal activity between left and right sides of the cord is provided by commissural interneurons (CINs), whose axons cross the midline. In this study, we construct and analyse two computational models of spinal locomotor circuits consisting of left and right rhythm generators interacting bilaterally via several neuronal pathways mediated by different CINs. The CIN populations incorporated in the models include the genetically identified inhibitory (V0D ) and excitatory (V0V ) subtypes of V0 CINs and excitatory V3 CINs. The model also includes the ipsilaterally projecting excitatory V2a interneurons mediating excitatory drive to the V0V CINs. The proposed network architectures and CIN connectivity allow the models to closely reproduce and suggest mechanistic explanations for several experimental observations. These phenomena include: different speed-dependent contributions of V0D and V0V CINs and V2a interneurons to left-right alternation of neural activity, switching gaits between the left-right alternating walking-like activity and the left-right synchronous hopping-like pattern in mutants lacking specific neuron classes, and speed-dependent asymmetric changes of flexor and extensor phase durations. The models provide insights into the architecture of spinal network and the organization of parallel inhibitory and excitatory CIN pathways and suggest explanations for how these pathways maintain alternating and synchronous gaits at different locomotor speeds. The models propose testable predictions about the neural organization and operation of mammalian locomotor circuits.

A closed-loop model of the respiratory system: focus on hypercapnia and active expiration.

Breathing is a vital process providing the exchange of gases between the lungs and atmosphere. During quiet breathing, pumping air from the lungs is mostly performed by contraction of the diaphragm during inspiration, and muscle contraction during expiration does not play a significant role in ventilation. In contrast, during intense exercise or severe hypercapnia forced or active expiration occurs in which the abdominal "expiratory" muscles become actively involved in breathing. The mechanisms of this transition remain unknown. To study these mechanisms, we developed a computational model of the closed-loop respiratory system that describes the brainstem respiratory network controlling the pulmonary subsystem representing lung biomechanics and gas (O2 and CO2) exchange and transport. The lung subsystem provides two types of feedback to the neural subsystem: a mechanical one from pulmonary stretch receptors and a chemical one from central chemoreceptors. The neural component of the model simulates the respiratory network that includes several interacting respiratory neuron types within the Bötzinger and pre-Bötzinger complexes, as well as the retrotrapezoid nucleus/parafacial respiratory group (RTN/pFRG) representing the central chemoreception module targeted by chemical feedback. The RTN/pFRG compartment contains an independent neural generator that is activated at an increased CO2 level and controls the abdominal motor output. The lung volume is controlled by two pumps, a major one driven by the diaphragm and an additional one activated by abdominal muscles and involved in active expiration. The model represents the first attempt to model the transition from quiet breathing to breathing with active expiration. The model suggests that the closed-loop respiratory control system switches to active expiration via a quantal acceleration of expiratory activity, when increases in breathing rate and phrenic amplitude no longer provide sufficient ventilation. The model can be used for simulation of closed-loop control of breathing under different conditions including respiratory disorders.

Rhythmic bursting in the pre-Bötzinger complex: mechanisms and models.

The pre-Bötzinger complex (pre-BötC), a neural structure involved in respiratory rhythm generation, can generate rhythmic bursting activity in vitro that persists after blockade of synaptic inhibition. Experimental studies have identified two mechanisms potentially involved in this activity: one based on the persistent sodium current (INaP) and the other involving calcium (ICa) and/or calcium-activated nonspecific cation (ICAN) currents. In this modeling study, we investigated bursting generated in single neurons and excitatory neural populations with randomly distributed conductances of INaP and ICa. We analyzed the possible roles of these currents, the Na(+)/K(+) pump, synaptic mechanisms, and network interactions in rhythmic bursting generated under different conditions. We show that a population of synaptically coupled excitatory neurons with randomly distributed INaP- and/or ICAN-mediated burst generating mechanisms can operate in different oscillatory regimes with bursting dependent on either current or independent of both. The existence of multiple oscillatory regimes and their state dependence may explain rhythmic activities observed in the pre-BötC under different conditions.

Effects of glycinergic inhibition failure on respiratory rhythm and pattern generation.

Inhibitory interactions between neurons of the respiratory network are involved in rhythm generation and pattern formation. Using a computational model of brainstem respiratory networks, we investigated the possible effects of suppressing glycinergic inhibition on the activity of different respiratory neuron types. Our study revealed that progressive suppression of glycinergic inhibition affected all neurons of the network and disturbed neural circuits involved in termination of inspiration. Causal was a dysfunction of postinspiratory inhibition targeting inspiratory neurons, which often led to irregular preterm reactivation of these neurons, producing double or multiple short-duration inspiratory bursts. An increasing blockade of glycinergic inhibition led to apneustic inspiratory activity. Similar disturbances of glycinergic inhibition also occur during hypoxia. A clear difference in prolonged hypoxia, however, is that the rhythm terminates in expiratory apnea. The critical function of glycinergic inhibition for normal respiratory rhythm generation and the consequences of its reduction, including in pathological conditions, are discussed.