Eupnea, Tachypnea, and Autoresuscitation in a Closed‐Loop Respiratory Control Model

How sensory information influences the dynamics of rhythm generation varies across systems, and general principles for understanding this aspect of motor control are lacking. Determining the origin of respiratory rhythm generation is challenging because the mechanisms in a central circuit considered in isolation may be different from those in the intact organism. We analyze a closed-loop respiratory control model incorporating a central pattern generator (CPG), the Butera-Rinzel-Smith (BRS) model, together with lung mechanics, oxygen handling, and chemosensory components. We show that 1) embedding the BRS model neuron in a control loop creates a bistable system; 2) although closed-loop and open-loop (isolated) CPG systems both support eupnea-like bursting activity, they do so via distinct mechanisms; 3) chemosensory feedback in the closed loop improves robustness to variable metabolic demand; 4) the BRS model conductances provide an autoresuscitation mechanism for recovery from transient interruption of chemosensory feedback; and 5) the in vitro brain stem CPG slice responds to hypoxia with transient bursting that is qualitatively similar to in silico autoresuscitation. Bistability of bursting and tonic spiking in the closed-loop system corresponds to coexistence of eupnea-like breathing, with normal minute ventilation and blood oxygen level and a tachypnea-like state, with pathologically reduced minute ventilation and critically low blood oxygen. Disruption of the normal breathing rhythm, through either imposition of hypoxia or interruption of chemosensory feedback, can push the system from the eupneic state into the tachypneic state. We use geometric singular perturbation theory to analyze the system dynamics at the boundary separating eupnea-like and tachypnea-like outcomes.NEW & NOTEWORTHY A common challenge facing rhythmic biological processes is the adaptive regulation of central pattern generator (CPG) activity in response to sensory feedback. We apply dynamical systems tools to understand several properties of a closed-loop respiratory control model, including the coexistence of normal and pathological breathing, robustness to changes in metabolic demand, spontaneous autoresuscitation in response to hypoxia, and the distinct mechanisms that underlie rhythmogenesis in the intact control circuit vs. the isolated, open-loop CPG.

[1]  David Terman,et al.  Mathematical foundations of neuroscience , 2010 .

[2]  U. Bässler On the definition of central pattern generator and its sensory control , 1986, Biological Cybernetics.

[3]  J. Widdicombe Pulmonary and respiratory tract receptors. , 1982, The Journal of experimental biology.

[4]  Perry L. Miller,et al.  Twenty years of ModelDB and beyond: building essential modeling tools for the future of neuroscience , 2016, Journal of Computational Neuroscience.

[5]  M. Khoo,et al.  A model-based evaluation of the single-breath CO2 ventilatory response test. , 1990, Journal of applied physiology.

[6]  Christopher Jones,et al.  Geometric singular perturbation theory , 1995 .

[7]  J. C. Smith,et al.  Models of respiratory rhythm generation in the pre-Bötzinger complex. II. Populations Of coupled pacemaker neurons. , 1999, Journal of neurophysiology.

[8]  J. Ramirez,et al.  Respiratory rhythm generation in mammals: synaptic and membrane properties. , 1997, Respiration physiology.

[9]  F. Verhulst,et al.  Averaging Methods in Nonlinear Dynamical Systems , 1985 .

[10]  John Rinzel,et al.  A Formal Classification of Bursting Mechanisms in Excitable Systems , 1987 .

[11]  Peter J. Thomas,et al.  Spontaneous autoresuscitation in a model of respiratory control , 2012, 2012 Annual International Conference of the IEEE Engineering in Medicine and Biology Society.

[12]  Alona Ben-Tal,et al.  A model for control of breathing in mammals: coupling neural dynamics to peripheral gas exchange and transport. , 2008, Journal of theoretical biology.

[13]  Alona Ben-Tal,et al.  Integrative approaches for modeling regulation and function of the respiratory system , 2013, Wiley interdisciplinary reviews. Systems biology and medicine.

[14]  M. Alexander,et al.  Principles of Neural Science , 1981 .

[15]  Christopher G. Wilson,et al.  Functional Imaging, Spatial Reconstruction, and Biophysical Analysis of a Respiratory Motor Circuit Isolated In Vitro , 2008, The Journal of Neuroscience.

[16]  Fred S. Grodins,et al.  Control Theory and Biological Systems , 1963 .

[17]  Neil S. Cherniack,et al.  Introduction of respiratory pattern generators into models of respiratory control , 2005, Respiratory Physiology & Neurobiology.

[18]  John Kattwinkel,et al.  Clinical Associations with Immature Breathing in Preterm Infants , 2016, Pediatric Research.

[19]  Alessandra Conversi,et al.  Comparative Analysis , 2009, Encyclopedia of Database Systems.

[20]  Jonathan E. Rubin,et al.  Synchronized Activity and Loss of Synchrony Among Heterogeneous Conditional Oscillators , 2002, SIAM J. Appl. Dyn. Syst..

[21]  Jonathan E. Rubin,et al.  Averaging, Folded Singularities, and Torus Canards: Explaining Transitions between Bursting and Spiking in a Coupled Neuron Model , 2015, SIAM J. Appl. Dyn. Syst..

[22]  Bard Ermentrout,et al.  Simulating, analyzing, and animating dynamical systems - a guide to XPPAUT for researchers and students , 2002, Software, environments, tools.

[23]  Richard Bertram,et al.  Multi-timescale systems and fast-slow analysis. , 2017, Mathematical biosciences.

[24]  Neil Fenichel Geometric singular perturbation theory for ordinary differential equations , 1979 .

[25]  P. Bressloff,et al.  Bursting: The genesis of rhythm in the nervous system , 2005 .

[26]  Eugene M. Izhikevich,et al.  Neural excitability, Spiking and bursting , 2000, Int. J. Bifurc. Chaos.

[27]  G. Ermentrout,et al.  Multiple rhythmic states in a model of the respiratory central pattern generator. , 2009, Journal of neurophysiology.

[28]  Robert Clewley,et al.  A Computational Tool for the Reduction of Nonlinear ODE Systems Possessing Multiple Scales , 2005, Multiscale Model. Simul..

[29]  L. Perko Differential Equations and Dynamical Systems , 1991 .

[30]  J. Feldman,et al.  Sodium and Calcium Current-Mediated Pacemaker Neurons and Respiratory Rhythm Generation , 2005, The Journal of Neuroscience.

[31]  J. C. Smith,et al.  Pre-Bötzinger complex: a brainstem region that may generate respiratory rhythm in mammals. , 1991, Science.

[32]  P. Gaspard Measurement of the instability rate of a far-from-equilibrium steady state at an infinite period bifurcation , 1990 .

[33]  Yangyang Wang,et al.  Multiple timescale mixed bursting dynamics in a respiratory neuron model , 2016, Journal of Computational Neuroscience.

[34]  Alla Borisyuk,et al.  UNDERSTANDING NEURONAL DYNAMICS BY GEOMETRICAL DISSECTION OF MINIMAL MODELS , 2005 .

[35]  Jeffrey C. Smith,et al.  Neuronal pacemaker for breathing visualized in vitro , 1999, Nature.

[36]  Eugene M. Izhikevich,et al.  Dynamical Systems in Neuroscience: The Geometry of Excitability and Bursting , 2006 .

[37]  M. Pembrey,et al.  The estimation of fetal haemoglobin in healthy adults by radioimmunoassay , 1983, British journal of haematology.

[38]  Sergey N Markin,et al.  Spatial organization and state-dependent mechanisms for respiratory rhythm and pattern generation. , 2007, Progress in brain research.

[39]  G. Ermentrout,et al.  Analysis of neural excitability and oscillations , 1989 .

[40]  H. Coleridge,et al.  Pulmonary reflexes: neural mechanisms of pulmonary defense. , 1994, Annual review of physiology.

[41]  E. Herlenius,et al.  Il-1β and prostaglandin E2 attenuate the hypercapnic as well as the hypoxic respiratory response via prostaglandin E receptor type 3 in neonatal mice. , 2014, Journal of applied physiology.

[42]  S. Wiggins Normally Hyperbolic Invariant Manifolds in Dynamical Systems , 1994 .

[43]  A. Paintal,et al.  Vagal sensory receptors and their reflex effects. , 1973, Physiological reviews.

[44]  J. Rubin,et al.  Geometric Singular Perturbation Analysis of Neuronal Dynamics , 2002 .

[45]  Manu S. Madhav,et al.  Feedback control as a framework for understanding tradeoffs in biology. , 2014, Integrative and comparative biology.

[46]  N. Prabhakar,et al.  Comparative analysis of neonatal and adult rat carotid body responses to chronic intermittent hypoxia. , 2008, Journal of applied physiology.

[47]  J C Smith,et al.  Spatial and functional architecture of the mammalian brain stem respiratory network: a hierarchy of three oscillatory mechanisms. , 2007, Journal of neurophysiology.

[48]  Limei Cheng,et al.  An integrative model of respiratory and cardiovascular control in sleep-disordered breathing , 2010, Respiratory Physiology & Neurobiology.

[49]  R J Full,et al.  Templates and anchors: neuromechanical hypotheses of legged locomotion on land. , 1999, The Journal of experimental biology.

[50]  D. Richter,et al.  L‐type Ca2+ channels in inspiratory neurones of mice and their modulation by hypoxia , 1998, The Journal of physiology.

[51]  D. Richter,et al.  Intracellular signalling pathways modulate KATP channels in inspiratory brainstem neurones and their hypoxic activation: involvement of metabotropic receptors, G-proteins and cytoskeleton , 2000, Brain Research.

[52]  R. Bertram,et al.  Topological and phenomenological classification of bursting oscillations , 1995 .

[53]  W. Guntheroth Cheyne-Stokes respiration: hypoxia plus a deep breath that interrupts hypoxic drive, initiating cyclic breathing. , 2011, Medical hypotheses.

[54]  Jonathan E. Rubin,et al.  Interactions of persistent sodium and calcium-activated nonspecific cationic currents yield dynamically distinct bursting regimes in a model of respiratory neurons , 2011, Journal of Computational Neuroscience.

[55]  Bruce G Lindsey,et al.  Computational models and emergent properties of respiratory neural networks. , 2012, Comprehensive Physiology.

[56]  Jeffery T. Erickson,et al.  Autoresuscitation responses to hypoxia-induced apnea are delayed in newborn 5-HT-deficient Pet-1 homozygous mice. , 2009, Journal of applied physiology.

[57]  A. Rudolph,et al.  Perinatal oxygen delivery and cardiac function. , 1985, Advances in pediatrics.

[58]  F. Eldridge,et al.  The North Carolina Respiratory Model , 1996 .

[59]  A. J. Bart,et al.  Mathematical analysis and digital simulation of the respiratory control system. , 1967, Journal of applied physiology.

[60]  Alla Borisyuk,et al.  The Dynamic Range of Bursting in a Model Respiratory Pacemaker Network , 2005, SIAM J. Appl. Dyn. Syst..

[61]  D. Richter,et al.  Oscillations and hypoxic changes of mitochondrial variables in neurons of the brainstem respiratory centre of mice , 2001, The Journal of physiology.

[62]  H. Koepchen,et al.  A Systems View of Respiratory Regulation , 1996 .

[63]  Hillel J. Chiel,et al.  Robustness, flexibility, and sensitivity in a multifunctional motor control model , 2016, Biological Cybernetics.

[64]  D. Richter,et al.  Hypoxia activates ATP‐dependent potassium channels in inspiratory neurones of neonatal mice , 1998, The Journal of physiology.

[65]  Nahum Shimkin,et al.  Nonlinear Control Systems , 2008 .

[66]  Mathieu Desroches,et al.  A showcase of torus canards in neuronal bursters , 2012, Journal of mathematical neuroscience.

[67]  Kendall F. Morris,et al.  A Joint Computational Respiratory Neural Network-Biomechanical Model for Breathing and Airway Defensive Behaviors , 2012, Front. Physio..

[68]  John Kattwinkel,et al.  Quantification of periodic breathing in premature infants , 2015, Physiological measurement.

[69]  J C Smith,et al.  Respiratory rhythm generation in neonatal and adult mammals: the hybrid pacemaker-network model. , 2000, Respiration physiology.

[70]  E Roth,et al.  A comparative approach to closed-loop computation , 2014, Current Opinion in Neurobiology.

[71]  N. Prabhakar,et al.  Altered carotid body function by intermittent hypoxia in neonates and adults: Relevance to recurrent apneas , 2007, Respiratory Physiology & Neurobiology.

[72]  John W. Clark,et al.  Dissection and reduction of a modeled bursting neuron , 1996, Journal of Computational Neuroscience.

[73]  Nancy J Kopell,et al.  New dynamics in cerebellar Purkinje cells: torus canards. , 2008, Physical review letters.

[74]  J. Neubauer,et al.  Heme oxygenase is necessary for the excitatory response of cultured neonatal rat rostral ventrolateral medulla neurons to hypoxia. , 2009, American journal of physiology. Regulatory, integrative and comparative physiology.

[75]  Ilya A. Rybak,et al.  A Closed-Loop Model of the Respiratory System: Focus on Hypercapnia and Active Expiration , 2014, PloS one.

[76]  J. Feldman,et al.  Understanding the rhythm of breathing: so near, yet so far. , 2013, Annual review of physiology.

[77]  John B. West,et al.  Respiratory Physiology - the Essentials , 1979 .

[78]  R. W. Jones,et al.  Respiratory responses to CO2 inhalation; a theoretical study of a nonlinear biological regulator. , 1954, Journal of applied physiology.

[79]  Michael P. Hlastala,et al.  Physiology of respiration , 1996 .

[80]  J. C. Smith,et al.  Models of respiratory rhythm generation in the pre-Bötzinger complex. I. Bursting pacemaker neurons. , 1999, Journal of neurophysiology.

[81]  Mark A Kramer,et al.  An elementary model of torus canards. , 2011, Chaos.

[82]  R. Stornetta,et al.  GABAergic pump cells of solitary tract nucleus innervate retrotrapezoid nucleus chemoreceptors , 2007, Journal of neurophysiology.

[83]  R. Martin,et al.  Physiologic basis for intermittent hypoxic episodes in preterm infants. , 2012, Advances in experimental medicine and biology.

[84]  S. Givan,et al.  Intermittent severe hypoxia induces plasticity within serotonergic and catecholaminergic neurons in the neonatal rat ventrolateral medulla. , 2016, Journal of applied physiology.