Transient Receptor Potential Channels TRPM4 and TRPC3 Critically Contribute to Respiratory Motor Pattern Formation but not Rhythmogenesis in Rodent Brainstem Circuits

Interdependence of cellular and network properties in respiratory rhythm generation

How breathing is generated by the preBötzinger complex (preBötC) remains divided between two ideological frameworks, and a persistent sodium current (INaP) lies at the heart of this debate. Although INaP is widely expressed, the pacemaker hypothesis considers it essential because it endows a small subset of neurons with intrinsic bursting or "pacemaker" activity. In contrast, burstlet theory considers INaP dispensable because rhythm emerges from "preinspiratory" spiking activity driven by feed-forward network interactions. Using computational modeling, we find that small changes in spike shape can dissociate INaP from intrinsic bursting. Consistent with many experimental benchmarks, conditional effects on spike shape during simulated changes in oxygenation, development, extracellular potassium, and temperature alter the prevalence of intrinsic bursting and preinspiratory spiking without altering the role of INaP. Our results support a unifying hypothesis where INaP and excitatory network interactions, but not intrinsic bursting or preinspiratory spiking, are critical interdependent features of preBötC rhythmogenesis.

Inspiratory and sigh breathing rhythms depend on distinct cellular signalling mechanisms in the preBötzinger complex

Breathing behaviour involves the generation of normal breaths (eupnoea) on a timescale of seconds and sigh breaths on the order of minutes. Both rhythms emerge in tandem from a single brainstem site, but whether and how a single cell population can generate two disparate rhythms remains unclear. We posit that recurrent synaptic excitation in concert with synaptic depression and cellular refractoriness gives rise to the eupnoea rhythm, whereas an intracellular calcium oscillation that is slower by orders of magnitude gives rise to the sigh rhythm. A mathematical model capturing these dynamics simultaneously generates eupnoea and sigh rhythms with disparate frequencies, which can be separately regulated by physiological parameters. We experimentally validated key model predictions regarding intracellular calcium signalling. All vertebrate brains feature a network oscillator that drives the breathing pump for regular respiration. However, in air-breathing mammals with compliant lungs susceptible to collapse, the breathing rhythmogenic network may have refashioned ubiquitous intracellular signalling systems to produce a second slower rhythm (for sighs) that prevents atelectasis without impeding eupnoea. KEY POINTS: A simplified activity-based model of the preBötC generates inspiratory and sigh rhythms from a single neuron population. Inspiration is attributable to a canonical excitatory network oscillator mechanism. Sigh emerges from intracellular calcium signalling. The model predicts that perturbations of calcium uptake and release across the endoplasmic reticulum counterintuitively accelerate and decelerate sigh rhythmicity, respectively, which was experimentally validated. Vertebrate evolution may have adapted existing intracellular signalling mechanisms to produce slow oscillations needed to optimize pulmonary function in mammals.

Interdependence of cellular and network properties in respiratory rhythmogenesis

How breathing is generated by the preBötzinger Complex (preBötC) remains divided between two ideological frameworks, and the persistent sodium current (INaP) lies at the heart of this debate. Although INaP is widely expressed, the pacemaker hypothesis considers it essential because it endows a small subset of neurons with intrinsic bursting or "pacemaker" activity. In contrast, burstlet theory considers INaP dispensable because rhythm emerges from "pre-inspiratory" spiking activity driven by feed-forward network interactions. Using computational modeling, we discover that changes in spike shape can dissociate INaP from intrinsic bursting. Consistent with many experimental benchmarks, conditional effects on spike shape during simulated changes in oxygenation, development, extracellular potassium, and temperature alter the prevalence of intrinsic bursting and pre-inspiratory spiking without altering the role of INaP. Our results support a unifying hypothesis where INaP and excitatory network interactions, but not intrinsic bursting or pre-inspiratory spiking, are critical interdependent features of preBötC rhythmogenesis.

Role of NaV1.6-mediated persistent sodium current and bursting-pacemaker properties in breathing rhythm generation

Inspiration is the inexorable active phase of breathing. The brainstem pre-Bötzinger complex (preBötC) gives rise to inspiratory neural rhythm, but its underlying cellular and ionic bases remain unclear. The long-standing "pacemaker hypothesis" posits that the persistent Na+ current (INaP) that gives rise to bursting-pacemaker properties in preBötC interneurons is essential for rhythmogenesis. We tested the pacemaker hypothesis by conditionally knocking out and knocking down the Scn8a (Nav1.6 [voltage-gated sodium channel 1.6]) gene in core rhythmogenic preBötC neurons. Deleting Scn8a substantially decreases the INaP and abolishes bursting-pacemaker activity, which slows inspiratory rhythm in vitro and negatively impacts the postnatal development of ventilation. Diminishing Scn8a via genetic interference has no impact on breathing in adult mice. We argue that the Scn8a-mediated INaP is not obligatory but that it influences the development and rhythmic function of the preBötC. The ubiquity of the INaP in respiratory brainstem interneurons could underlie breathing-related behaviors such as neonatal phonation or rhythmogenesis in different physiological conditions.

Modulation of inspiratory burst duration and frequency by bombesin in vitro

Mammalian respiratory rhythm-generating circuits in the brainstem are subject to neuromodulation by multiple peptidergic afferent inputs controlling circuit behavior and outputs. Although functionally important, actions of neuropeptide modulators have not been fully characterized. We analyzed at cellular and circuit levels two inspiratory patterns intrinsically generated by the preBötzinger complex (preBötC) and their modulation by the neuropeptides bombesin and substance P (SP) in neonatal rat medullary slices in vitro. We found that, in recordings of hypoglossal nerve and preBötC neuron inspiratory activity, some inspiratory bursts occurring spontaneously under basal conditions have a biphasic shape with longer duration than normal inspiratory bursts and occur at a lower frequency. This biphasic burst pattern has been proposed to represent inspiratory activity underling periodic sighs. Bath-applied bombesin or SP decreased the period and increased the duration of both normal inspiratory and biphasic bursts and their underlying synaptic drives. The ratio of the biphasic long-duration burst period to the normal inspiratory burst period and the ratio of their burst durations remained the same before and after peptidergic modulation. Bombesin increased the frequency of the inspiratory rhythm in a Ca²⁺-independent manner and the frequency of long-duration bursts in a Ca²⁺-dependent manner. This finding suggests that period and burst duration coupling are due to intrinsic mechanisms controlling simultaneously timing and burst termination within the inspiratory rhythm-generating network. We propose a model in which signaling cascades activated by bombesin and SP modulate mechanisms controlling inspiratory burst frequency and duration to coordinate preBötC circuit behavioral outputs.

Predictions and experimental tests of a new biophysical model of the mammalian respiratory oscillator

Previously our computational modeling studies (Phillips et al., 2019) proposed that neuronal persistent sodium current (INaP) and calcium-activated non-selective cation current (ICAN) are key biophysical factors that, respectively, generate inspiratory rhythm and burst pattern in the mammalian preBötzinger complex (preBötC) respiratory oscillator isolated in vitro. Here, we experimentally tested and confirmed three predictions of the model from new simulations concerning the roles of INaP and ICAN: (1) INaP and ICAN blockade have opposite effects on the relationship between network excitability and preBötC rhythmic activity; (2) INaP is essential for preBötC rhythmogenesis; and (3) ICAN is essential for generating the amplitude of rhythmic output but not rhythm generation. These predictions were confirmed via optogenetic manipulations of preBötC network excitability during graded INaP or ICAN blockade by pharmacological manipulations in slices in vitro containing the rhythmically active preBötC from the medulla oblongata of neonatal mice. Our results support and advance the hypothesis that INaP and ICAN mechanistically underlie rhythm and inspiratory burst pattern generation, respectively, in the isolated preBötC.

Inspiratory rhythm generation is stabilized by Ih

Cellular and network properties must be capable of generating rhythmic activity that is both flexible and stable. This is particularly important for breathing, a rhythmic behavior that dynamically adapts to environmental, behavioral, and metabolic changes from the first to the last breath. The pre-Bötzinger complex (preBötC), located within the ventral medulla, is responsible for producing rhythmic inspiration. Its cellular properties must be tunable, flexible as well as stabilizing. Here, we explore the role of the hyperpolarization-activated, nonselective cation current (I_h) for stabilizing PreBötC activity during opioid exposure and reduced excitatory synaptic transmission. Introducing I_h into an in silico preBötC network predicts that loss of this depolarizing current should significantly slow the inspiratory rhythm. By contrast, in vitro and in vivo experiments revealed that the loss of I_h minimally affected breathing frequency, but destabilized rhythmogenesis through the generation of incompletely synchronized bursts (burstlets). Associated with the loss of I_h was an increased susceptibility of breathing to opioid-induced respiratory depression or weakened excitatory synaptic interactions, a paradoxical depolarization at the cellular level, and the suppression of tonic spiking. Tonic spiking activity is generated by nonrhythmic excitatory and inhibitory preBötC neurons, of which a large percentage express I_h. Together, our results suggest that I_h is important for maintaining tonic spiking, stabilizing inspiratory rhythmogenesis, and protecting breathing against perturbations or changes in network state.NEW & NOTEWORTHY The I_h current plays multiple roles within the preBötC. This current is important for promoting intrinsic tonic spiking activity in excitatory and inhibitory neurons and for preserving rhythmic function during conditions that dampen network excitability, such as in the context of opioid-induced respiratory depression. We therefore propose that the I_h current expands the dynamic range of rhythmogenesis, buffers the preBötC against network perturbations, and stabilizes rhythmogenesis by preventing the generation of unsynchronized bursts.

A glibenclamide-sensitive TRPM4-mediated component of CA1 excitatory postsynaptic potentials appears in experimental autoimmune encephalomyelitis

The transient receptor potential melastatin 4 (TRPM4) channel contributes to disease severity in the murine experimental autoimmune encephalomyelitis (EAE) model of multiple sclerosis and to neuronal cell death in models of excitotoxicity and traumatic brain injury. As TRPM4 is activated by intracellular calcium and conducts monovalent cations, we hypothesized that TRPM4 may contribute to and boost excitatory synaptic transmission in CA1 pyramidal neurons of the hippocampus. Using single-spine calcium imaging and electrophysiology, we found no effect of the TRPM4 antagonists 9-phenanthrol and glibenclamide on synaptic transmission in hippocampal slices from healthy mice. In contrast, glibenclamide but not 9-phenanthrol reduced excitatory synaptic potentials in slices from EAE mice, an effect that was absent in slices from EAE mice lacking TRPM4. We conclude that TRPM4 plays little role in basal hippocampal synaptic transmission, but a glibenclamide-sensitive TRPM4-mediated contribution to excitatory postsynaptic responses is upregulated at the acute phase of EAE.

Putting the theory into 'burstlet theory' with a biophysical model of burstlets and bursts in the respiratory preBötzinger complex

Inspiratory breathing rhythms arise from synchronized neuronal activity in a bilaterally distributed brainstem structure known as the preBötzinger complex (preBötC). In in vitro slice preparations containing the preBötC, extracellular potassium must be elevated above physiological levels (to 7-9 mM) to observe regular rhythmic respiratory motor output in the hypoglossal nerve to which the preBötC projects. Reexamination of how extracellular K+ affects preBötC neuronal activity has revealed that low-amplitude oscillations persist at physiological levels. These oscillatory events are subthreshold from the standpoint of transmission to motor output and are dubbed burstlets. Burstlets arise from synchronized neural activity in a rhythmogenic neuronal subpopulation within the preBötC that in some instances may fail to recruit the larger network events, or bursts, required to generate motor output. The fraction of subthreshold preBötC oscillatory events (burstlet fraction) decreases sigmoidally with increasing extracellular potassium. These observations underlie the burstlet theory of respiratory rhythm generation. Experimental and computational studies have suggested that recruitment of the non-rhythmogenic component of the preBötC population requires intracellular Ca2+ dynamics and activation of a calcium-activated nonselective cationic current. In this computational study, we show how intracellular calcium dynamics driven by synaptically triggered Ca2+ influx as well as Ca2+ release/uptake by the endoplasmic reticulum in conjunction with a calcium-activated nonselective cationic current can reproduce and offer an explanation for many of the key properties associated with the burstlet theory of respiratory rhythm generation. Altogether, our modeling work provides a mechanistic basis that can unify a wide range of experimental findings on rhythm generation and motor output recruitment in the preBötC.

Advancing respiratory-cardiovascular physiology with the working heart-brainstem preparation over 25 years

Twenty‐five years ago, a new physiological preparation called the working heart–brainstem preparation (WHBP) was introduced with the claim it would provide a new platform allowing studies not possible before in cardiovascular, neuroendocrine, autonomic and respiratory research. Herein, we review some of the progress made with the WHBP, some advantages and disadvantages along with potential future applications, and provide photographs and technical drawings of all the customised equipment used for the preparation. Using mice or rats, the WHBP is an in situ experimental model that is perfused via an extracorporeal circuit benefitting from unprecedented surgical access, mechanical stability of the brain for whole cell recording and an uncompromised use of pharmacological agents akin to in vitro approaches. The preparation has revealed novel mechanistic insights into, for example, the generation of distinct respiratory rhythms, the neurogenesis of sympathetic activity, coupling between respiration and the heart and circulation, hypothalamic and spinal control mechanisms, and peripheral and central chemoreceptor mechanisms. Insights have been gleaned into diseases such as hypertension, heart failure and sleep apnoea. Findings from the in situ preparation have been ratified in conscious in vivo animals and when tested have translated to humans. We conclude by discussing potential future applications of the WHBP including two‐photon imaging of peripheral and central nervous systems and adoption of pharmacogenetic tools that will improve our understanding of physiological mechanisms and reveal novel mechanisms that may guide new treatment strategies for cardiorespiratory diseases.

Transcriptomes of electrophysiologically recorded Dbx1-derived respiratory neurons of the preBötzinger complex in neonatal mice

Breathing depends on interneurons in the preBötzinger complex (preBötC) derived from Dbx1-expressing precursors. Here we investigate whether rhythm- and pattern-generating functions reside in discrete classes of Dbx1 preBötC neurons. In a slice model of breathing with ~ 5 s cycle period, putatively rhythmogenic Type-1 Dbx1 preBötC neurons activate 100-300 ms prior to Type-2 neurons, putatively specialized for output pattern, and 300-500 ms prior to the inspiratory motor output. We sequenced Type-1 and Type-2 transcriptomes and identified differential expression of 123 genes including ionotropic receptors (Gria3, Gabra1) that may explain their preinspiratory activation profiles and Ca2+ signaling (Cracr2a, Sgk1) involved in inspiratory and sigh bursts. Surprisingly, neuropeptide receptors that influence breathing (e.g., µ-opioid and bombesin-like peptide receptors) were only sparsely expressed, which suggests that cognate peptides and opioid drugs exert their profound effects on a small fraction of the preBötC core. These data in the public domain help explain the neural origins of breathing.

Respiratory rhythm and pattern generation: Brainstem cellular and circuit mechanisms

Breathing movements in mammals are driven by rhythmic neural activity automatically generated within spatially and functionally organized brainstem neural circuits comprising the respiratory central pattern generator (CPG). This chapter reviews up-to-date experimental information and theoretical studies of the cellular and circuit mechanisms of respiratory rhythm and pattern generation operating within critical components of this CPG in the lower brainstem. Over the past several decades, there have been substantial advances in delineating the spatial architecture of essential medullary regions and their regional cellular and circuit properties required to understand rhythm and pattern generation mechanisms. A fundamental concept is that the circuits in these regions have rhythm-generating capabilities at multiple cellular and circuit organization levels. The regional cellular properties, circuit organization, and control mechanisms allow flexible expression of neural activity patterns for a repertoire of respiratory behaviors under various physiologic conditions that are dictated by requirements for homeostatic regulation and behavioral integration. Many mechanistic insights have been provided by computational modeling studies driven by experimental results and have advanced understanding in the field. These conceptual and theoretical developments are discussed.

Pharmacological Modulation and (Patho)Physiological Roles of TRPM4 Channel-Part 2: TRPM4 in Health and Disease

Transient receptor potential melastatin 4 (TRPM4) is a unique member of the TRPM protein family and, similarly to TRPM5, is Ca²⁺ sensitive and permeable for monovalent but not divalent cations. It is widely expressed in many organs and is involved in several functions; it regulates membrane potential and Ca²⁺ homeostasis in both excitable and non-excitable cells. This part of the review discusses the currently available knowledge about the physiological and pathophysiological roles of TRPM4 in various tissues. These include the physiological functions of TRPM4 in the cells of the Langerhans islets of the pancreas, in various immune functions, in the regulation of vascular tone, in respiratory and other neuronal activities, in chemosensation, and in renal and cardiac physiology. TRPM4 contributes to pathological conditions such as overactive bladder, endothelial dysfunction, various types of malignant diseases and central nervous system conditions including stroke and injuries as well as in cardiac conditions such as arrhythmias, hypertrophy, and ischemia-reperfusion injuries. TRPM4 claims more and more attention and is likely to be the topic of research in the future.

Trpm5 channels encode bistability of spinal motoneurons and ensure motor control of hindlimbs in mice

Bistable motoneurons of the spinal cord exhibit warmth-activated plateau potential driven by Na⁺ and triggered by a brief excitation. The thermoregulating molecular mechanisms of bistability and their role in motor functions remain unknown. Here, we identify thermosensitive Na⁺-permeable Trpm5 channels as the main molecular players for bistability in mouse motoneurons. Pharmacological, genetic or computational inhibition of Trpm5 occlude bistable-related properties (slow afterdepolarization, windup, plateau potentials) and reduce spinal locomotor outputs while central pattern generators for locomotion operate normally. At cellular level, Trpm5 is activated by a ryanodine-mediated Ca²⁺ release and turned off by Ca²⁺ reuptake through the sarco/endoplasmic reticulum Ca²⁺-ATPase (SERCA) pump. Mice in which Trpm5 is genetically silenced in most lumbar motoneurons develop hindlimb paresis and show difficulties in executing high-demanding locomotor tasks. Overall, by encoding bistability in motoneurons, Trpm5 appears indispensable for producing a postural tone in hindlimbs and amplifying the locomotor output.

Dynamics of ramping bursts in a respiratory neuron model

Intensive computational and theoretical work has led to the development of multiple mathematical models for bursting in respiratory neurons in the pre-Bötzinger Complex (pre-BötC) of the mammalian brainstem. Nonetheless, these previous models have not captured the pre-inspiratory ramping aspects of these neurons' activity patterns, in which relatively slow tonic spiking gradually progresses to faster spiking and a full-blown burst, with a corresponding gradual development of an underlying plateau potential. In this work, we show that the incorporation of the dynamics of the extracellular potassium ion concentration into an existing model for pre-BötC neuron bursting, along with some parameter adjustments, suffices to induce this ramping behavior. Using fast-slow decomposition, we show that this activity can be considered as a form of parabolic bursting, but with burst termination at a homoclinic bifurcation rather than as a SNIC bifurcation. We also investigate the parameter-dependence of these solutions and show that the proposed model yields a greater dynamic range of burst frequencies, durations, and duty cycles than those produced by other models in the literature.

TRPM4 mediates a subthreshold membrane potential oscillation in respiratory chemoreceptor neurons that drives pacemaker firing and breathing

Brainstem networks that control regular tidal breathing depend on excitatory drive, including from tonically active, CO₂/H⁺-sensitive neurons of the retrotrapezoid nucleus (RTN). Here, we examine intrinsic ionic mechanisms underlying the metronomic firing activity characteristic of RTN neurons. In mouse brainstem slices, large-amplitude membrane potential oscillations are evident in synaptically isolated RTN neurons after blocking action potentials. The voltage-dependent oscillations are abolished by sodium replacement; blocking calcium channels (primarily L-type); chelating intracellular Ca²⁺; and inhibiting TRPM4, a Ca²⁺-dependent cationic channel. Likewise, oscillation voltage waveform currents are sensitive to calcium and TRPM4 channel blockers. Extracellular acidification and serotonin (5-HT) evoke membrane depolarization that augments TRPM4-dependent oscillatory activity and action potential discharge. Finally, inhibition of TRPM4 channels in the RTN of anesthetized mice reduces central respiratory output. These data implicate TRPM4 in a subthreshold oscillation that supports the pacemaker-like firing of RTN neurons required for basal, CO₂-stimulated, and state-dependent breathing.

Transient Receptor Potential Canonical (TRPC) Channels: Then and Now

Twenty-five years ago, the first mammalian Transient Receptor Potential Canonical (TRPC) channel was cloned, opening the vast horizon of the TRPC field. Today, we know that there are seven TRPC channels (TRPC1-7). TRPCs exhibit the highest protein sequence similarity to the Drosophila melanogaster TRP channels. Similar to Drosophila TRPs, TRPCs are localized to the plasma membrane and are activated in a G-protein-coupled receptor-phospholipase C-dependent manner. TRPCs may also be stimulated in a store-operated manner, via receptor tyrosine kinases, or by lysophospholipids, hypoosmotic solutions, and mechanical stimuli. Activated TRPCs allow the influx of Ca2+ and monovalent alkali cations into the cytosol of cells, leading to cell depolarization and rising intracellular Ca2+ concentration. TRPCs are involved in the continually growing number of cell functions. Furthermore, mutations in the TRPC6 gene are associated with hereditary diseases, such as focal segmental glomerulosclerosis. The most important recent breakthrough in TRPC research was the solving of cryo-EM structures of TRPC3, TRPC4, TRPC5, and TRPC6. These structural data shed light on the molecular mechanisms underlying TRPCs' functional properties and propelled the development of new modulators of the channels. This review provides a historical overview of the major advances in the TRPC field focusing on the role of gene knockouts and pharmacological tools.

Critical Components for Spontaneous Activity and Rhythm Generation in Spinal Cord Circuits in Culture

Neuronal excitability contributes to rhythm generation in central pattern generating networks (CPGs). In spinal cord CPGs, such intrinsic excitability partly relies on persistent sodium currents (INaP), whereas respiratory CPGs additionally depend on calcium-activated cation currents (ICAN). Here, we investigated the contributions of INaP and ICAN to spontaneous rhythm generation in neuronal networks of the spinal cord and whether they mainly involve Hb9 neurons. We used cultures of ventral and transverse slices from the E13-14 embryonic rodent lumbar spinal cord on multielectrode arrays (MEAs). All cultures showed spontaneous bursts of network activity. Blocking synaptic excitation with the AMPA receptor antagonist CNQX suppressed spontaneous network bursts and left asynchronous intrinsic activity at about 30% of the electrodes. Such intrinsic activity was completely blocked at all electrodes by both the INaP blocker riluzole as well as by the ICAN blocker flufenamic acid (FFA) and the more specific TRPM4 channel antagonist 9-phenanthrol. All three antagonists also suppressed spontaneous bursting completely and strongly reduced stimulus-evoked bursts. Also, FFA reduced repetitive spiking that was induced in single neurons by injection of depolarizing current pulses to few spikes. Other antagonists of unspecific cation currents or calcium currents had no suppressing effects on either intrinsic activity (gadolinium chloride) or spontaneous bursting (the TRPC channel antagonists clemizole and ML204 and the T channel antagonist TTA-P2). Combined patch-clamp and MEA recordings showed that Hb9 interneurons were activated by network bursts but could not initiate them. Together these findings suggest that both INaP through Na+-channels and ICAN through putative TRPM4 channels contribute to spontaneous intrinsic and repetitive spiking in spinal cord neurons and thereby to the generation of network bursts.

Complex bursting dynamics in an embryonic respiratory neuron model

Pre-Bötzinger complex (pre-BötC) network activity within the mammalian brainstem controls the inspiratory phase of the respiratory rhythm. While bursting in pre-BötC neurons during the postnatal period has been extensively studied, less is known regarding inspiratory pacemaker neuron behavior at embryonic stages. Recent data in mouse embryo brainstem slices have revealed the existence of a variety of bursting activity patterns depending on distinct combinations of burst-generating I_NaP and I_CAN conductances. In this work, we consider a model of an isolated embryonic pre-BötC neuron featuring two distinct bursting mechanisms. We use methods of dynamical systems theory, such as phase plane analysis, fast-slow decomposition, and bifurcation analysis, to uncover mechanisms underlying several different types of intrinsic bursting dynamics observed experimentally including several forms of plateau bursts, bursts involving depolarization block, and various combinations of these patterns. Our analysis also yields predictions about how changes in the balance of the two bursting mechanisms contribute to alterations in an inspiratory pacemaker neuron activity during prenatal development.

Evaluating the Burstlet Theory of Inspiratory Rhythm and Pattern Generation

The preBötzinger complex (preBötC) generates the rhythm and rudimentary motor pattern for inspiratory breathing movements. Here, we test “burstlet” theory (Kam et al., 2013a), which posits that low amplitude burstlets, subthreshold from the standpoint of inspiratory bursts, reflect the fundamental oscillator of the preBötC. In turn, a discrete suprathreshold process transforms burstlets into full amplitude inspiratory bursts that drive motor output, measurable via hypoglossal nerve (XII) discharge in vitro. We recap observations by Kam and Feldman in neonatal mouse slice preparations: field recordings from preBötC demonstrate bursts and concurrent XII motor output intermingled with lower amplitude burstlets that do not produce XII motor output. Manipulations of excitability affect the relative prevalence of bursts and burstlets and modulate their frequency. Whole-cell and photonic recordings of preBötC neurons suggest that burstlets involve inconstant subsets of rhythmogenic interneurons. We conclude that discrete rhythm- and pattern-generating mechanisms coexist in the preBötC and that burstlets reflect its fundamental rhythmogenic nature.

Effects of persistent sodium current blockade in respiratory circuits depend on the pharmacological mechanism of action and network dynamics

The mechanism(s) of action of most commonly used pharmacological blockers of voltage-gated ion channels are well understood; however, this knowledge is rarely considered when interpreting experimental data. Effects of blockade are often assumed to be equivalent, regardless of the mechanism of the blocker involved. Using computer simulations, we demonstrate that this assumption may not always be correct. We simulate the blockade of a persistent sodium current (I_NaP), proposed to underlie rhythm generation in pre-Bötzinger complex (pre-BötC) respiratory neurons, via two distinct pharmacological mechanisms: (1) pore obstruction mediated by tetrodotoxin and (2) altered inactivation dynamics mediated by riluzole. The reported effects of experimental application of tetrodotoxin and riluzole in respiratory circuits are diverse and seemingly contradictory and have led to considerable debate within the field as to the specific role of I_NaP in respiratory circuits. The results of our simulations match a wide array of experimental data spanning from the level of isolated pre-BötC neurons to the level of the intact respiratory network and also generate a series of experimentally testable predictions. Specifically, in this study we: (1) provide a mechanistic explanation for seemingly contradictory experimental results from in vitro studies of I_NaP block, (2) show that the effects of I_NaP block in in vitro preparations are not necessarily equivalent to those in more intact preparations, (3) demonstrate and explain why riluzole application may fail to effectively block I_NaP in the intact respiratory network, and (4) derive the prediction that effective block of I_NaP by low concentration tetrodotoxin will stop respiratory rhythm generation in the intact respiratory network. These simulations support a critical role for I_NaP in respiratory rhythmogenesis in vivo and illustrate the importance of considering mechanism when interpreting and simulating data relating to pharmacological blockade.

The application of pharmacological agents that affect transmembrane ionic currents in neurons is a commonly used experimental technique. A simplistic interpretation of experiments involving these agents suggests that antagonist application removes the impacted current and that subsequently observed changes in activity are attributable to the loss of that current’s effects. The more complex reality, however, is that different drugs may have distinct mechanisms of action, some corresponding not to a removal of a current but rather to a changing of its properties. We use computational modeling to explore the implications of the distinct mechanisms associated with two drugs, riluzole and tetrodotoxin, that are often characterized as sodium channel blockers. Through this approach, we offer potential explanations for disparate findings observed in experiments on neural respiratory circuits and show that the experimental results are consistent with a key role for the persistent sodium current in respiratory rhythm generation.

Robustness of respiratory rhythm generation across dynamic regimes

A central issue in the study of the neural generation of respiratory rhythms is the role of the intrinsic pacemaking capabilities that some respiratory neurons exhibit. The debate on this issue has occurred in parallel to investigations of interactions among respiratory network neurons and how these contribute to respiratory behavior. In this computational study, we demonstrate how these two issues are inextricably linked. We use simulations and dynamical systems analysis to show that once a conditional respiratory pacemaker, which can be tuned across oscillatory and non-oscillatory dynamic regimes in isolation, is embedded into a respiratory network, its dynamics become masked: the network exhibits similar dynamical properties regardless of the conditional pacemaker node's tuning, and that node's outputs are dominated by network influences. Furthermore, the outputs of the respiratory central pattern generator as a whole are invariant to these changes of dynamical properties, which ensures flexible and robust performance over a wide dynamic range.

Biophysical mechanisms in the mammalian respiratory oscillator re-examined with a new data-driven computational model

An autorhythmic population of excitatory neurons in the brainstem pre-Bötzinger complex is a critical component of the mammalian respiratory oscillator. Two intrinsic neuronal biophysical mechanisms—a persistent sodium current (INaP) and a calcium-activated non-selective cationic current (ICAN)—were proposed to individually or in combination generate cellular- and circuit-level oscillations, but their roles are debated without resolution. We re-examined these roles in a model of a synaptically connected population of excitatory neurons with ICAN and INaP. This model robustly reproduces experimental data showing that rhythm generation can be independent of ICAN activation, which determines population activity amplitude. This occurs when ICAN is primarily activated by neuronal calcium fluxes driven by synaptic mechanisms. Rhythm depends critically on INaP in a subpopulation forming the rhythmogenic kernel. The model explains how the rhythm and amplitude of respiratory oscillations involve distinct biophysical mechanisms.

Modeling breathing rhythms

Computational models are helping researchers to understand how certain properties of neurons contribute to respiratory rhythms.

Trpm4 ion channels in pre-Bötzinger complex interneurons are essential for breathing motor pattern but not rhythm

Inspiratory breathing movements depend on pre-Bötzinger complex (preBötC) interneurons that express calcium (Ca2+)-activated nonselective cationic current (ICAN) to generate robust neural bursts. Hypothesized to be rhythmogenic, reducing ICAN is predicted to slow down or stop breathing; its contributions to motor pattern would be reflected in the magnitude of movements (output). We tested the role(s) of ICAN using reverse genetic techniques to diminish its putative ion channels Trpm4 or Trpc3 in preBötC neurons in vivo. Adult mice transduced with Trpm4-targeted short hairpin RNA (shRNA) progressively decreased the tidal volume of breaths yet surprisingly increased breathing frequency, often followed by gasping and fatal respiratory failure. Mice transduced with Trpc3-targeted shRNA survived with no changes in breathing. Patch-clamp and field recordings from the preBötC in mouse slices also showed an increase in the frequency and a decrease in the magnitude of preBötC neural bursts in the presence of Trpm4 antagonist 9-phenanthrol, whereas the Trpc3 antagonist pyrazole-3 (pyr-3) showed inconsistent effects on magnitude and no effect on frequency. These data suggest that Trpm4 mediates ICAN, whose influence on frequency contradicts a direct role in rhythm generation. We conclude that Trpm4-mediated ICAN is indispensable for motor output but not the rhythmogenic core mechanism of the breathing central pattern generator.

Kinetic properties of persistent Na+ current orchestrate oscillatory bursting in respiratory neurons

The rhythmic pattern of breathing depends on the pre-Bötzinger complex (preBötC) in the brainstem, a vital circuit that contains a population of neurons with intrinsic oscillatory bursting behavior. Here, we investigate the specific kinetic properties that enable voltage-gated sodium channels to establish oscillatory bursting in preBötC inspiratory neurons, which exhibit an unusually large persistent Na⁺ current (I_NaP). We first characterize the kinetics of I_NaP in neonatal rat brainstem slices in vitro, using whole-cell patch-clamp and computational modeling, and then test the contribution of I_NaP to rhythmic bursting in live neurons, using the dynamic clamp technique. We provide evidence that subthreshold activation, persistence at suprathreshold potentials, slow inactivation, and slow recovery from inactivation are kinetic features of I_NaP that regulate all aspects of intrinsic rhythmic bursting in preBötC neurons. The slow and cumulative inactivation of I_NaP during the burst active phase controls burst duration and termination, while the slow recovery from inactivation controls the duration of the interburst interval. To demonstrate this mechanism, we develop a Markov state model of I_NaP that explains a comprehensive set of voltage clamp data. By adding or subtracting a computer-generated I_NaP from a live neuron via dynamic clamp, we are able to convert nonbursters into intrinsic bursters, and vice versa. As a control, we test a model with inactivation features removed. Adding noninactivating I_NaP into nonbursters results in a pattern of random transitions between sustained firing and quiescence. The relative amplitude of I_NaP is the key factor that separates intrinsic bursters from nonbursters and can change the fraction of intrinsic bursters in the preBötC. I_NaP could thus be an important target for regulating network rhythmogenic properties.