Hypothesis and Theory ARTICLE
Anesthetic action on extra-synaptic receptors: effects in neural population models of EEG activity
- 1INRIA CR Nancy - Grand Est, Team Neurosys, Villers-les-Nancy, France
- 2Department of Anaesthesiology, Waikato Clinical School, University of Auckland, Hamilton, New Zealand
The role of extra-synaptic receptors in the regulation of excitation and inhibition in the brain has attracted increasing attention. Because activity in the extra-synaptic receptors plays a role in regulating the level of excitation and inhibition in the brain, they may be important in determining the level of consciousness. This paper reviews briefly the literature on extra-synaptic GABA and NMDA receptors and their affinity to anesthetic drugs. We propose a neural population model that illustrates how the effect of the anesthetic drug propofol on GABAergic extra-synaptic receptors results in changes in neural population activity and the electroencephalogram (EEG). Our results show that increased tonic inhibition in inhibitory cortical neurons cause a dramatic increase in the power of both δ− and α− bands. Conversely, the effects of increased tonic inhibition in cortical excitatory neurons and thalamic relay neurons have the opposite effect and decrease the power in these bands. The increased δ-activity is in accord with observed data for deepening propofol anesthesia; but is absolutely dependent on the inclusion of extrasynaptic (tonic) GABA action in the model.
General anesthesia is used daily to enable surgery, but its underlying mechanisms of action are still largely a mystery. In recent decades there have been successful efforts to reveal the drug action on single receptors (Franks and Lieb, 1994; Alkire et al., 2008; Brickley and Mody, 2012), however their effect on neural populations, networks of neural populations, and brain areas, still remains unsolved. To explain the underlying neural mechanism during the loss of consciousness, two prominent hypotheses are the loss of integration information, developed by Tononi (Tononi, 2004; Murphy et al., 2011; Boly et al., 2012), and a sharp phase transition of the brain activity involving a drop of neural activity, put forward by Steyn-Ross et al. (Steyn-Ross et al., 2004; Friedman et al., 2010). These hypotheses are not mutually exclusive. For instance, a recent experimental study on the effects of propofol on neural activity measured at various spatial scales (Lewis et al., 2012) has revealed both decreased functional connectivity between brain areas and a dramatic drop of neuron firing rates after loss of consciousness. A large amount of experimental literature has revealed characteristic spectral signal changes in electroencephalographic data (EEG) and Local Field Potentials (LFPs) during general anesthesia (Cimenser et al., 2011; Lewis et al., 2012; Purdon et al., 2012; Sellers et al., 2013; Vizuete et al., 2014). Moreover, several previous theoretical studies have proposed neural models to explain certain EEG signal features observed during anesthesia (Steyn-Ross et al., 1999, 2013; Bojak and Liley, 2005; Wilson et al., 2006; Foster et al., 2008; McCarthy et al., 2008; Ching et al., 2010, 2012; Hutt, 2013; Liley and Walsh, 2013; Hutt et al., 2013; Hutt and Buhry, 2014). Although these studies may incorporate realistic neurobiological details of the brains' network topology and neuronal function, they have simplified dramatically the anesthetic action by considering only synaptic excitatory and inhibitory receptors. There is a growing amount of experimental research that has revealed the importance of extra-synaptic receptors (ESR) for neural interactions in general (Brickley and Mody, 2012; Hardingham and Bading, 2012), and for anesthesia especially, see Alkire et al. (2008); Hutt (2012) and references therein.
To elucidate the role of ESR in the context of anesthesia, one approach might be to do a theoretical study of a realistic neural population model which reproduces the characteristic signal features observed in EEG. To perform such a theoretical study, it is necessary to incorporate physiological properties of extrasynaptic receptors into neural population models.
Gamma-aminobutyric acid (GABA) receptors are a large and important class of inotropic receptors. These receptors are located in the neuron's membrane and respond to the neurotransmitter GABA by opening Cl− channels and inducing an inward hyperpolarizing membrane current. This response may either be: phasic at synaptic receptors or, tonic at ESR which lie distant from synaptic terminals (Kaneda et al., 1995; Brickley et al., 1996; Semyanov et al., 2003, 2004; Yeung et al., 2003; Belelli et al., 2009). The phasic response evolves on a time scale of 10–200 ms whereas tonic response evolves on a much longer time scale (Hamann et al., 2002; Cavalier et al., 2005).
The precise biochemical origin of tonic inhibition is still heavily debated (Farrant and Nusser, 2005; Bright et al., 2011). A rather simple and intuitive model explains the tonic current as a spillover of excess neurotransmitter from synapses. This is due to incomplete GABA uptake by nearby synaptic GABAA-receptors. The remaining neurotransmitter is thus able to diffuse to more distant GABAA-receptors via extracellular space (Nusser et al., 1997; Semyanov et al., 2004; Farrant and Nusser, 2005; Bright et al., 2007). This spillover may explain the longer time scale of tonic responses found experimentally. In addition, this explanation implies that even small concentrations of neurotransmitters are sufficient to generate tonic activity because of the high sensitivity of ESRs.
The effect of ESRs on the dendritic activity has not attracted much attention. This may because there are only a relatively small number of such receptors as compared to synaptic receptors (Kopanitsa, 1997; Farrant and Nusser, 2005). Moreover, only recently have experimental studies been able to classify and localize different sub-types of GABAA receptors (Semyanov et al., 2004; Farrant and Nusser, 2005). GABAA receptors are pentameric ligand-gated ion channels and it has been found that δ−sub units of GABAA receptors occur exclusively at ESRs (Nusser et al., 1998; Wei et al., 2003; Farrant and Nusser, 2005; Belelli et al., 2009; Ye et al., 2013). This indicates a specific role of these receptors for the neural information processing in general with specific implications in diseases (Brickley and Mody, 2012) and consciousness (Kopanitsa, 1997).
Tonic inhibition induced by extra-synaptic GABAA-receptors represents a persistent increase in the cell membrane's conductance. On the single neuron level, this diminishes the membrane time constant and, consequently, reduces the size and duration of excitatory post-synaptic potentials propagating on the dendrite. Hence tonic inhibition reduces the excitability of the membrane and increases the effective firing threshold (Farrant and Nusser, 2005). At the neural population level, ESRs affect the excitability of interneuron-pyramidal cell networks and thus modify network oscillations (Semyanov et al., 2003). Kopanitsa (1997) argues that the sustained spatially widespread tonic inhibition is energetically more effective for the system to diminish neural population activity than short-lasting local phasic inhibition, since lower neurotransmitter concentrations are sufficient. The critical factor in this mechanism is the the relatively high sensitivity of ESRs to modulations by anesthetic agents (Yeung et al., 2003; Farrant and Nusser, 2005; Orser, 2006; Houston et al., 2012). The brain areas that have been shown to be affected by anesthetic-induced tonic inhibition are the hippocampus (Bai et al., 2001), brain stem (McDougall et al., 2008), cerebellum (Houston et al., 2012), and the thalamus (Belelli et al., 2009). Since these areas are supposed to play a role in general anesthesia (Alkire et al., 2008), ESRs may mediate clinical anesthetic effects, such as hypnosis and amnesia (Kretschmannova et al., 2013). Thus, it is reasonable to argue that GABAA ESRs set the background inhibition of neural populations and the brain network and mediate slow consciousness phenomena, such as loss of consciousness, sleep or arousal (Kopanitsa, 1997).
Converse to GABAergic receptors, N-methyl-D-aspartate (NMDA) receptors respond to the neurotransmitter glutamate by excitatory inward Na+ and Ca2+ currents and K+ outward currents. The response of NMDA receptors to glutamate depends on their spatial location with respect to synaptic terminals and the presence of co-agonists. A recent experimental study has revealed that the population of NMDA receptors, which are close to synaptic terminals, are primarily activated by the co-agonist d-serine in the presence of glutamate; while extra-synaptic NMDA receptors (more distant from the synaptic terminals) respond primarily to glutamate and the co-agonist glycine (Mothet et al., 2000; Papouin et al., 2012). D-serine and glycine are endogenous amino acids found naturally in the brain (d-serine is a derivative of glycine). Similar to GABAergic ESRs, it has been shown that there exists a significant ambient glutamate concentration which induces a tonic excitatory current (Sah et al., 1989; Fleming et al., 2011). This current is evoked primarily at extra-synaptic NMDA receptors (Le Meur et al., 2007) and may be regulated by other cells, such as neighboring astrocytes (Panatier et al., 2006; Fleming et al., 2011) which control glutamate uptake and also synthesize d-serine (Wolosker et al., 1999).
Commonly-used GABAergic anesthetic drugs directly modify the corresponding receptors. However, various anesthetics are also known to affect the endogenous co-agonists of NMDA receptors (Martin et al., 1995; Daniels and Roberts, 1998; Papouin et al., 2012). Hence, the possible anesthetic effect on NMDA receptors is more complex and indirect than for GABAergic ESRs. There is a large class of NMDA receptor antagonists, that inhibit the excitatory action of NMDA receptors. These anesthetics induce so-called dissociative anesthesia (Pender, 1970) leading to amnesia and analgesia without depressing respiration, but also characterized by distorted perceptions of sight and sound and feelings of dissociation from the environment. An example of a dissociative anesthetic drug is the inhalational anesthetic xenon which—amongst other actions—binds primarily to the extra-synaptic glycine site of NMDA receptors (Dickinson et al., 2007) and attenuates long-term potentiation present in the hippocampus by reducing extrasynaptic receptor currents (Kratzer et al., 2012).
To understand how the anesthetic effect of ESR activity on the microscopic single neuron scale could lead to changes in EEG and behavior that can be observed at macroscopic scales, it is necessary to establish a bridge between the two scales. This bridge may be formulated as a dynamical theoretical model. Neural population models represent a good candidate for a dynamic description of neural activity at an intermediate mesoscopic scale (Coombes, 2006; Nunez and Srinivasan, 2006; Bressloff, 2012). These models describe properties of ensembles of neurons, such as the mean firing rate and the mean dendritic current (Hutt, 2009), whilst their output variables can be strongly linked to macroscopic experimental quantities such as Local Field Potentials (LFPs) and EEG (Wright and Kydd, 1992; Nunez and Srinivasan, 2006). An increasing number of theoretical studies have used neural population models to describe signal features in LFPs and EEG observed during anesthesia (Foster et al., 2008; Hutt et al., 2013). Most of these studies take into account anesthetic action on excitatory and/or inhibitory synapses (Steyn-Ross et al., 2004; Liley and Bojak, 2005; Hutt and Longtin, 2009; Ching et al., 2010; Hindriks and van Putten, 2012) while few consider ESRs (Talavera et al., 2009). This link between the synaptic receptor properties in an ensemble of neurons and the average population dynamics is straight-forward, since classical neural population models already involve the average synaptic response function. The situation is different for ESRs, since their action is not incorporated into the classical models. A very recent work has filled this gap (Hutt and Buhry, 2014). This theoretical work demonstrated a method to include mathematically extra-synaptic GABAA receptor action in neural population models; which enables researchers to study how changing anesthetic ESR action modifies spectral features in the EEG, which might then be observed experimentally.
The current work uses a thalamo-cortical neural population model involving anesthetic synaptic inhibition with a well-established connection topology; and then extends this model by including the effects of extra-synaptic GABAergic receptor action in the presence of the anesthetic drug propofol. With the help of this model, we demonstrate the role of extra-synaptic GABAergic inhibition, and the importance of tonic inhibition in the cortical inhibitory neuronal population, in explaining experimental EEG power spectra.
2. Materials and Methods
2.1. EEG Data
We re-analyzed previously-obtained experimental data from subjects that had been given a short propofol anesthetic. The details of the methods can be found in Johnson et al. (2003). In brief, after obtaining regional ethical committee approval and written informed consent, five healthy subjects (mean age 27.7 yrs, four males) were studied. They were on no psychoactive drugs and had been starved for at least 6 h prior to the study. They were monitored and managed as per clinical anesthesia, according to the Australia and New Zealand College of Anesthetists best practice guidelines. The induction consisted of an intravenous infusion of propofol at 1500 mg/hr until the subject no longer responded to verbal command. Typically this occurred about 5 min into the infusion. The estimated effect-site concentrations of propofol were calculated using standard population-based pharmacokinetics models.
The EEG was acquired using the Electrical Geodesics 128 channel Ag/AgCl electrode system (Eugene, CO, USA) referenced to Cz. Electrode impedances were below 30 KOhm (100 MOhm input impedance amplifier). The sampling frequency was 250 Hz, with a 0.1–100 Hz analog band pass filter, and A–D conversion was at 12 bits precision. The EEG data were re-referenced to a grand mean, and band-pass filtered using 3-rd order Butterworth filters 0.2–45 Hz to eliminate line-noise. An additional Whittaker filter was applied to reduce movement and blink artifacts. The power in each frequency was obtained applying a short-time Fourier transform with a moving window of 60 s and 54 s overlap. The power spectra have been computed 1 min before infusion start (t = 1 min) and 4 min after infusion (t = 5 min). For visualization reasons, these power spectra at different time instances have been smoothed by a running average over frequencies with a 1 Hz window and a 0.017 Hz frequency step.
2.2. Thalamo-Cortical Model
The body of the model (Robinson et al., 2001; Rennie et al., 2002) is based on a population-level description of a single thalamo-cortical module comprising four populations of neurons, namely excitatory (E) and inhibitory (I) cortical population, a population built of thalamo-cortical relay neurons (S) and of thalamic reticular neurons (R), as shown in Figure 1. The details of the model and the nominal parameter values are taken from a previous work (Robinson et al., 2001, 2002). This model is based on the original idea of Lopes da Silva et al., stating that the α-rhythm represents the noisy thalamic input signal band-pass filtered by feedback-connected cortical and thalamic neural populations (Lopes da Silva et al., 1974). Here we just briefly describe the key concepts of the model. The average soma membrane potential denoted by Va, for a = E, I, S, R is modeled by
where ⊗ represents the temporal convolution, hb(t) = Hbb(t) where b(t) denotes the mean synaptic response function defined by
and α and βb (with units Hz) are the synaptic rise and decay rates of the synaptic response, respectively. The synaptic decay rates and synaptic response functions depend on the source neurons of type b only and are independent of the target neurons. The constant pre-factor Hb defines the response function amplitude subject to the anesthetic concentration. Here, we assume identical excitatory synaptic receptors with constant rise and decay rate. Inhibitory synaptic receptors are also assumed to exhibit identical constant rise and decay rates while their decay rates depend on the anesthetic concentration. This strong approximation is taken from a previous study (Hindriks and van Putten, 2012) to be able to compare our results, while preliminary studies with more realistic parameters show similar results (not shown).
Figure 1. The schematic of a thalamo-cortical module. The blue arrows indicate excitatory connections while the red connections with filled circle ends denote inhibitory connections. The symbols E, I, S, and R denote the excitatory and inhibitory cortical neurons, thalamo-cortical relay, and thalamic reticular neurons, respectively. In addition thalamocortical and corticothalamic connections exhibit the same nonzero time delay τ.
The constants νab are the strengths of the connections from population of type b to population of type a (in mVs), and ϕb is the average firing rate of the population b (in Hz). The connections between cortex and thalamus are associated with a same nonzero time delay, τab = τ, while the delay term is assumed to be zero within the cortex and within the thalamus (Victor et al., 2011).
By virtue of long-range axonal projections of excitatory cortical neurons and by assuming the spatially-homogeneous dynamics on the cortex, the average firing rate ϕE obeys the damped oscillator equation
where the operator D is defined as
and γ is the cortical damping rate. It is assumed that the spatial spread of activity is very fast in other populations and the activity variable can be approximated by a sigmoidal function as ϕb = S(Vb), for b = I, S, R. Conversely to the original model (Robinson et al., 2001, 2002; Victor et al., 2011) we use a more realistic transfer function derived from properties of type I-neurons given by Hutt and Buhry (2014)
with σ = 10 mV and ρ = 0.08 mV−1, where σ is related to the standard deviation of firing thresholds, the parameter ρ < ∞ reflects the properties of type I-neurons, Smax is the maximum population firing rate, and θa is the mean firing threshold of neurons in population a. In contrast to the standard transfer function given in Robinson et al. (2001), the transfer function in Equation (6) is not anti-symmetric to its inflection point anymore (Hutt, 2012) and exhibits a larger non-linear gain (slope) for large potentials Va > θa compared to small potentials Va < θa. This asymmetry results from the firing properties of type-I neurons, see Hutt (2012); Hutt and Buhry (2014) for more details. For ρ → ∞, the sigmoid function becomes the conventional antisymmetric transfer function.
The external input to the system is considered as a non-specific input to thalamo-cortical relay neurons as
where 〈ϕN〉 indicates its mean value and ξ(t) is a zero average Gaussian white noise and κ is the noise intensity.
The power spectrum characterizes small fluctuations about the resting state of the system defined by dVa(t)/dt = 0. Following Robinson et al. (2001); Nunez and Srinivasan (2006), it is assumed that the activity of excitatory cortical neurons generates the EEG, and due to the specific choice of external input to thalamo-cortical relay neurons, the power spectrum of the EEG is related to the Greens function of linear deviations about the resting state by Hutt (2013)
in which PE(ω) depends just on one matrix component of the matrix Greens function (ω), see the Supplementary Material for its definition. We point out that the subsequent power spectrum analysis is based on Equation (8) and changing a system parameter, such as the factor p, changes the resting state, the corresponding non-linear gains and consequently the power spectrum. In addition, the power spectrum analysis is valid only if the resting state is stable and hence the fluctuations do not diverge. We have taken care of this additional condition and all given parameters guarantee the existence and stability of the resting state.
2.3. Effect of Propofol on Neural Populations
In order to mimic anesthetic action, we consider the general anesthetic propofol which affects synaptic and extra-synaptic GABAergic receptors. We assume that the decay rate of inhibitory synapses is identical in all neural populations under study, and decreases with increasing propofol concentration in accordance with experimental findings (Kitamura et al., 2002; Hutt and Longtin, 2009). Mathematically, such a dependence on the anesthetic concentration can be taken into account by a concentration factor p ≥ 1 and βb = β0/p while increasing p reflects an increase of the on-site concentration of propofol (Foster et al., 2008; Hindriks and van Putten, 2012; Hutt et al., 2013). Since propofol has been shown to retain the amplitude of inhibitory synaptic response functions (Kitamura et al., 2002), one can define Hb = Γ(α, β0)/Γ(α, βb) for b = I, R, where
i.e., Γ(α, βb) = b(t0) is the peak amplitude of b(t) at time t0 = ln(α/βb)/(α − βb). Thereby the maximum height of hb(t) is hmax = Γ(α, β0) which is independent of the action of propofol (Hindriks and van Putten, 2012). Moreover, since it is assumed that propofol does not act on excitatory synaptic transmission, Hb = 1 and hb(t) = b(t) for b = E, S. The GABAergic ESR tonic inhibition can be represented in the model as a constant shift of the firing threshold in neural population models (Hutt and Buhry, 2014). For simplicity, we assume a linear relationship between the anesthetic concentration parameter p and the extra-synaptic threshold shift
with the unique firing threshold θ0 = 15 mV identical for all populations in the absence of propofol and the extra-synaptic anesthetic sensitivity ka > 0. Here, (p − 1)ka is the tonic inhibition induced by extra-synaptic action which depends linearly on the propofol concentration. Future experimental studies may motivate a more realistic relationship of threshold shift and the anesthetic concentration parameter. Summarizing, synaptic and extra-synaptic inhibition, and hence anesthetic action, is present in the cortical populations E and I and in the thalamic population of relay neurons S.
2.4. Power Spectrum
The present study examines the effect of tonic inhibition in various populations E, I, S on the power spectrum of neural activity in cortical excitatory neurons, i.e., population E. We will focus on the power in the δ− and α− frequency ranges in the interval [0.5 Hz−4Hz] and [8 Hz−12 Hz], respectively.
The subsequent analysis reveals power peaks in these frequency ranges, whose magnitude changes with the level of tonic inhibition. These power peaks exhibit a maximum of power, expressed mathematically as a local maximum of the function PE(f) where PE is taken from Equation (8). The local maximum at frequency f0 is defined as dPE/df = 0, d2PE/df2 < 0 computed at f0. If there is a local maximum of power in the δ− frequency range, then δ−activity is present, whereas a missing local maximum in the δ−frequency range indicates missing δ−activity. Since the magnitude and frequency of power peaks change with the propofol concentration and extra-synaptic threshold, the concentration factor p and the extra-synaptic anesthetic sensitivity ka are the parameters of the power spectrum, i.e., PE = PE(p, ka, f).
To illustrate the usefulness of this parametrization, let us assume a factor ka0 for which no δ−power peak exists in the power spectrum PE(p, ka0, f), and ka1, ka1 > ka0 is the extra-synaptic anesthetic sensitivity leading to a spectral δ−power peak in PE(p, ka1, f) with dPE(p, ka1, fmax)/dfmax = 0 where fmax is a frequency in the δ−frequency range. Mathematically, then the continuity of all model functions and variables guarantee that there is a threshold for the emergence of δ−activity at a certain extra-synaptic anesthetic sensitivity ka, thr with ka0 ≤ ka, thr ≤ ka1. Consequently, if a threshold extra-synaptic anesthetic sensitivity for δ−activity exists, then the variation of model variables about this critical point guarantees the emergence of δ−activity. This mathematical reasoning allows us to investigate conditions under which δ−activity may emerge.
Figure 2 illustrates how the EEG power spectrum depends on the concentration of propofol for a single subject. After starting the infusion at t = 0 min, the estimated propofol effect-site concentration increases gradually with time (Figure 2A); resulting in increased power in the δ− and α−frequency ranges (Figure 2B). Over the period of the spectrogram the subject has become progressively more sedated; until a t = 5 min the subject no longer responds to verbal command but would still be responsive to nociceptive stimuli. Figure 2C shows the power spectra in the awake and sedation conditions. We observe a power enhancement primarily in the δ− and α−frequency ranges.
Figure 2. Electroencephalographic data observed under anesthesia sedation in a single subject while increasing the propofol concentration. (A) Blood plasma concentration of propofol with respect to administration time. (B) Spectrogram of frontal EEG power. The vertical lines denote time windows well before the administration (left line) and at about 5 min after the start of propofol infusion (right side); (C) Power spectra computed before the infusion of propofol (black) and 5 min after the start of infusion (red).
To understand how propofol might enhance δ− and α−power, we study the power spectrum of our theoretical model for different anesthetic concentration levels and examine the impact of adding tonic inhibition via extra-synaptic GABAA receptors. Figure 3A shows the interaction between propofol and tonic inhibition in the cortical inhibitory neuronal population. If we set the tonic inhibition to zero (kI = 0 mV), we observe a decrease in spectral power as propofol concentrations increase (i.e., the power moves from the black line to the blue line in the figure). If we set the tonic inhibition to (p − 1) · 15 mV we see the opposite effect—there is an increase of δ− and α−power (black line to red line), with increasing propofol concentration.
Figure 3. The theoretical EEG power spectrum in the baseline and in the sedation condition with and without tonic inhibition in the cortical inhibitory neurons I in (A) and the thalamic relay neurons S in (B) and corresponding simulated EEG time-seris. In (A) the administration of propofol without tonic inhibition (blue line) attenuates the power spectrum compared to the baseline condition (black line) while the tonic inhibition (red line) increases the global power and generates oscillatory activity in the δ-frequency range. In (B) increasing the anesthetic concentration yields a global power decrease in the sedation condition without tonic inhibition (blue line) and a further power decrease in the presence of tonic inhibition (red line). In (A) and (B), the black lines indicate the EEG-spectral power in the baseline condition (p = 1), and the blue and red lines show the power spectrum in anesthesia condition (p = 1.125) in the absence (ka = 0) and in the presence (ka = 15 mV) of tonic inhibition, respectively. (C) The simulated EEG time-series (ϕE(t) defined in Equation (3)) in the absence of extra-synaptic effects, i.e., kE = kI = kS = 0 mV. (D) The EEG time-series in the presence of extra-synaptic action in cortical inhibitory neurons with kI = 15 mV, kE = kS = 0 mV. The tonic inhibition changes the EEGs from low-amplitude, high-frequency pattern to high-amplitude, low-frequency pattern. In addition, the strength of cortical self-inhibition is νii = −1.8 mVs.
Previous studies have indicated that extra-synaptic inhibition in thalamic relay neurons may control the level of inhibition in the brain (Brickley and Mody, 2012). However, Figure 3B reveals that adding a nonzero tonic inhibition in the thalamic relay neurons causes a decrease in the spectral power, similar to the previous case of absent tonic inhibition in the inhibitory cortical neurons.
It is well-known that GABAergic anesthetics change the EEG from high frequency-low amplitude signals to low frequency-high amplitude signals (Gugino et al., 2001; Feshchenko et al., 2004). Figures 3C,D show simulated time series in the absence and presence of tonic inhibition in cortical inhibitory cells reproducing this experimental finding.
Our results elucidates that tonic inhibition in cortical interneurons and thalamic relay neurons affect the cortical power spectrum differently. This finding is similar to results of a previous computational neural population study of a cortico-thalamic feedback single-neuron model (Talavera et al., 2009). Figure 4 shows how the resting membrane potential (Figure 4A) and the non-linear gain (Figure 4B) in the cortical excitatory population change with differing extra-synaptic anesthetic sensitivity in cortical inhibitory neurons (kI) and in the thalamic relay neurons (kS). We observe that both the resting potential and the non-linear gain of cortical excitatory neurons increase when the cortical inhibitory extra-synaptic anesthetic sensitivity kI increases, whereas resting potential and non-linear gain of cortical excitatory neurons decrease when the extra-synaptic anesthetic sensitivity in thalamic relay neuron kS increases. Since the non-linear gain is proportional to the systems responsiveness to external stimuli, the power enhancement in population I may be explained by the augmented responsiveness of the cortical excitatory neurons. This responsiveness depends on the sub-circuit in which the neurons are involved. Since relay neurons are part of the thalamo-cortical feedback loop, while cortical inhibitory neurons contribute to the cortical loop, the cell types respond differently to the thalamic input. Essentially assuming tonic inhibition in the population of cortical excitatory neurons E, the study reveals a similar propofol concentration dependence of the power spectrum, the resting state potential and the non-linear gain as for the thalamic tonic inhibition S. This shows the unique tonic inhibition effect in the cortical inhibitory neurons.
Figure 4. Increasing the tonic inhibition (factor ka for a = I and S) affects the resting state of excitatory cortical neurons V*E (A) and the corresponding non-linear cortical gain function (B). Here the anesthetic concentration factor is identical in the populations a = E, I and S to p = 1.125. In addition, the strength of cortical self-inhibition is νii = −1.8 mVs.
Figure 3A shows the power spectrum for single values of the extra-synaptic sensitivity kI, for single values of the concentration factor p and fixed strength of cortical self-inhibition νii, while Figure 4 gives more details on the role of extra-synaptic sensitivity for fixed values of the concentration factor p and fixed cortical self-inhibition. To understand better the interplay between tonic inhibition, synaptic inhibition and the strength of cortical self-inhibition, Figure 5 shows the parameter pairs of synaptic inhibition p and the threshold of extra-synaptic sensitivity kI, thr at different self-inhibition levels, for which a peak in the δ−frequency range emerges. Recall that the kI, thr is the critical (smallest) value of extra-synaptic sensitivity in cortical inhibitory neurons kI, that lead to dPE/df = 0, d2PE/df2 < 0 computed at fmax ∈ δ−range, cf. the subsection on the power spectrum in Section 2. Parameter values beyond the respective curves lead to δ−activity power peaks. We observe that δ−activity always emerges for sufficiently strong tonic inhibition (large extra-synaptic sensitivity kI) and sufficiently strong self-inhibition νii, while the weaker the self-inhibition is the larger is the necessary extra-synaptic sensitivity or the synaptic inhibition to generate δ−activity. Even for vanishing cortical self-inhibition (νii = 0), mathematical analysis (not shown) reveals that there is still a δ−peak in the power spectrum for large enough synaptic or tonic inhibition (for p or kI large enough).
Figure 5. Parameter space for δ-power peak. The lines give the smallest (threshold) value of the extra-synaptic sensitivity kI, thr that induces δ− oscillations in the sedation condition with respect to the concentration factor p for different values of self-inhibitory connections νii. The weaker the cortical self-inhibition (the smaller |νii|), the higher the necessary level of propofol concentration (larger p) and the tonic inhibition [larger (p − 1) · kI] to induce δ−activity.
Moreover, Figure 5 reveals a minimum tonic inhibition level (minimum value of kI) beneath which no δ−power peak emerges, irrespective of the level of synaptic inhibition (p). This result indicates a major role of tonic inhibition in the generation of δ−activity, since it may support δ−activity even if the synaptic inhibition level is not sufficient to support it.
In the sedation phase, for modest concentrations of propofol, the EEG power spectrum exhibits an increase in the δ− and α−frequency ranges (Figure 2) as found experimentally in the induction phase of propofol anesthesia (San-Juan et al., 2010). One possible explanation for these phenomena is by postulating stronger GABAergic potentiation within cortical inhibitory neurons than within cortical pyramidal neurons (Hindriks and van Putten, 2012). We hypothesize, that cortical GABAergic self-inhibition plays a decisive role. Figure 3 reveals that the power surge in these frequency ranges might also result from extra-synaptic tonic inhibition active in cortical inhibitory neurons. Tonic inhibition increases the firing threshold and hence diminishes the output of inhibitory neurons to excitatory cortical neurons, which then allows increased excitation in the excitatory population and a power surge in the EEG. Conversely tonic inhibition in the thalamic relay cells does not induce power surge in EEG since augmented inhibition in the thalamic relay cells yields diminished excitation in cortical excitatory neurons, leading to a decrease in EEG power. This interpretation is corroborated by Figure 4 which demonstrates augmented and diminished non-linear gain in cortical excitatory neurons assuming tonic inhibition in inhibitory and thalamic relay population, respectively. This reflects enhanced and weakened response to the noisy thalamic external input, see previous theoretical studies (Hindriks and van Putten, 2012; Hutt, 2013) for a similar line of argument.
Figure 3 clearly reveals the emergence of δ−activity caused by extra-synaptic tonic inhibition which is affirmed by the existence of a minimum level of extra-synaptic inhibition shown in Figure 5. Conversely, α−activity appears to be much less sensitive to tonic inhibition since it is present for all tonic inhibition levels. One interpretation may be the generation of α−activity by the cortico-thalamic feedback as hypothesized theoretically (Robinson et al., 2004) while δ−activity results from the cortical interaction of excitatory and inhibitory neurons. The exact origins of propofol-induced α− and δ−activity are not known for certain. We find that the α−oscillations arise from thalamocortical resonances. These oscillations are commonly synchronous across widespread cortical regions and are not easily generated in isolated cortical tissue (Contreras et al., 1996; Destexhe et al., 1999). This affirms the original model of Lopes da Silva et al., (1974). However, our model results are equivalent to results of other models describing α−activity by purely cortical interactions. We are not aware of a methodology ruling out one or the other model and this is not the aim of the present work. Our work just reveals the additional possibility that the thalamus serves as a possible (indirect) source of α−activity. Similarly, the origin of δ−activity is not clear, but slow activity does increase at higher concentrations of propofol—which may be associated with decreasing α−waves as observed during desflurane general anesthesia (Mulholland et al., 2014). This is in keeping with δ−waves becoming more pronounced as the cortico-thalamic systems becomes increasingly hyperpolarized. However, there is a lot of variability between patients as regards the relative power of α− and δ−activity during general anesthesia; which would suggest that the true explanation is more complex, and requires recognition of other factors such as the one presented in this paper—the influence of the propofol on extra-synaptic inhibition.
Although anesthetic action on synaptic and extra-synaptic GABAergic receptors is different, both actions diminish neural activity and hence increase inhibition. Figure 5 elucidates that strong enough extra-synaptic or synaptic inhibition induce δ−activity. Hence, one may argue that the level of inhibition plays an important role while its origin, i.e., synaptic or extra-synaptic, plays a secondary role. This interpretation corroborates the idea of the balance of excitation and inhibition as the major mechanism in general anesthesia. This interpretation is in good accordance to previous experimental findings on the important role of the balance of excitation and inhibition in brain network under anesthesia (Okun and Lampl, 2008; Taub et al., 2013). Such global concepts as excitation-inhibition balance are attractive to describe complex processes in general anesthesia. For instance, anesthetics alter arousal in several pathways, such as the cholinergic pathway (Brown et al., 2010) and the orexinergic pathway which has been identified to activate a complex functional network controlling, inter alia, the emergence from unconsciousness (Kelz et al., 2008).
Our theoretical study considers the anesthetic propofol and its corresponding action at synaptic and extra-synaptic GABAergic receptors only, whereas it is known that propofol induces inhibition at various other receptors as well (Alkire et al., 2008; Nguyen et al., 2009) including minor effects on NMDA-receptors and voltage-gated potassium channels (Alkire et al., 2008). Propofol also potentiates glycine receptors which are are found all over the central nervous system and have a major role in regulating inhibition, e.g., in the brain stem (Lynch, 2004).
Similar to extra-synaptic inhibition resulting from ambient GABA concentrations, the presence of ambient concentrations of glycine close to NMDA-receptors entails tonic depolarization. This tonic excitation diminishes the firing threshold of neurons and hence may counteract inhibition. The present work considers tonic inhibition only and neglects tonic excitation effect. Although it would be important to study tonic excitation effects, this additional study would exceed the major aim of the manuscript, namely demonstrating the fundamental effect of tonic anesthetic action.
In addition, by virtue of the focus on extra-synaptic action, the model proposed neglects known anesthetic effects on different receptors and ion channels, although they have been shown experimentally, e.g., Grasshoff et al. (2006); Alkire et al. (2008) and references therein, and theoretically (Bojak et al., 2013) to affect EEG activity. Specifically, the latter work of Bojak et al. (2013) considers anesthetic effects on hyperpolarization-activated cyclic nucleotide-gated potassium channel 1 (HCN1) subunits which, effectively, increase the mean firing threshold in neural populations and strongly resembles the tonic inhibition induced by extra-synaptic GABA-receptors.
The model network topology includes a single module of a closed thalamo-cortical feedback loop (Granger and Hearn, 2007) comprising two thalamic nuclei and cortical excitatory and inhibitory neurons. This model represents a first approximation of brain networks since it neglects brain stem activity including the reticular activating system (RAS) (Magoun, 1952) which has significant modulating effects on attention, arousal and consciousness. Future work will include structures of the brain stem, propofol action on glycine receptors, and will take into account the RAS—since its neural structures involved exhibit strong extra-synaptic inhibition (Fiset et al., 1999; Franks, 2008; Vanini and Baghdoyan, 2013). The model also neglects the cholinergic pathway originating from the basal forebrain (Laalou et al., 2008) which is known to co-regulate the level of consciousness (Brown et al., 2010).
Essentially, our theoretical model assumes population coding implying rate-coding response activity of neuron populations subjected to external thalamic noise. The model does not consider specific single neuron dynamics found experimentally under anesthetic conditions. For instance, it has been hypothesized that, at certain levels of anesthetic concentration, thalamic neurons switch their activity from tonic firing to bursting and induce loss of consciousness (Alkire et al., 2000).
In spite of these limitations, our model reproduces qualitatively the action of propofol on EEG and reveals the possible impact of extra-synaptic GABAergic receptors on the EEG power. To our knowledge, the present work is the first to link extra-synaptic GABAergic action and experimental EEG. Future work will refine the model involving additional receptor action, e.g., tonic excitation caused by ambient glycine concentrations, and sub-cortical brain structures.
Meysam Hashemi has chosen the model and performed the model analysis, Axel Hutt has conceived the study and Jamie Sleigh has acquired and analyzed the experimental data. The authors have written the manuscript to equal parts and have approved the final version.
The research resulting to the presented work has received funding from the European Research Council under the European Unions Seventh Framework Programme (FP7/2007-2013) / ERC grant agreement no. 257253.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
The Supplementary Material for this article can be found online at: http://www.frontiersin.org/journal/10.3389/fnsys.2014.00232/abstract
Alkire, M., Haier, R. J., and Fallon, J. H. (2000). Toward a unified theory of narcosis: brain imaging evidence for a thalamocortical switch as the neurophysiologic basis of anesthetic-induced unconsciousness. Conscious. Cogn. 9, 370–386. doi: 10.1006/ccog.1999.0423
Bai, D., Zhu, G., Pennefather, P., Jackson, M. F., MacDonald, J. F., and Orser, B. (2001). Distinct functional and pharmacological properties of tonic and quantal inhibitory postsynaptic currents mediated by γ-aminobutyric acid a receptors in hippocampal neurons. Mol. Pharmacol. 59, 814–824. doi: 10.1124/mol.59.4.814
Belelli, D., Harrison, N. L., Maguire, J., MacDonald, R. L., Walker, M. C., and Cope, D. W. (2009). Extra-synaptic GABAA receptors: form, pharmacology, and function. J. Neurosci. 29, 12757–12763. doi: 10.1523/JNEUROSCI.3340-09.2009
Boly, M., Moran, R., Murphy, M., Boveroux, P., Bruno, M. A., Noirhomme, Q., et al. (2012). Connectivity changes underlying spectral EEG changes during propofol-induced loss of consciousness. J. Neurosci. 32, 7082–7090. doi: 10.1523/JNEUROSCI.3769-11.2012
Brickley, S. G., Cull-Candy, S. G., and Farrant, M. (1996). Development of a tonic form of synaptic inhibition in rat cerebellar granule cells resulting from persistent activation of GABAA receptors. J. Physiol. 497, 753–759.
Bright, D. P., Aller, M., and Brickley, S. (2007). Synaptic relase generates a tonic GABAA receptor-mediated conductance that modulates burst precision in thalamic relay neurons. J. Neurosci. 27, 2560–2569. doi: 10.1523/JNEUROSCI.5100-06.2007
Bright, D. P., Renzi, M., Bartram, J., McGee, T. P., MacKenzie, G., Hosie, A. M., et al. (2011). Profound desensitization by ambient GABA limits activation of delta-containing GABAA receptors during spillover. J. Neurosci. 31, 753–763. doi: 10.1523/JNEUROSCI.2996-10.2011
Cavalier, P., Hamann, M., Rossi, D., Mobbs, P., and Attwell, D. (2005). Tonic excitation and inhibition of neurons: ambient transmitter sources and computational consequences. Prog. Biophys. Mol. Biol. 87, 3–16. doi: 10.1016/j.pbiomolbio.2004.06.001
Ching, S., Cimenser, A., Purdon, P. L., Brown, E. N., and Kopell, N. J. (2010). Thalamocortical model for a propofol-induced-rhythm associated with loss of consciousness. Proc. Natl. Acad. Sci. U.S.A. 107, 22665–22670. doi: 10.1073/pnas.1017069108
Ching, S., Purdon, P. L., Vijayand, S., Kopell, N. J., and Brown, E. N. (2012). A neurophysiologicalmetabolic model for burst suppression. Proc. Natl. Acad. Sci. U.S.A. 109, 3095–3100. doi: 10.1073/pnas.1121461109
Cimenser, A., Purdon, P. L., Pierce, E. T., Walsh, J. L., Salazar-Gomez, A. F., Harrell, P. G., et al. (2011). Tracking brain states under general anesthesia by using global coherence analysis. Proc. Natl. Acad. Sci. U.S.A. 108, 8832–8837. doi: 10.1073/pnas.1017041108
Contreras, D., Destexhe, A., Sejnowski, T., and Steriade, M. (1996). Control of spatiotemporal coherence of a thalamic oscillation by corticothalamic feedback. Science 274, 771–774. doi: 10.1126/science.274.5288.771
Destexhe, A., Contreras, D., and Steriade, M. (1999). Mechanisms underlying the synchronizing action of corticothalamic feedback through inhibition of thalamic relay cells. J. Neurophysiol. 79, 999–1016.
Dickinson, R., Peterson, B., Banks, P., Simillis, C., Martin, J., Valenzuela, C., et al. (2007). Competitive inhibition at the glycine site of the n-methyl-d-aspartate receptor by the anesthetics xenon and isoflurane: evidence from molecular modeling and electrophysiology. Anesthesiology 107, 756–767. doi: 10.1097/01.anes.0000287061.77674.71
Fiset, P., Paus, T., Daloze, T., Plourde, G., Meuret, P., Bonhomme, V., et al. (1999). Brain mechanisms of propofol-induced loss of consciousness in humans: a positron emission tomographic study. J. Neurosci. 19, 5506–5513.
Fleming, T., Scott, V., Naskar, K., Joe, N., Brown, C., and Stern, J. (2011). State-dependent changes in astrocyte regulation of extrasynaptic nmda receptor signalling in neurosecretory neurons. J. Physiol. 589, 3929–3941. doi: 10.1113/jphysiol.2011.207340
Friedman, E. B., Sun, Y., Moore, J. T., Hung, H.-T., Meng, Q. C., Perera, P., et al. (2010). A conserved behavioral state barrier impedes transitions between anesthetic-induced unconsciousness and wakefulness: evidence for neural inertia. PLoS ONE 5:e11903. doi: 10.1371/journal.pone.0011903
Gugino, L. D., Chabot, R. J., Prichep, L. S., John, E. R., Formanek, V., and Aglio, L. S. (2001). Quantitative EEG changes associated with loss and return of conscious- ness in healthy adult volunteers anaesthetized with propofol or sevoflurane. Br. J. Anaesth. 87, 421–428. doi: 10.1093/bja/87.3.421
Hamann, M., Rossi, D., and Attwell, D. (2002). Tonic and spillover inhibition of grnule cells control information flow through cerebellar cortex. Neuron 33, 625–633. doi: 10.1016/S0896-6273(02)00593-7
Houston, C., McGee, T., MacKenzie, G., Troyano-Cuturi, K., Mateos Rodriguez, P., Kutsarova, E., et al. (2012). Are extrasynaptic gabaa receptors important targets for sedative/hypnotic drugs? J. Neurosci. 32, 3887–3897. doi: 10.1523/JNEUROSCI.5406-11.2012
Hutt, A. (2012). The population firing rate in the presence of GABAergic tonic inhibition in single neurons and application to general anaesthesia. Cogn. Neurodyn. 6, 227–237. doi: 10.1007/s11571-011-9182-9
Hutt, A. (2013). The anaesthetic propofol shifts the frequency of maximum spectral power in EEG during general anaesthesia: analytical insights from a linear model. Front. Comput. Neurosci. 7:2. doi: 10.3389/fncom.2013.00002
Hutt, A., and Buhry, L. (2014). Study of GABAergic extra-synaptic tonic inhibition in single neurons and neural populations by traversing neural scales: application to propofol-induced anaesthesia. J. Comput. Neurosci. 37, 417–437. doi: 10.1007/s10827-014-0512-x
Johnson, B., Sleigh, J., Kirk, I., and Williams, M. (2003). High-density EEG mapping during general anaesthesia with xenon and propofol: a pilot study. Anaesth. Intensive Care. 31, 155–163. Available online at: http://www.aaic.net.au/Document/D=2002375
Kelz, M., Sun, Y., Chen, J., Cheng Meng, Q., Moore, J., Veasey, S., et al. (2008). An essential role for orexins in emergence from general anesthesia. Proc. Natl. Acad. Sci. U.S.A. 105, 1309–1314. doi: 10.1073/pnas.0707146105
Kitamura, A., Marszalec, W., Yeh, J., and Narahashi, T. (2002). Effects of halothane and propofol on excitatory and inhibitory synaptic transmission in rat cortical neurons. J. Pharmacol. 304, 162–171. doi: 10.1124/jpet.102.043273
Kratzer, S., Mattusch, C., Kochs, E., Eder, M., Haseneder, R., and Rammes, G. (2012). Xenon attenuates hippocampal long-term potentiation by diminishing synaptic and extrasynaptic n-methyl-d- aspartate receptor currents. Anesthesiology 116, 673–682. doi: 10.1097/ALN.0b013e3182475d66
Kretschmannova, K., Hines, R. M., Revilla-Sanchez, R., Terunuma, M., Tretter, V., Jurd, R., et al. (2013). Enhanced tonic inhibition influences the hypnotic and amnestic actions of the intravenous anesthetics etomidate and propofol. J. Neurosci. 33, 7264–7273. doi: 10.1523/JNEUROSCI.5475-12.2013
Laalou, F., de Vasconcelos, A., Oberling, P., Jeltsch, H., Cassel, J., and Pain, L. (2008). Involvement of the basal cholinergic forebrain in the mediation of general (propofol) anesthesia. Anesthesiology 108, 888–896. doi: 10.1097/ALN.0b013e31816d919b
Le Meur, K., Galante, M., Angulo, M.-C., and Audinat, E. (2007). Tonic activation of nmda receptors by ambient glutamate of non-synaptic origin in the rat hippocampus. J. Physiol. 580, 373–383. doi: 10.1113/jphysiol.2006.123570
Lewis, L., Weiner, V., Mukamel, E., Donoghue, J., Eskandar, E., Madsen, J., et al. (2012). Rapid fragmentation of neuronal networks at the onset of propofol-induced unconsciousness. Proc. Natl. Acad. Sci. U.S.A. 109, E3377–E3386. doi: 10.1073/pnas.1210907109
Liley, D., and Bojak, I. (2005). Understanding the transition to seizure by modeling the epileptiform activity of general anaesthetic agents. J. Clin. Neurophysiol. 22, 300–313. Available online at: https://wwws.lww.com/clinicalneurophys/Abstract/2005/10000/Understanding_the_Transition_to_Seizure_by.3.aspx
Martin, D., Plagenhoef, M., Abraham, J., Dennison, R., and Aronstam, R. (1995). Volatile anesthetics and glutamate activation of n-methyl-d-aspartate receptors. Biochem. Pharmacol. 49, 809–817. doi: 10.1016/0006-2952(94)00519-R
McCarthy, M. M., Brown, E. N., and Kopell, N. (2008). Potential network mechanisms mediating electroencephalographic beta rhythm changes during propofol-induced paradoxical excitation. J. Neurosci. 28, 13488–13504. doi: 10.1523/JNEUROSCI.3536-08.2008
McDougall, S. J., Bailey, T. W., Mendelowitz, D., and Andresen, M. C. (2008). Propofol enhances both tonic and phasic inhibitory currents in second-order neurons of the solitary tract nucleus (nts). Neuropharmacology 54, 552–563. doi: 10.1016/j.neuropharm.2007.11.001
Mothet, J., Parent, A., Wolosker, H., Brady, Jr., R, Linden, D., Ferris, C., et al. (2000). d-serine is an endogenous ligand for the glycine site of the n-methyl-d-aspartate receptor. Proc. Natl. Acad. Sci. U.S.A. 97, 4926–4931. doi: 10.1073/pnas.97.9.4926
Mulholland, C., Somogyi, A., Barratt, D., Coller, J., Hutchinson, M., Jacobson, G., et al. (2014). Association of innate immune single-nucleotide polymorphisms with the electroencephalogram during desflurane general anaesthesia. J. Mol. Neurosci. 52, 497–506. doi: 10.1007/s12031-013-0201-7
Murphy, M., Bruno, M.-A., Riedner, B. A., Boveroux, P., Noirhomme, Q., Landsness, E. C., et al. (2011). Propofol anesthesia and sleep: a high-density EEG study. Sleep 34, 283–291. Available online at: http://www.journalsleep.org/ViewAbstract.aspx?pid=28064
Nguyen, H., Li, K., da Graca, R., Delphin, E., Xiong, M., and Ye, J. (2009). Behavior and cellular evidence for propofol-induced hypnosis involving brain glycine receptors. Anesthesiology 110, 326–332. doi: 10.1097/ALN.0b013e3181942b5b
Orser, B. (2006). Extrasynaptic GABAA receptors are critical targets for sedative-hypnotic drugs. J. Clin. Sleep Med. 2, S12–S18. Available online at: http://www.aasmnet.org/jcsm/ViewAbstract.aspx?pid=26526
Panatier, A., Theodosis, D., Mothet, J.-P., Touquet, B., Pollegioni, L., Poulain, D., et al. (2006). Glia-derived d-serine controls nmda receptor activity and synaptic memory. Cell 125, 775–784. doi: 10.1016/j.cell.2006.02.051
Papouin, T., Ladépêche, L., Ruel, J., Sacchi, S., Labasque, M., Hanini, M., et al. (2012). Synaptic and extrasynaptic nmda receptors are gated by different endogeneous coagonists. Cell 150, 633–646. doi: 10.1016/j.cell.2012.06.029
Purdon, P. L., Pierce, E. T., Mukamel, E. A., Prerau, M. J., Walsh, J. L., Wong, K. F., et al. (2012). Electroencephalogram signatures of loss and recovery of consciousness from propofol. Proc. Natl. Acad. Sci. U.S.A 110, E1142–E1150. doi: 10.1073/pnas.1221180110
Robinson, P., Rennie, C. J., Rowe, D. L., and O'Connor, S. C. (2004). Estimation of multiscale neurophysiologic parameters by electroencephalographic means. Hum. Brain Mapp. 23, 53–72. doi: 10.1002/hbm.20032
Robinson, P., Loxley, P., O'Connor, S., and Rennie, C. (2001). Modal analysis of corticothalamic dynamics, electroencephalographic spectra and evoked potentials. Phys. Rev. E Stat. Nonlin. Soft. Matter Phys. 63:041909. doi: 10.1103/PhysRevE.63.041909
Robinson, P., Rennie, C., and Rowe, D. (2002). Dynamics of large-scale brain activity in normal arousal states and epileptic seizures. Phys. Rev. E Stat. Nonlin. Soft. Matter Phys. 65:041924. doi: 10.1103/PhysRevE.65.041924
Sellers, K. K., Bennett, D. V., Hutt, A., and Frohlich, F. (2013). Anesthesia differentially modulates spontaneous network dynamics by cortical area and layer. J. Neurophysiol. 110, 2739–2751. doi: 10.1152/jn.00404.2013
Semyanov, A., Walker, M. C., Kullmann, D. M., and Silver, R. A. (2004). Tonically active GABAA receptors: modulating gain and maintaining the tone. Trends Neurosci. 27, 262–269. doi: 10.1016/j.tins.2004.03.005
Steyn-Ross, M., Steyn-Ross, D., and Sleigh, J. (2004). Modelling general anaesthesia as a first-order phase transition in the cortex. Prog. Biophys. Mol. Biol. 85, 369–385. doi: 10.1016/j.pbiomolbio.2004.02.001
Steyn-Ross, M., Steyn-Ross, D., and Sleigh, J. (2013). Interacting turing-hopf instabilities drive symmetry-breaking transitions in a mean-field model of the cortex: a mechanism for the slow oscillation. Phys. Rev. X 3:021005. doi: 10.1103/PhysRevX.3.021005
Steyn-Ross, M., Steyn-Ross, D., Sleigh, J. W., and Liley, D. T. J. (1999). Theoretical electroencephalogram stationary spectrum for a white-noise-driven cortex: evidence for a general anesthetic-induced phase transition. Phys. Rev. E Stat. Phys. Plasmas Fluids Relat. Interdiscip. Topics 60, 7299–7311. doi: 10.1103/PhysRevE.60.7299
Taub, A., Katz, Y., and Lampl, I. (2013). Cortical balance of excitation and inhibition is regulated by the rate of synaptic activity. J. Neurosci. 33, 14359–14368. doi: 10.1523/JNEUROSCI.1748-13.2013
Victor, J., Drover, J., Conte, M., and Schiff, N. (2011). Mean-field modeling of thalamocortical dynamics and a model-driven approach to EEG analysis. Proc. Natl. Acad. Sci. U.S.A. 118, 15631–15638. doi: 10.1073/pnas.1012168108
Vizuete, J., Pillay, S., Ropella, K., and Hudetz, A. (2014). Graded defragmentation of cortical neuronal firing during recovery of consciousness in rats. Neuroscience 275, 340–351. doi: 10.1016/j.neuroscience.2014.06.018
Wei, W., Zhang, N., Peng, Z., Houser, C. R., and Mody, I. (2003). Perisynaptic localization of subunit-containing GABAA receptors and their activation by GABA spillover in the mouse dentate gyrus. J. Neurosci. 23, 10650–10661. Available online at: http://www.jneurosci.org/content/23/33/10650.short
Wilson, M., Sleigh, J., Steyn-Ross, A., and Steyn-Ross, M. (2006). General anesthetic-induced seizures can be explained by a mean-field model of cortical dynamics. Anesthesiology 104, 588–593. doi: 10.1097/00000542-200603000-00026
Wolosker, H., Blackshaw, S., and Snyder, S. (1999). Serine racemase: a glial enzyme synthesizing d-serine to regulate glutamate-n-methyl-d-aspartate neurotransmission. Proc. Natl. Acad. Sci. U.S.A. 96, 13409–13414. doi: 10.1073/pnas.96.23.13409
Ye, Z., McGee, T., Houston, C., and Brickley, S. (2013). The contribution of δ subunit-containing gabaa receptors to phasic and tonic conductance changes in cerebellum, thalamus and neocortex. Front. Neural Circuit. 7, 1–8. doi: 10.3389/fncir.2013.00203
Yeung, J. Y. T., Canning, K. J., Zhu, G., Pennefather, P., MacDonald, J. F., and Orser, B. A. (2003). Tonically activated GABAA receptors in hippocampal neurons are high-affinity, low-conductance sensors for extracellular gaba. Mol. Pharmacol. 63, 2–8. doi: 10.1124/mol.63.1.2
Keywords: GABA receptor, neural mass, propofol, power spectrum, general anesthesia
Citation: Hashemi M, Hutt A and Sleigh J (2014) Anesthetic action on extra-synaptic receptors: effects in neural population models of EEG activity. Front. Syst. Neurosci. 8:232. doi: 10.3389/fnsys.2014.00232
Received: 12 September 2014; Accepted: 19 November 2014;
Published online: 10 December 2014.
Edited by:Anthony G. Hudetz, Medical College of Wisconsin, USA
Reviewed by:Ingo Bojak, University of Reading, UK
Kingsley Paul Storer, Weill Cornell Medical College, USA
Copyright © 2014 Hashemi, Hutt and Sleigh. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Axel Hutt, INRIA CR Nancy - Grand, 615 Rue du Jardin Botanique, Villers-les-Nancy, France e-mail: email@example.com