Biophysical mechanisms in the mammalian respiratory oscillator re-examined with a new data-driven computational model
Abstract
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 () and a calcium-activated non-selective cationic current ()—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 and . This model robustly reproduces experimental data showing that rhythm generation can be independent of activation, which determines population activity amplitude. This occurs when is primarily activated by neuronal calcium fluxes driven by synaptic mechanisms. Rhythm depends critically on in a subpopulation forming the rhythmogenic kernel. The model explains how the rhythm and amplitude of respiratory oscillations involve distinct biophysical mechanisms.
https://doi.org/10.7554/eLife.41555.001Introduction
Defining cellular and circuit mechanisms generating the vital rhythm of breathing in mammals remains a fundamental unsolved problem of wide-spread interest in neurophysiology (Richter and Smith, 2014; Del Negro et al., 2018; Ramirez and Baertsch, 2018a), with potentially far-reaching implications for understanding mechanisms of oscillatory circuit activity and rhythmic motor pattern generation in neural systems (Marder and Calabrese, 1996; Buzsaki, 2006; Grillner, 2006; Kiehn, 2006). The brainstem pre-Bötzinger complex (pre-BötC) region (Smith et al., 1991) located in the ventrolateral medulla oblongata is established to contain circuits essential for respiratory rhythm generation (Smith et al., 2013; Del Negro et al., 2018), but the operational cellular biophysical and circuit synaptic mechanisms are continuously debated. Pre-BötC excitatory neurons and circuits have autorhythmic properties and drive motor circuits that can be isolated and remain rhythmically active in living rodent brainstem slices in vitro. Numerous experimental and theoretical analyses have focused on the rhythmogenic mechanisms operating in these in vitro conditions to provide insight into biophysical and circuit processes involved, with potential relevance for rhythm generation during breathing in vivo (Feldman and Del Negro, 2006; Lindsey et al., 2012; Richter and Smith, 2014; Ramirez and Baertsch, 2018b). The ongoing rhythmic activity in vitro has been suggested to arise from a variety of cellular and circuit biophysical mechanisms including from a subset(s) of intrinsically bursting neurons which, through excitatory synaptic interactions, recruit and synchronize neurons within the network (pacemaker-network models) (Butera et al., 1999b; Ramirez et al., 2004; Toporikova and Butera, 2011; Chevalier et al., 2016), or as an emergent network property through recurrent excitation (e.g. Rekling and Feldman, 1998; Jasinski et al., 2013) and/or synaptic depression (group pacemaker model) (Rubin et al., 2009a; Del Negro et al., 2010).
From these previous analyses, involvement of two possible cellular-level biophysical mechanisms have been proposed. One based on a slowly inactivating persistent sodium current () (Butera et al., 1999a), and the other on a calcium-activated non-selective cation current () coupled to intracellular calcium () dynamics (for reviews see Rybak et al., 2014; Del Negro et al., 2010), or a combination of both mechanisms (Thoby-Brisson and Ramirez, 2001; Jasinski et al., 2013; Peña et al., 2004). Despite the extensive experimental and theoretical investigations of these sodium- and calcium-based mechanisms, the actual roles of , and the critical source(s) of transients in the pre-BötC are still unresolved. Furthermore, in pre-BötC circuits, the process of rhythm generation must be associated with an amplitude of circuit activity sufficient to drive downstream circuits to produce adequate inspiratory motor output. Biophysical mechanisms involved in generating the amplitude of pre-BötC circuit activity have also not been established.
is proposed to mediate an essential oscillatory burst-generating mechanism since pharmacologically inhibiting abolishes intrinsic neuronal rhythmic bursting as well as pre-BötC circuit inspiratory activity and rhythmic inspiratory motor output in vitro (Koizumi et al., 2008; Toporikova et al., 2015), although some studies suggest that block of both and are necessary to disrupt rhythmogenesis in vitro (Peña et al., 2004). Theoretical models of cellular and circuit activity based on -dependent bursting mechanisms closely reproduce experimental observations such as voltage-dependent frequency control, spike-frequency adaptation during bursts, and pattern formation of inspiratory motor output (Butera et al., 1999b; Pierrefiche et al., 2004; Smith et al., 2007). This indicates the plausibility of -dependent rhythm generation.
In the pre-BötC, was originally postulated to underlie intrinsic pacemaker-like oscillatory bursting at the cellular level and contribute to circuit-level rhythm generation, since intrinsic bursting in a subset of neurons in vitro was found to be terminated by the inhibitor flufenamic acid (FFA) (Peña et al., 2004). Furthermore, inhibition of in the pre-BötC reduces the amplitude of the rhythmic depolarization (inspiratory drive potential) driving neuronal bursting and can eliminate inspiratory motor activity in vitro (Pace et al., 2007). became the centerpiece of the 'group pacemaker' model for rhythm generation, in which this conductance was proposed to be activated by inositol trisphosphate (IP3) receptor/ER-mediated intracellular calcium fluxes initiated via glutamatergic metabotropic receptor-mediated signaling in the pre-BötC excitatory circuits (Del Negro et al., 2010). The molecular correlate of was postulated to be the transient receptor potential channel M4 (TRPM4) (Mironov, 2008; Pace et al., 2007), one of the two known Ca2+-activated TRP channels (Guinamard et al., 2010; Ullrich et al., 2005), or alternatively, by the transient receptor potential channels C3/7 (TRPC3/7) (Ben-Mabrouk and Tryba, 2010); however, these latter channels are not known to be Ca2+-activated (Clapham, 2003). TRPM4 and TRPC3 have now been identified by immunolabeling and RNA expression profiling in pre-BötC inspiratory neurons in vitro (Koizumi et al., 2018).
Investigations into the sources of intracellular Ca2+ activating /TRPM4 suggested that (1) somatic calcium transients from voltage-gated sources do not contribute to the inspiratory drive potential (Morgado-Valle et al., 2008), (2) IP3/ER-mediated intracellular Ca2+ release does not contribute to inspiratory rhythm generation in vitro, and (3) in the dendrites calcium transients may be triggered by excitatory synaptic inputs and travel in a wave propagated to the soma (Mironov, 2008). Theoretical studies have demonstrated the plausibility of --dependent bursting (Rubin et al., 2009b; Toporikova and Butera, 2011); however, these models omit and/or depend on additional unproven mechanisms to generate intracellular calcium oscillations to provide burst termination, such as IP3-dependent calcium-induced calcium release (Toporikova and Butera, 2011), partial depolarization block of action potentials (Rubin et al., 2009a), and the Na+/K+ pump (Jasinski et al., 2013). Interestingly, pharmacological inhibition of /TRPM4 has been shown to produce large reductions in the amplitude of pre-BötC inspiratory neuron population activity with essentially no, or minor perturbations of inspiratory rhythm (Peña et al., 2004). These observations constrain the role of , and require theoretical re-examination of pre-BötC neuronal conductance mechanisms and network dynamics, particularly how rhythm generation mechanisms can be independent of -dependent mechanisms that regulate the amplitude of network activity.
In this theoretical study, we examine the role of in pre-BötC excitatory circuits by considering two plausible mechanisms of intracellular calcium fluxes: (1) from voltage-gated and (2) from synaptically activated sources. We deduce that is primarily activated by calcium transients that are coupled to rhythmic excitatory synaptic inputs originating from -dependent bursting inspiratory neurons. Additionally, we show that contributes to the inspiratory drive potential by mirroring the excitatory synaptic current. This concept is consistent with a mechanism underlying generation of the inspiratory drive potential involving a synaptic-based activation described in previous work (Pace et al., 2007; Rubin et al., 2009a). Our model explains the experimental observations obtained from in vitro neonatal rodent slices isolating the pre-BötC, showing large reductions in network activity amplitude by inhibiting /TRPM4 without perturbations of inspiratory rhythm generation in pre-BötC excitatory circuits in vitro. The model supports the concept that activation in a subpopulation of pre-BötC excitatory neurons is critically involved in amplifying synaptic drive from a subset of neurons whose rhythmic bursting is critically dependent on and forms the kernel for rhythm generation in vitro. The model suggests how the functions of generating the rhythm and amplitude of inspiratory oscillations in pre-BötC excitatory circuits are determined by distinct biophysical mechanisms.
Results
variation has opposite effects on amplitude and frequency of network bursting in the and models
Experimental work (Peña et al., 2004) has demonstrated that pharmacological inhibition of /TRPM4 in the pre-BötC from in vitro neonatal mouse/rat slice preparations, strongly reduces the amplitude of (or completely eliminates) the inspiratory hypoglossal (XII) motor output, as well as the amplitude of pre-BötC excitatory circuit activity that is highly correlated with the decline of XII activity, while having relatively little effect on inspiratory burst frequency. Here, we systematically examine in our model the relationship between conductance () and the amplitude and frequency of network activity for voltage-gated () and synaptically activated sources () of intracellular calcium in a heterogeneous network of 100 synaptically coupled single-compartment pre-BötC model excitatory neurons. In addition to voltage-gated and synaptically activated calcium currents, each model neuron incorporates voltage-gated action potential generating currents, as well as , , leak, and excitatory synaptic currents adapted from the conductance-based biophysical model of Jasinski et al. (2013) (see Materials and methods for a full model description). We note that in our full model, is activated by both voltage-gated and synaptic mechanisms, consistent with experimental results (Thoby-Brisson and Ramirez, 2001; Pace et al., 2007; Morgado-Valle et al., 2008). Initial separate consideration of the and provides a means to deduce the relative contribution of these two general sources of intracellular calcium to activation. We found that reduction of drives opposing effects on network activity amplitude (spike/s/neuron) and frequency that are dependent on the source of intracellular calcium transients (Figure 1). The network activity amplitude is a measure of the average neuronal population firing rate and is defined by the number of spikes generated by the network per 50 ms bin divided by the number of neurons in the network. In the network, where calcium influx is generated exclusively from voltage-gated calcium channels, increasing has no effect on amplitude but increases the frequency of network oscillations (Figure 1A, C, D). Conversely, in the network where calcium influx is generated exclusively by excitatory synaptic input, increasing strongly increases the amplitude and slightly decreases the oscillation frequency (Figure 1B, C, D).
Effects of subthreshold activation of on network frequency
In -dependent bursting neurons in the pre-BötC, bursting frequency depends on their excitability (i.e. baseline membrane potential) which can be controlled in different ways, for example by directly injecting a depolarizing current (Smith et al., 1991; Del Negro et al., 2005; Yamanishi et al., 2018) or varying the conductance and/or reversal potentials of some ionic channels (Butera et al., 1999a). Due to their relatively short duty cycle, the bursting frequency in these neurons is largely determined by the interburst interval, defined as the time between the end of one burst and the start of the next. During the burst, slowly inactivates (Butera et al., 1999a) resulting in burst termination and abrupt neuronal hyperpolarization. The interburst interval is then determined by the amount of time required for to recover from inactivation and return the membrane potential back to the threshold for burst initiation. This process is governed by the kinetics of inactivation gating variable . Higher neuronal excitability reduces the value of required to initiate bursting. Consequently, the time required to reach this value is decreased, which results in a shorter interburst interval and increased frequency.
To understand how changing affects network bursting frequency, we quantified the values of averaged over all rhythm-generating pacemaker neurons immediately preceding each network burst and, also, the average values between the bursts in the and networks (Figure 2). In the network, as modeled remains residually activated between the bursts thus creating the background calcium concentration which partially activates . Therefore, between the bursts functions as a depolarizing leak current. Consistently, we found that in the network increasing increases (Figure 2A), progressively depolarizing the network, which reduces the threshold for burst initiation (Figure 2B) and, thus, increases network oscillation frequency (Figure 1D).
In the model, the intracellular calcium depletes entirely during the interburst interval. Consequently, increasing has no effect on (Figure 2A) and frequency is essentially unaffected (Figure 1D).
Changes in network activity amplitude are driven by recruitment of neurons
As previously stated, the network activity amplitude is defined as the total number of spikes produced by the network per a time bin divided by the number of neurons in the network. Consequently, changes in amplitude can only occur by increasing the number of neurons participating in bursts (recruitment) and/or increasing the firing rate of the recruited neurons. To analyze changes in amplitude, we quantified the number of recruited neurons (Figure 3A) and the average spike frequency in recruited neurons (Figure 3B) as a function of for both network models. In the network, increasing increases the number of recruited neurons (Figure 3A), but decreases the average spiking frequency in recruited neurons (Figure 3B) which, together result in no change in amplitude (Figure 1C). In the network, increasing strongly increases the number of recruited neurons (Figure 3A) and increases the spike frequency of recruited neurons (Figure 3B) resulting in a large increase in network activity amplitude (Figure 3C).
Manipulating in the model is qualitatively equivalent to changing the strength of synaptic interactions
Since changes in in the model primarily affect network activity amplitude through recruitment of neurons, and the network activity amplitude strongly depends on the strength of synaptic interactions, we next examined the relationship between , synaptic strength and network activity amplitude and frequency (Figure 4). Synaptic strength is defined as the number of neurons multiplied by the synaptic connection probability multiplied by the average weight of synaptic connections (), where the weight of synaptic connections ranges from to We found that the effects of varying or the synaptic strength on network activity amplitude and frequency are qualitatively equivalent in the network which is indicated by symmetry of the heat plots (across the X=Y line) in Figure 4A,B. This symmetry results from the fact that the effective strength of synaptic interactions in the network is roughly proportional to a product of the synaptic strength and . A transition from bursting to tonic spiking occurs when this effective excitation exceeds a certain critical value. This is why the bifurcation curve corresponding to a transition from rhythmic bursting to tonic spiking (a boundary between yellow and black in Figure 4A) looks like a hyperbola ().
We further investigated and compared the effect of reducing or the synaptic strength on network activity amplitude and frequency as well as the effects on the recruitment of neurons not involved in rhythm generation (Figure 4C–F). To make this comparison, we picked a starting point in the 2D parameter space of and synaptic strength where the network is bursting. Then in separate simulations, we gradually reduced either or the synaptic strength to zero. We show that reducing either or the synaptic strength have very similar effects on network activity amplitude and frequency (Figure 4C,D). Furthermore, de-recruitment of neurons in both cases is nearly identical (Figure 4E,F). Reducing either or the synaptic strength decreases the excitatory input to the neurons during network oscillations which is a major component of the inspiratory drive potential. Therefore, in the network, manipulations of will affect the strength of the inspiratory drive potential in the rhythmic inspiratory neurons in a way that is equivalent to changing the synaptic strength of the network. In contrast, manipulations of in the network will only slightly affect the inspiratory drive potential due to changes in the average firing rate of active neurons (see Figure 3B).
Robustness of amplitude and frequency effects
We also examined if the effects are conserved in both the and networks over a range of network parameters. To test this, we investigated the dependence of network activity amplitude and frequency on and average synaptic strength for and networks with high () and low () connection probabilities, and high (), medium () and low () strengths of calcium sources (Figure 5 and 6). We found that changing the synaptic connection probability and changing the strength of the calcium sources has no effect on the general relationship between and the amplitude or frequency of bursts in the or networks. In other words, the general effect of increasing on amplitude and frequency is conserved in both networks regardless of the synaptic connection probability or strength of the calcium sources. Increasing the strength of the calcium sources does, however, affects the range of possible values where both networks produce rhythmic activity.
To summarize, in the model, increasing increases frequency, through increased excitability but has no effect on amplitude. In contrast, in the model, increasing slightly decreases frequency and increases amplitude. In this case, increasing acts as a mechanism to increase the inspiratory drive potential and recruit previously silent neurons. Additionally, these features of the and models are robust and conserved across a wide range of network parameters.
Intracellular calcium transients activating primarily result from synaptically activated sources
In experiments where /TRPM4 was blocked by bath application of FFA or 9-phenanthrol (Koizumi et al., 2018) in vitro, the amplitude of network oscillations was strongly reduced and their frequency remained unchanged or was reduced insignificantly. Our model revealed that the effects of blockade on amplitude and frequency depend on the source(s) of intracellular calcium (see Figures 1 and 2). If the calcium influx is exclusively voltage-gated, our model predicts that blockade will have no effect on amplitude but reduce the frequency. In contrast, if the calcium source is exclusively synaptically gated, our model predicts that blocking will strongly reduce the amplitude and slightly increase the frequency. Therefore, a multi-fold decrease in amplitude, seen experimentally, is consistent with the synaptically driven calcium influx mechanisms, while nearly constant bursting frequency may be due to calcium influx through both voltage- and synaptically gated channels. Following the predictions above, to reproduce experimental data, we incorporated both mechanisms in the full model and inferred their individual contributions by finding the best fit. We found that the best match is observed (Figure 7) if synaptically mediated and voltage-gated calcium influxes comprise about 95% and 5% of the total calcium influx, respectively. We note that some experiments have shown a small (~20%) reduction of inspiratory burst frequency accompanying larger reductions in pre-BötC population activity amplitude with block (Peña et al., 2004). In our model, such perturbations of frequency can occur if the contribution of voltage-gated calcium influxes is larger than indicated above, or if neuronal background calcium concentration which partially activates is higher than specified in the model.
-dependent and --sensitive intrinsic bursting
In our model, we included , as well as voltage-gated and synaptic mechanisms of influx. Activation of by is the equivalent mechanism used in computational group-pacemaker models (Rubin et al., 2009a; Song et al., 2015). Rhythmic burst generation and termination in our model, however, are dependent on (Butera et al., 1999a). We investigated the sensitivity of intrinsic bursting in our model to and calcium channel blockade (Figure 8). Intrinsic bursting was identified in neurons by zeroing the synaptic weights to simulate synaptic blockade. and blockade was simulated by setting and to . We found that after decoupling the network () a subset of neurons remained rhythmically active (7% in this simulation) and that these were all neurons with a high conductance. In these rhythmically active neurons, bursting was abolished in all neurons by blockade. Interestingly, blockade applied before block abolished intrinsic bursting in three of the seven neurons and block applied afterwards abolished intrinsic bursting in the remaining four neurons. Although only one rhythmogenic (-based) mechanism exists in this model, bursting in a subset of these intrinsically bursting neurons is calcium sensitive, consistent with experimental observations of calcium-sensitive intrinsic bursters (Thoby-Brisson and Ramirez, 2001; Del Negro et al., 2005; Peña et al., 2004). In calcium-sensitive bursters, Ca2+ blockade in our model abolishes bursting by reducing the intracellular calcium concentration and, hence, activation, which ultimately reduces excitability. We note that in our model the numbers of -dependent and calcium-sensitive intrinsic bursters will vary depending on the mean and width of the distribution and the background intracellular calcium concentration.
The rhythmogenic kernel
Our simulations have shown that the primary role of is amplitude but not oscillation frequency modulation with little or no effect on network activity frequency. Here we examined the neurons that remain active and maintain rhythm after blockade (Figure 9). We found that the neurons that remain active are primarily neurons with the highest and that bursting in these neurons is dependent on . Some variability exists and neurons with relatively low value can remain active due to synaptic interactions while a neuron with a slightly higher without sufficient synaptic input may become silent. These neurons that remain active after compete blockade of form a -dependent kernel of a rhythm generating circuit.
Intracellular calcium transients and network activity during inhibition of /TRPM4
Dynamic calcium imaging has been used to assess activity of the population of pre-BöC excitatory neurons as well as individual pre-BötC neurons during pharmacological inhibition of /TRPM4 in vitro (Koizumi et al., 2018). These experiments indicate that network output activity amplitude and pre-BötC excitatory population-level intracellular calcium transients are highly correlated while the network oscillation frequency is not significantly perturbed. Interestingly, during /TRPM4 block, changes in calcium transients of individual neurons can differ significantly from the average population-level calcium transients. To assess if our model is consistent with these experimental results and gain additional insight into intracellular calcium dynamics during network activity, we analyzed simultaneous changes in the amplitude of network neuronal spiking activity, the average intracellular calcium concentration ( of all network neurons, as well as of individual neurons, with different network connection probabilities () during simulated block (Figure 10). We found that regardless of , the network activity amplitude and average intracellular calcium concentration are highly correlated (Figure 10A,B). has no effect on the relationship between amplitude, calcium transients at the network level, or network oscillation frequency provided that the synaptic strength remains constant (). does, however, affect the change in the peak in individual neurons. In a network with a high connection probability () the synaptic current/calcium transient is nearly identical for all neurons and therefore the change in during blockade is approximately the same for each neuron (Figure 10C). In a sparsely connected network, as proposed for the connectivity of the pre-BötC network (Carroll et al., 2013; Carroll and Ramirez, 2013) the synaptic current and calcium influx are more variable and reflect the heterogeneity in spiking frequency of the pre-synaptic neurons (Figure 10). Interestingly, in a network with low connection probability (), the peak transient in some neurons increases when is blocked (Figure 10E), consistent with our experimental results (Koizumi et al., 2018).
Discussion
Establishing cellular and circuit mechanisms generating the rhythm and amplitude of respiratory oscillations in the mammalian brainstem pre-BötC has remained an unsolved problem of wide-spread interest in neurophysiology since this structure, essential for breathing to support mammalian life, was discovered nearly three decades ago (Smith et al., 1991). Our objective in this theoretical study was to re-examine and further define contributions of two of the main currently proposed neuronal biophysical mechanisms operating in pre-BötC excitatory circuits, specifically mechanisms involving activated by neuronal calcium fluxes (Thoby-Brisson and Ramirez, 2001; Peña et al., 2004; Del Negro et al., 2005; Mironov, 2008; Rubin et al., 2009a ) and voltage-dependent in the circuit neurons (Butera et al., 1999a; Del Negro et al., 2002; Koizumi and Smith, 2008; Koizumi and Smith, 2008). While these sodium- and calcium-based mechanisms have been studied extensively over the past two-decades and shown experimentally to be integrated in pre-BötC circuits, their actual roles in circuit operation are continuously debated and unresolved (Rybak et al., 2014). Both mechanisms have been proposed to be fundamentally involved in rhythm generation either separately or in combination, as plausibly shown from previous theoretical modeling studies (Butera et al., 1999a; Jasinski et al., 2013; Rubin et al., 2009a; Toporikova and Butera, 2011). Furthermore, the process of rhythm generation in pre-BötC circuits must be associated with an amplitude of excitatory circuit activity sufficient to drive downstream circuits to produce adequate respiratory motor output. Biophysical mechanisms involved in generating excitatory population activity amplitude have also not been established. Our analysis is motivated by the experimental observations obtained from neonatal rodent slices isolating pre-BötC circuits in vitro that inhibition of the endogenously active /TRPM4 strongly reduces the amplitude of network oscillations within pre-BötC circuits but has a small effect on oscillation frequency (Peña et al., 2004; Koizumi et al., 2018). These findings challenge currently proposed -based models for rhythm generation in the isolated pre-BötC and indicate a functional organization of pre-BötC circuits, in terms of oscillatory frequency and amplitude generation, that needs to be defined.
We accordingly analyzed the role of and possible sources of intracellular calcium transients activating this conductance and found that the effect of simulated blockade on amplitude and frequency is highly dependent on the source(s) of intracellular calcium, which is also a central issue to be resolved. In the case where is the primary intracellular calcium source, blockade generates a large reduction in network activity amplitude. In contrast, when is the only intracellular calcium source, blockade has little effect on network activity amplitude and primarily affects the population bursting frequency that is caused by decreased excitability. Additionally, we show that activation of by functions as a mechanism to augment the inspiratory drive potential and amplitude of population activity, and that this effect is similar to increasing the synaptic coupling strength within the network. Therefore, in the case of , blockade of reduces the inspiratory drive potential causing de-recruitment of non-pacemaker rhythmic neurons and reduction of network activity amplitude. In a model where is activated by both and with contributions of 5% and 95% respectively, we show that simulated blockade of generates a large reduction in network population activity amplitude and a slight decrease in frequency. This closely reproduces experimental blockade of /TRPM4 by either 9-phenanthrol or FFA (Figure 7). Finally, we showed that the change in the peak calcium transients for individual neurons during blockade, particularly at relatively low network connection probabilities (), are consistent with experimental data.
Role of in the pre-BötC respiratory network
The hypothesis that is involved in generation of the inspiratory rhythm is based on experimental observations from in vitro mouse medullary slice preparations (Pace et al., 2007; Mironov, 2008; Del Negro et al., 2005; Del Negro et al., 2010; Peña et al., 2004; Thoby-Brisson and Ramirez, 2001), and in silico modeling studies (Jasinski et al., 2013; Rybak et al., 2003; Toporikova and Butera, 2011). Theories of -dependent bursting rely on intracellular signaling mechanisms that have not been well defined.
Two models of -dependent rhythmic bursting in vitro have been proposed and are referred to as the 'dual pacemaker' and 'group pacemaker' models. In the dual pacemaker model, two types of pacemaker neurons are proposed that are either -dependent (riluzole sensitive) or -dependent (sensitive) intrinsic bursters (inspiratory 'pacemaker' neurons) (see Rybak et al., 2014 for review). In this model network, oscillations are thought to originate from these pacemaker neurons which through excitatory synaptic interactions synchronize bursting and drive activity of other rhythmic inspiratory neurons within the pre-BötC. Although pacemaker neurons sensitive to neuronal flux blockade through have been reported (Peña et al., 2004; Thoby-Brisson and Ramirez, 2001, the source and mechanism driving intracellular oscillations has not been definitely delineated. Computational models of -dependent pacemaker neurons rely on mechanisms for burst initiation and termination, for example IP3-dependent oscillations (Toporikova and Butera, 2011; Del Negro et al., 2010), that have been questioned from recent negative experimental results (Beltran-Parrazal et al., 2012; Toporikova et al., 2015). In favor of the dual pacemaker concept, experimental evidence has been presented that pharmacological block of both and are necessary to disrupt rhythmogenesis under normoxic conditions in vitro (Peña et al., 2004). Other experimental results suggest that blocking alone in vitro is sufficient (Koizumi and Smith, 2008; Toporikova et al., 2015), and we have discussed from a theoretical perspective how such discrepancies may depend on network excitability state and connectivity (Jasinski et al., 2013). Our present model defines more clearly parameters such as calcium flux sources and background calcium levels that can influence the relative contributions of in various network states to be further explored in model simulation studies.
In versions of the 'group pacemaker' model (Rubin et al., 2009a; Del Negro et al., 2010; Feldman and Del Negro, 2006), network oscillations are initiated through recurrent synaptic excitation that trigger postsynaptic influx. Subsequent activation generates membrane depolarization (inspiratory drive potential) to drive neuronal bursting. Synaptically triggered influx and the contribution of to the inspiratory drive potential of individual pre-BötC neurons are experimentally supported (Mironov, 2008; Pace et al., 2007); however, the mechanism of burst termination remains unclear. Again, the computational group-pacemaker models that have been explored (Rubin et al., 2009a) rely on speculative mechanisms for burst termination that need to be tested experimentally, and in some cases lack key biophysical features of the pre-BötC neurons such as voltage-dependent frequency control and expression of .
In our model, we showed that blockade of either or synaptic interactions produce qualitatively equivalent effects on network population activity amplitude and frequency when the calcium transients are primarily generated from synaptic sources (Figure 4). Consequently, our model predicts that blockade of or synaptic interactions in the isolated pre-BötC in vitro will produce comparable effects on amplitude and frequency. This is the case as Johnson et al. (1994) showed that gradual blockade of synaptic interactions by low calcium solution significantly decreases network activity amplitude while having little effect on frequency, similar to the experiments where the channel TRPM4 is blocked with 9-phenanthrol (Koizumi et al., 2018). We note that complete blockade of in our model can ultimately abolish synchronized network oscillations due to weakened excitatory synaptic transmission, which results in neuronal de-recruitment and desynchronization of the network, particularly when synaptic strength is low (Figure 4). Thus, plays a critical role in network activity synchronization that determines the ability of the pre-BötC excitatory network to produce rhythmic output, depending on synaptic strength.
Overall, our new model simulations for the isolated pre-BötC excitatory network suggest that the role of /TRPM4 activation is to amplify excitatory synaptic drive in generating the amplitude of inspiratory population activity, essentially independent of the biophysical mechanism generating inspiratory rhythm. We note that the recent experiments have also shown that in the more intact brainstem respiratory network that ordinarily generates patterns of inspiratory and expiratory activity, endogenous activation of /TRPM4 appears to augment the amplitude of both inspiratory and expiratory population activity, and hence these channels are fundamentally involved in inspiratory-expiratory pattern formation (Koizumi et al., 2018).
Intracellular calcium dynamics and network activity during inhibition of /TRPM4
We analyzed the correlation between calcium transients and inspiratory activity of individual inspiratory neurons as well as the entire network, particularly since dynamic calcium imaging has been utilized to assess activity of individual cells and populations of pre-BötC excitatory neurons in vitro during pharmacological inhibition of /TRPM4 (Koizumi et al., 2018). We show that intracellular calcium transients are highly correlated with network and cellular activity across the duration of an blockade simulation, consistent with experimental observations.
Additionally, we examined the relative change in the peak calcium transients in single neurons as a function of conductance. We show that in a subset of neurons the peak calcium transient increases with reduced . This result is surprising but is also supported by the calcium imaging data (Koizumi et al., 2018). This occurs in neurons that receive most of their synaptic input from pacemaker neurons and our analyses suggest this is possible in sparse networks, that is with low connection probability. In pacemaker neurons, blockade leads to a reduction of their excitability resulting in an increased value of inactivation gating variable at the burst onset. Thus, during the burst, the peak action potential frequency and the synaptic output from these neurons is increased with blockade. Consequently, neurons that receive synaptic input from pacemaker neurons will see an increase in their peak calcium transients. In most neurons, however, synaptic input is received primarily from non-pacemaker rhythmic neurons. Since blockade de-recruits these non-pacemaker neurons, the synaptic input and subsequent calcium influx in most of these cells decreases. Therefore, our model predicts that in a sparse network, which has been proposed for pre-BötC network connectivity (Carroll et al., 2013; Carroll and Ramirez, 2013), blocking results in very diverse responses at the cellular level with an overall tendency to reduce intracellular calcium transients such that the amplitude of these transients averaged over the entire population decreases during blockade, while their burst frequency is essentially unchanged, as found from the calcium imaging experiments (Koizumi et al., 2018).
Synaptic calcium sources
Our model suggests that calcium transients in the pre-BötC are coupled to excitatory synaptic input, that is pre-synaptic glutamate release and binding to post-synaptic glutamate receptors triggers calcium entry. The specific mechanisms behind this process are unclear; however, this is likely dependent on specific types of ionotropic or metabotropic glutamate receptors.
There are three subtypes of ionotropic glutamate receptors, N-methyl-D-aspartate (NMDA), Kainate (KAR), and a-amino-3-hydroxy-5-methyl-4-isoxazolepropionic acid (AMPA), all of which are expressed in the pre-BötC (Paarmann et al., 2000) and have varying degrees of calcium permeability. NMDA and AMPA are unlikely candidates for direct involvement in synaptically mediated calcium influx in the pre-BötC. Pharmacological blockade of NMDA receptors does not consistently effect the amplitude or frequency of XII motor output (Lieske and Ramirez, 2006a; Morgado-Valle and Feldman, 2007; Pace et al., 2007) and AMPA receptors in the pre-BötC show high expression of the subunit GluR2, which renders the AMPA ion channel pore impermeable to Ca2+ (Paarmann et al., 2000). It is possible, however, that AMPA mediated depolarization may trigger calcium influx indirectly through the voltage-gated calcium channel activation on the post-synaptic terminal. The latter may contribute to synaptically triggered calcium influx as blockade of P/Q-type (but not L- N-type) calcium channels reduces XII motor output driven from the pre-BötC in the in vitro mouse slice preparations from normal animals (Lieske and Ramirez, 2006a; Koch et al., 2013).
Calcium permeability through KAR receptors is dependent on subunit expression. The KAR subunit GluK3 is highly expressed in the pre-BötC (Paarmann et al., 2000) and is calcium permeable (Perrais et al., 2009) making it a possible candidate for synaptically mediated calcium entry. Furthermore, GluK3 is insensitive to tonic glutamate release and only activated by large glutamate transients (Perrais et al., 2009). Consequently, GluK3 may only be activated when receiving synaptic input from a bursting presynaptic neuron which would presumably generate large glutamate transients. The role of GluK3 in the pre-BötC has not been investigated.
Metabotropic glutamate receptors (mGluR) indirectly activate ion channels through G-protein mediated signaling cascades. Group 1 mGluRs which include mGluR1 and mGluR5 are typically located on post-synaptic terminals (Shigemoto et al., 1997) and activation of group 1 mGluRs is commonly associated with calcium influx through calcium permeable channels (Berg et al., 2007; Endoh, 2004; Mironov, 2008) and calcium release from intracellular calcium stores (Pace et al., 2007).
In the pre-BötC, mGluR1/5 are thought to contribute to calcium influx by triggering the release of calcium from intracellular stores (Pace et al., 2007) and/or the activation of the transient receptor potential C3 (TRPC3) channel (Ben-Mabrouk and Tryba, 2010). Blockade of mGluR1/5 reduces the inspiratory drive potential in pre-BötC neurons (Pace et al., 2007) without significant perturbation of inspiratory frequency (Pace et al., 2007; Lieske and Ramirez, 2006a), which is consistent with the effects of /TRPM4 blockade (Koizumi et al., 2018). TRPC3 is a calcium permeable channel (Thebault et al., 2005) that is associated with calcium signaling (Hartmann et al., 2011), store-operated calcium entry (Kwan et al., 2004), and synaptic transmission (Hartmann et al., 2011). TRPC3 is activated by diacylglycerol (DAG) (Clapham, 2003), which is formed after synaptic activation of mGluR1/5. TRPC3, which is highly expressed in pre-BötC glutamatergic inspiratory neurons, is often co-expressed with TRPM4 (Koizumi et al., 2018) and was hypothesized to underlie activation in the pre-BötC (Ben-Mabrouk and Tryba, 2010) and other brain regions (Amaral and Pozzo-Miller, 2007; Zitt et al., 1997). Furthermore, TRPC3 and have been shown to underlie slow excitatory post-synaptic current (sEPSC) (Hartmann et al., 2008; Hartmann et al., 2011). This is consistent with our model since activation is dependent on synaptically triggered calcium entry, and the calcium dynamics are slower than the fast AMPA based current . Therefore, in our model, decays relatively slowly and, hence, can be treated as a sEPSC.
The effect of TRPC3 blockade by 3-pyrazole on pre-BötC network activity amplitude is remarkably similar to that with blockade of TRPM4 (Koizumi et al., 2018). This suggests that the /TRPM4 activation may be dependent on/coupled to TRPC3. A possible explanation is that TRPC3, which is calcium permeable, mediates synaptically triggered calcium entry. It is also likely that TRPC3 plays a role in maintaining background calcium concentration levels. We tested this hypothesis by simulating the blockade of synaptically-triggered calcium influx while simultaneously lowering the background calcium concentration (Figure 11). These simulations generated large reductions in activity amplitude with no effect on frequency which are consistent with data from experiments where TRPC3 is blocked using 3-pyrazole (Koizumi et al., 2018). This indirectly suggests that TRPC3 may be critical for synaptically-triggered calcium entry and subsequent activation.
-dependent rhythmogenic kernel
is a conductance present ubiquitously in pre-BötC inspiratory neurons, and is established to underlie intrinsic oscillatory neuronal bursting in the absence of excitatory synaptic interactions in neurons with sufficiently high conductance densities (Del Negro et al., 2002; Koizumi and Smith, 2008; Koizumi and Smith, 2008; Yamanishi et al., 2018). Accordingly, we randomly incorporated this conductance in our model excitatory neurons from a uniform statistical distribution to produce heterogeneity in conductance density across the population. Our simulations indicate that the circuit neurons mostly with relatively high conductance values underlie rhythm generation and remain active after compete blockade of in our model network, thus forming a -dependent rhythmogenic kernel, including some neurons with intrinsic oscillatory bursting behavior when synaptically uncoupled.
As noted above, the rhythmogenic properties of individual neurons depend on whether their conductance is high enough and, therefore, the number of intrinsic bursters in the model is defined by the width of this conductance distribution over the population. However, the critical value of conductance and intrinsic bursting properties in general were also shown theoretically and experimentally to critically depend on the conductances of other ionic channels (e.g. leak and delayed rectifier potassium conductances) and extracellular ion concentrations such as potassium concentration (Bacak et al., 2016b; Koizumi and Smith, 2008; Rybak et al., 2003). Therefore, we believe that in reality the specific composition of the rhythmogenic kernel and its oscillatory capabilities strongly depend on the existing combination of neuronal conductances and on the in vitro experimental conditions.
Recently, it has become apparent that there is functional heterogeneity within pre-BötC excitatory circuits, including distinct subpopulations of neurons involved in generating periodic sighs (Toporikova et al., 2015; Li et al., 2016), arousal (Yackle et al., 2017), and the subpopulations generating regular inspiratory activity. Activity of the normal inspiratory and sigh-generating subpopulations in the pre-BötC isolated in vitro is proposed to be dependent on activation of (Toporikova et al., 2015). Our experimental and modeling results suggest that within the normal inspiratory population, there are subpopulations distinguished by their role in rhythm versus amplitude generation due to biophysical properties: there is a /TRPM4-dependent recruitable population of excitatory neurons for burst amplitude generation and the -dependent rhythmogenic kernel population. The spatial arrangements of these two synaptically interconnected excitatory populations within the pre-BötC are currently unknown, and it remains an important experimental problem to identify the cells constituting the rhythmogenic kernel and their biophysical properties. This should now be possible, since our analysis and experimental results suggest that the rhythmically active neurons of the kernel population can be revealed and studied after pharmacologically inhibiting the /TRPM4-dependent inspiratory burst-generating population.
A 'burstlet theory' for emergent network rhythms has recently been proposed to account for inspiratory rhythm and pattern generation in the isolated pre-BötC in vitro (Kam et al., 2013; Del Negro et al., 2018). This theory postulates that a subpopulation of excitatory neurons generating small amplitude oscillations (burstlets) functions as the inspiratory rhythm generator that drives neurons that generate the larger amplitude, synchronized inspiratory population bursts. This concept emphasizes that subthreshold neuronal membrane oscillations need to be considered and that there is a neuronal subpopulation that functions to independently form the main inspiratory bursts. This is similar to our concept of distinct excitatory subpopulations generating the rhythm versus the amplitude of inspiratory oscillations. Biophysical mechanisms generating rhythmic burstlets and the large amplitude inspiratory population bursts in the burstlet theory are unknown, and the general problem of understanding the dynamic interplay of circuit interactions and cellular biophysical processes in the generation of population-level bursting activity has been highlighted (Richter and Smith, 2014; Ramirez and Baertsch, 2018b). We have identified a major Ca2+-dependent conductance mechanism for inspiratory burst amplitude (pattern) generation and show theoretically how this mechanism may be coupled to excitatory synaptic interactions and is independent of the rhythm-generating mechanism. We also note that a basic property of is its ability to generate subthreshold oscillations and promote burst synchronization (Butera et al., 1999b; Bacak et al., 2016a). However, in contrast to our proposal for the mechanisms operating in the kernel rhythm-generating subpopulation, with its favorable voltage-dependent and kinetic autorhythmic properties– is not proposed to be a basic biophysical mechanism for rhythm generation in the burstlet theory (Del Negro et al., 2018).
We emphasize that the above discussions regarding the role of pertain to the excitatory circuits in the isolated pre-BötC including in more mature rodent experimental preparations in situ where inspiratory rhythm generation has also been shown to be dependent on (Smith et al., 2007). The analysis is more complex when the pre-BötC is embedded within interacting respiratory circuits in the intact nervous system generating the full complement of inspiratory and expiratory phase activity (Lindsey et al., 2012; Ramirez and Baertsch, 2018a; Richter and Smith, 2014), where rhythmogenesis is tightly controlled by inhibitory circuit interactions, including via local inhibitory circuits in the pre-BötC (Harris et al., 2017; Ausborn et al., 2018; Baertsch et al., 2018; Ramirez and Baertsch, 2018a), and the contribution of kinetic properties alone in setting the timing of inspiratory oscillations is diminished (Smith et al., 2007; Richter and Smith, 2014; Rubin et al., 2009b). Extending our analysis to consider inhibitory circuit interactions with the excitatory subpopulations generating oscillation frequency and amplitude that we propose will provide additional insight into biophysical mechanisms controlling these two processes.
Conclusions
Based on our computational model, distinct biophysical mechanisms are involved in generating the rhythm and amplitude of inspiratory oscillations in the isolated pre-BötC excitatory circuits. According to this model, inspiratory rhythm generation arises from a group of -dependent excitatory neurons, including cells with intrinsic oscillatory bursting properties, that form a rhythmogenic kernel. Rhythmic synaptic drive from these neurons triggers post-synaptic calcium transients, activation, and subsequent membrane depolarization which drives rhythmic bursting in the rest of the population of inspiratory neurons. We showed that activation of by synaptically-driven calcium influx functions as a mechanism that amplifies the excitatory synaptic input to generate the inspiratory drive potential and population activity amplitude in these non-rhythmogenic neurons. Consequently, reduction of causes a robust decrease in overall network activity amplitude via de-recruitment of these burst amplitude-generating neurons without substantial perturbations of the inspiratory rhythm. Thus, plays a critical role in generating the amplitude of rhythmic population activity, which is consistent with the results from experimental inhibition of /TRPM4 channels (Peña et al., 2004; Koizumi et al., 2018). Our model provides a theoretical explanation for these experimental results and new insights into the biophysical operation of pre-BötC excitatory circuits. The theoretical framework that we have developed here should provide the bases for further exploration of biophysical mechanisms operating in the mammalian respiratory oscillator.
Materials and methods
Model description
Request a detailed protocolThe model describes a network of synaptically coupled excitatory neurons. Simulated neurons are comprised of a single compartment described using a Hodgkin Huxley formalism. For each neuron, the membrane potential is given by the following current balance equation:
where is the membrane capacitance, , , , , , and are ionic currents through sodium, potassium, leak, persistent sodium, calcium activated non-selective cation, voltage-gated calcium, and synaptic channels, respectively. Description of these currents, synaptic interactions, and parameter values are taken from Jasinski et al. (2013). The channel currents are defined as follows:
where is the maximum conductance, is the reversal potential, and are voltage dependent gating variables for channel activation and inactivation, respectively, and . The parameters and are given in Table 1.
For , , and , the dynamics of voltage-dependent gating variables , and are defined by the following differential equation:
where steady state activation/inactivation and time constant are given by:
For the voltage-gated potassium channel, steady state activation and time constant are given by:
where
The parameters , , , , , , , , , and are given in Table 1. activation is dependent on the intracellular calcium concentration and is given by:
The parameters and n, given in Table 1, represent the half-activation calcium concentration and the Hill Coefficient, respectively.
Calcium enters the neurons through voltage-gated calcium channels () and/or synaptic channels (), where a percentage () of the synaptic current () is assumed to consist of Ca2+ ions. A calcium pump removes excess calcium with a time constant and sets the minimum calcium concentration . The dynamics of is given by the following differential equation:
The parameters is a conversion factor relating current and rate of change in , see Table 1 for parameter values.
The synaptic conductance of the ith neuron () in the population is described by the following equation:
where is the weight of the synaptic connection from cell to cell , C is a connectivity matrix ( if neuron makes a synapse on neuron , and otherwise), is the Heaviside step function, is time, is the exponential decay constant and is the time at which an action potential is generated in neuron and reaches neuron .
To account for heterogeneity of neuron properties within the network, the persistent sodium current conductance, , for each neuron was assigned randomly based on a uniform distribution over the range which is consistent with experimental measurements (Rybak et al., 2003; Koizumi and Smith, 2008; Koizumi and Smith, 2008). We also uniformly distributed over the range , however, simulation results did not depend on whether we used such a distribution, or assigned for all neurons to the same value of , which is the mean of this distribution. In simulations where was varied, we multiplied for each neuron by the same factor. This factor was used as a control parameter for all such simulations, and shown as a percentage of the baseline or as the mean values for the population in figures. The weight of each synaptic connection was uniformly distributed over the range where ranged from to depending on the network connectivity and specific simulation. The elements of the network connectivity matrix, , are randomly assigned values of or such that the probability of any connection between neuron and neuron being is equal to the network connection probability . We varied the connection probability over the range , however, a value of was used in most simulations.
Data analysis and definitions
Request a detailed protocolThe time of an action potential was defined as when the membrane potential of a neuron crosses in a positive direction. The network activity amplitude and frequency were determined by identifying peaks and calculating the inverse of the interpeak interval in histograms of network spiking. Network histograms of the population activity were calculated as the number of action potentials generated by all neurons per 50 ms bin per neuron with units of . The number of recruited neurons is defined as the peak number of neurons that spiked at least once per bin during a network burst. The average spike frequency of recruited neurons is defined as the number of action potentials per bin per recruited neuron with units of . The average network resting membrane potential was defined as the average minimum value of in a window following a network burst. The average inactivation of the persistent sodium current at the start of each burst was defined by the maximum of the average value of in a 500-ms window before the peak of each network burst. The average inactivation of the persistent sodium current at the end of each burst was defined by the maximum of the average value of in a 500-ms window after the peak of each network burst. Synaptic strength is defined as the number of neurons in the network multiplied by the connection probability multiplied by the average weight of synaptic connections (). Pacemaker neurons were defined as neurons that continue bursting intrinsically after complete synaptic blockade. Follower neurons were defined as neurons that become silent after complete synaptic blockade. The inspiratory drive potential is defined as the envelope of depolarization that occurs in neurons during the inspiratory phase of the network oscillations (Morgado-Valle et al., 2008).
Characterization of in regulating network activity amplitude and frequency in and models
Request a detailed protocolTo characterize the role of in regulation of network activity amplitude and frequency we slowly increased the conductance () in our simulations from zero until the network transitioned from a rhythmic bursting to a tonic (non-bursting) firing regime. To ensure that the effect(s) are robust, these simulations were repeated over a wide range of synaptic weights, synaptic connection probabilities, and strengths of the intracellular calcium transients from or sources. Changes in network activity amplitude were further examined by plotting the number of recruited neurons and the average action potential frequency of recruited neurons versus .
Simulated pharmacological manipulations
Request a detailed protocolIn simulations that are compared with experimental data, both and calcium sources are included. Pharmacological blockade of was simulated by varying the conductance, according to a decaying exponential function
.
The percent block , decay constant and the maximum conductance were adjusted to match the experimental changes in network amplitude. The synaptic weight of the network was chosen such that at the network activity amplitude was close to 20% of maximum. To reduce the computational time, the duration of block simulations was one tenth of the total of experimental durations. For comparison, the plots of normalized change in amplitude and frequency of the simulations were stretched over the same time-period as experimental data. Increasing the simulation time had no effect on our results (data not shown).
Comparison with calcium imaging data
Request a detailed protocolTo allow comparisons with network and cellular calcium imaging data, we analyzed rhythmic calcium transients from our simulations. Single cell calcium signals are represented by . The network calcium signal was calculated as the average intracellular calcium concentration in the network ().
Integration methods
Request a detailed protocolAll simulations were performed locally on an 8-core Linux-based operating system or on the high-performance computing cluster Biowulf at the National Institutes of Health. Simulation software was custom written in C++. Numerical integration was performed using the exponential Euler method with a fixed step-size () of . In all simulations, the first 50 s of simulation time was discarded to allow for the decay of any initial condition-dependent transients.
Data availability
All data in this study are generated by computational simulations. All model parameters and equations are included in the manuscript and source code is included with this submission.
References
-
Organization of the core respiratory network: insights from optogenetic and modeling studiesPLOS Computational Biology 14:e1006148.https://doi.org/10.1371/journal.pcbi.1006148
-
Modeling the effects of extracellular potassium on bursting properties in pre-Bötzinger complex neuronsJournal of Computational Neuroscience 40:231–245.https://doi.org/10.1007/s10827-016-0594-8
-
Substance P modulation of TRPC3/7 channels improves respiratory rhythm regularity and ICAN-dependent pacemaker activityEuropean Journal of Neuroscience 31:1219–1232.https://doi.org/10.1111/j.1460-9568.2010.07156.x
-
Models of respiratory rhythm generation in the pre-Bötzinger complex. I. bursting pacemaker neuronsJournal of Neurophysiology 82:382–397.https://doi.org/10.1152/jn.1999.82.1.382
-
BookRhythms of the BrainOxford University Press.https://doi.org/10.1093/acprof:oso/9780195301069.001.0001
-
Patterns of inspiratory phase-dependent activity in the in vitro respiratory networkJournal of Neurophysiology 109:285–295.https://doi.org/10.1152/jn.00619.2012
-
Cycle-by-cycle assembly of respiratory network activity is dynamic and stochasticJournal of Neurophysiology 109:296–305.https://doi.org/10.1152/jn.00830.2011
-
Sodium and calcium current-mediated pacemaker neurons and respiratory rhythm generationJournal of Neuroscience 25:446–453.https://doi.org/10.1523/JNEUROSCI.2237-04.2005
-
Looking for inspiration: new perspectives on respiratory rhythmNature Reviews Neuroscience 7:232–241.https://doi.org/10.1038/nrn1871
-
Different roles for inhibition in the rhythm-generating respiratory networkJournal of Neurophysiology 118:2070–2088.https://doi.org/10.1152/jn.00174.2017
-
mGluR1/TRPC3-mediated synaptic transmission and calcium signaling in mammalian central neuronsCold Spring Harbor Perspectives in Biology 3:a006726.https://doi.org/10.1101/cshperspect.a006726
-
Pacemaker behavior of respiratory neurons in medullary slices from neonatal ratJournal of Neurophysiology 72:2598–2608.https://doi.org/10.1152/jn.1994.72.6.2598
-
Locomotor circuits in the mammalian spinal cordAnnual Review of Neuroscience 29:279–306.https://doi.org/10.1146/annurev.neuro.29.051605.112910
-
Stable respiratory activity requires both P/Q-type and N-type voltage-gated calcium channelsJournal of Neuroscience 33:3633–3645.https://doi.org/10.1523/JNEUROSCI.6390-11.2013
-
Pattern-specific synaptic mechanisms in a multifunctional network. I. Effects of alterations in synapse strengthJournal of Neurophysiology 95:1323–1333.https://doi.org/10.1152/jn.00505.2004
-
Computational models and emergent properties of respiratory neural networksComprehensive Physiology 2:1619–1670.https://doi.org/10.1002/cphy.c110016
-
Principles of rhythmic motor pattern generationPhysiological Reviews 76:687–717.https://doi.org/10.1152/physrev.1996.76.3.687
-
Somatic Ca2+ transients do not contribute to inspiratory drive in preBötzinger Complex neuronsThe Journal of Physiology 586:4531–4540.https://doi.org/10.1113/jphysiol.2008.154765
-
Ionic currents and endogenous rhythm generation in the pre-Bötzinger complex: modelling and in vitro studiesAdvances in Experimental Medicine and Biology 551:121–126.https://doi.org/10.1007/0-387-27023-X_19
-
Pacemaker neurons and neuronal networks: an integrative viewCurrent Opinion in Neurobiology 14:665–674.https://doi.org/10.1016/j.conb.2004.10.011
-
The dynamic basis of respiratory rhythm generation: one breath at a timeAnnual Review of Neuroscience 41:475–499.https://doi.org/10.1146/annurev-neuro-080317-061756
-
Multiple rhythmic states in a model of the respiratory central pattern generatorJournal of Neurophysiology 101:2146–2165.https://doi.org/10.1152/jn.90958.2008
-
Sodium currents in neurons from the rostroventrolateral medulla of the ratJournal of Neurophysiology 90:1635–1642.https://doi.org/10.1152/jn.00150.2003
-
Rhythmic bursting in the pre-Bötzinger complex: mechanisms and modelsProgress in Brain Research 209:1–23.https://doi.org/10.1016/B978-0-444-63274-6.00001-1
-
Differential presynaptic localization of metabotropic glutamate receptor subtypes in the rat hippocampusThe Journal of Neuroscience 17:7503–7522.https://doi.org/10.1523/JNEUROSCI.17-19-07503.1997
-
Brainstem respiratory networks: building blocks and microcircuitsTrends in Neurosciences 36:152–162.https://doi.org/10.1016/j.tins.2012.11.004
-
Receptor-operated Ca2+ entry mediated by TRPC3/TRPC6 proteins in rat prostate smooth muscle (PS1) cell lineJournal of Cellular Physiology 204:320–328.https://doi.org/10.1002/jcp.20301
-
Two types of independent bursting mechanisms in inspiratory neurons: an integrative modelJournal of Computational Neuroscience 30:515–528.https://doi.org/10.1007/s10827-010-0274-z
-
Persistent sodium current kinetics orchestrate oscillatory burst firing in the pre-Bötzinger complexThe Journal of General Physiology 150:1523–1540.https://doi.org/10.1085/jgp.201812100
-
Expression of TRPC3 in chinese hamster ovary cells results in calcium-activated cation currents not related to store depletionThe Journal of Cell Biology 138:1333–1341.https://doi.org/10.1083/jcb.138.6.1333
Article and author information
Author details
Funding
The Jayne Koskinas Ted Giovanis Foundation for Health and Policy
- Ryan S Phillips
National Institute of Neurological Disorders and Stroke (Intramural Research Program of the National Institutes of Health)
- Ryan S Phillips
- Tibin T John
- Hidehiko Koizumi
- Jeffrey C Smith
National Institutes of Health (R01 AT008632)
- Yaroslav I Molkov
National Institutes of Health (U01 EB021960)
- Yaroslav I Molkov
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
This work was supported in part by the Jayne Koskinas and Ted Giovanis Foundation for Health and Policy, the Intramural Research Program of the National Institutes of Health (NIH), National Institute of Neurological Disorders and Stroke (NINDS), and NIH Grants R01 AT008632 and U01 EB021960.
Copyright
© 2019, Phillips et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 1,266
- views
-
- 217
- downloads
-
- 25
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Computational and Systems Biology
- Neuroscience
Computational models are helping researchers to understand how certain properties of neurons contribute to respiratory rhythms.
-
- Computational and Systems Biology
Degree distributions in protein-protein interaction (PPI) networks are believed to follow a power law (PL). However, technical and study bias affect the experimental procedures for detecting PPIs. For instance, cancer-associated proteins have received disproportional attention. Moreover, bait proteins in large-scale experiments tend to have many false-positive interaction partners. Studying the degree distributions of thousands of PPI networks of controlled provenance, we address the question if PL distributions in observed PPI networks could be explained by these biases alone. Our findings are supported by mathematical models and extensive simulations and indicate that study bias and technical bias suffice to produce the observed PL distribution. It is, hence, problematic to derive hypotheses about the topology of the true biological interactome from the PL distributions in observed PPI networks. Our study casts doubt on the use of the PL property of biological networks as a modeling assumption or quality criterion in network biology.