Clusters of cooperative ion channels enable a membrane-potential-based mechanism for short-term memory

  1. Paul Pfeiffer  Is a corresponding author
  2. Alexei V Egorov
  3. Franziska Lorenz
  4. Jan-Hendrik Schleimer
  5. Andreas Draguhn
  6. Susanne Schreiber  Is a corresponding author
  1. Humboldt-Universität zu Berlin, Germany
  2. Heidelberg University, Germany

Abstract

Across biological systems, cooperativity between proteins enables fast actions, supra-linear responses, and long-lasting molecular switches. In the nervous system, however, the function of cooperative interactions between voltage-dependent ionic channels remains largely unknown. Based on mathematical modeling, we here demonstrate that clusters of strongly cooperative ion channels can plausibly form bistable conductances. Consequently, clusters are permanently switched on by neuronal spiking, switched off by strong hyperpolarization, and remain in their state for seconds after stimulation. The resulting short-term memory of the membrane potential allows to generate persistent firing when clusters of cooperative channels are present together with non-cooperative spike-generating conductances. Dynamic clamp experiments in rodent cortical neurons confirm that channel cooperativity can robustly induce graded persistent activity – a single-cell based, multistable mnemonic firing mode experimentally observed in several brain regions. We therefore propose that ion channel cooperativity constitutes an efficient cell-intrinsic implementation for short-term memories at the voltage level.

Introduction

Cooperative molecular interactions are ubiquitous in biology and guide cellular processes from sensing to memory formation (Bray et al., 1998; Burrill and Silver, 2010). They are found not only in the simplest organisms like bacteria, but also in higher organisms including mammals (Stefan and Le Novère, 2013). Evidence increases that also ion channels of excitable membranes in the heart and the nervous system, including the mammalian brain, can exhibit cooperative properties (Choi, 2014; Molina et al., 2006; Grage et al., 2011; Kim et al., 2014; Gianoli et al., 2017; Navedo et al., 2010; Dixon et al., 2015; Moreno et al., 2016; Clatot et al., 2017; Vivas et al., 2017). Nevertheless, the majority of neuron models relies on the independent gating assumption: channels communicate indirectly via the common membrane potential, but do not directly influence each other. Thus, it is an open question how cooperative channels affect the electrical dynamics - and therefore the computations - of a neuron.

Experimentally, ion channel cooperativity has been studied in various channel types with key roles in the nervous system such as potassium (Molina et al., 2006; Kim et al., 2014), sodium (Iwasa et al., 1986; Undrovinas et al., 1992; Clatot et al., 2017), HCN (Dekker and Yellen, 2006) and calcium channels (Dixon et al., 2015; Moreno et al., 2016). However, most studies have focused on the demonstration of cooperative gating in small ensembles of channels and avoided the complex dynamics of the neuron as a whole. An exception is the study of cooperative calcium channels in rodent hippocampal neurons, where optical control of channel coupling has revealed that the spontaneous firing rate rises when channels cooperate (Moreno et al., 2016). With such control of channel interactions, future experiments have the means to test hypothesis on the function of cooperativity.

So far, a few computational studies predict effects of cooperativity. Along these lines, cooperative sodium channels have been suggested to underlie the rapid initiation and low variability of spiking onset in cortical neurons (Naundorf et al., 2006; Huang et al., 2012; for an alternative mechanism see Yu et al., 2008) and mild cooperative interactions in potassium and calcium channels have been shown to modify the steepness of a channel’s activation curve and modulate neural excitability (Zarubin et al., 2012). Furthermore, cardiac alternans, a pathological condition of unstable heart contractions, has been linked to strong degrees of coupling among cooperative calcium channels (Sato et al., 2018). These studies demonstrate that cooperative channels can significantly alter cellular firing properties with a range of effects from advantageous to pathological depending on the interplay with other currents in the cell.

Here, we show in simulations and mathematical analysis that small clusters of cooperative channels with simple and generic activation dynamics can induce a multistability of the membrane potential and demonstrate that this multistability enables a form of cellular memory. Our central observation is that the mutually enhancing nature of cooperative gating favors joint opening and closing of channels, and, more importantly, results in a hysteresis of their gating behavior. We demonstrate that such cooperative channels - when arranged in clusters and located in membranes with ‘normal’ independent conductances that mediate spiking - can induce graded persistent neural activity: spiking that persists after transient suprathreshold depolarization and represents successive inputs with increasing persistent firing rates. Taken together, we propose that cooperativity of a few ion channels is an elegant way to implement a cell-intrinsic memory directly reflected in a cell’s firing. In principle, ion channel cooperativity could thus efficiently complement network-based mechanisms of persistent activation and thereby contribute to decision making (accumulation of evidence) and working memory. Overall, the emergence of memory in a cluster of coupled, yet in isolation memoryless channels suggests a more general design principle: cooperativity serves to dynamically build functionally rich macrochannels from simpler channels.

In the following, we first dissect the mechanism and compare the gating of independent versus cooperative ion channel clusters. We then add clusters of cooperative channels to a simple neuron model with non-cooperative, spike-generating sodium and potassium channels, proving the ability of ion channel cooperativity to mediate graded persistent activity. Finally, we use the dynamic clamp technique to experimentally endow perirhinal cortex neurons with ‘virtual’ cooperative ion channels and show that stable persistent activity can be robustly induced, rendering ion channel cooperativity an efficient and plausible mechanism for cellular voltage memory.

Results

Although the mechanisms of cooperative interactions among ion channels are still not fully understood on the molecular level, it can be assumed that ion channels need to be in spatial proximity to directly interact and gate in a cooperative manner (Gutkin and Ermentrout, 2006). Matching this assumption, cooperativity is often found in channels that form small clusters (Molina et al., 2006; Choi, 2014; Navedo et al., 2010; Moreno et al., 2016). For our study, we therefore assume that the cooperative channels are distributed in multiple clusters and limit cooperative interactions to channels within the same cluster. In this regard, we deviate from previous theoretical studies, where all-to-all interactions between channels in a membrane have been investigated (Naundorf et al., 2006; Zarubin et al., 2012; Huang et al., 2012). To understand the underlying principles, we begin with a comparison of gating properties in a cluster of cooperative versus non-cooperative ion channels.

Gating of a cluster of cooperative channels

A cluster is assumed to comprise S voltage-gated channels of the same type, whose individual activation curve is a monotonically increasing, yet saturating function of voltage. To ensure the generality of the analysis, we make no further assumptions on other channel properties like their ionic nature and therefore simply refer to them as ‘channels’. Cooperative interactions are chosen to be enhancing, that is the opening of a channel in the cluster increases the opening probabilities of all other S-1 channels. This increase in opening probability is implemented phenomenologically by shifting the channels’ activation curves to lower voltages for each new opening of a neighbor. The size of the shift j induced by one new opening is a measure of the coupling between channels in the cluster (Naundorf et al., 2006; Zarubin et al., 2012; Huang et al., 2012). An important determinant for the cluster dynamics is the maximal shift J=(S-1)j, which a channel experiences, when all neighbors are in the open state (Figure 1A). Variation of this interaction strength from independent (J = 0 mV) to strongly cooperative (J = 70 mV) yields the activation curve of a cluster of S channels. Clearly, an isolated channel and channels in an independent ensemble have the same activation curve. Mild cooperative interactions, however, increase the channel activity and make the activation curve steeper. Beyond a critical coupling, the activation curve ‘bends over’, which characterizes the regime of strong cooperativity (Figure 1B).

Gating of a cluster of cooperative channels.

(A) Model of cooperative gating: the activation of a channel in a cluster depends on the state of the surrounding channels. Opening of a neighboring channel leads to a shift of the activation curve by the coupling strength j. When all neighbors are closed, the activation coincides with the one of an isolated channel, whereas when all neighbors are open, the activation undergoes the maximal shift J. (B) Cluster activation: cooperative interactions inside a cluster increase the channel activity, so that the activation curve becomes steeper (weak coop.). Above a critical coupling strength, the activation starts to ‘bend over’ (strong coop.); whether a channel is open or closed is determined mostly by the state of its surrounding channels. Activation curve based on a self-consistency relation, for details see Materials and methods. (C) Top: Simulation of a voltage clamp experiment with an isolated two-state channel shows fast, random switching between the conducting (black) and non-conducting (white) state. Operated at the half-activation potential V0.5, the channel spends about half of the time in the open state. Middle: In a cluster of independent channels, their asynchronous gating results in a fluctuation of the cluster conductance around half of the total conductance. Bottom: Cooperative channels have a strong preference for the state of the surrounding channels. In a cluster, they open and close in synchrony acting like a macrochannel. Its slow switching frequency demonstrates the stability of the open and closed cluster state. As the cooperative channels are more active compared to the independent ones, they are clamped at more hyperpolarized voltages: −1 mV (independent) and 36 mV (cooperative). The details of the jump process simulation are given in the section Materials and methods and parameters are summarized in Table 1.

Table 1
Cluster and cooperative channel parameters.

N: number of clusters, S: cluster size, J: overall coupling, j: channel-channel coupling, V0.5: half activation, k: width of the channel activation, τ: time constant of the channel with maximum at attained at Vm and width σ, gcoop: single channel conductance and Ecoop: reversal potential.

FigureNSJjV0.5kτVmσgcoopEcoop
mVmVmVmVmsmVmVpSmV
Figure 1 (independent)1600−1150.5−130--
Figure 1 (weak coop.)1622.54.5−1150.5−130--
Figure 1 (strong coop.)167014−1150.5−130--
Figure 2 (independent)1600−1150.5−130--
Figure 2 (cooperative)1610020−1150.5−130--
Figure 31510025−1150.5−130--
Figure 410088011.4−3010120−30202.5100
Figure 510088011.4−3010120−30202.5100
Figure 6c908101.514.5−1015100−10301100
Figure 3—figure supplement 11510025−1150.5−130--
Figure 4—figure supplement 1A10088011.4−3010120−30202.5100
Figure 4—figure supplement 1B10088011.4−3010120−302010−100
Figure 4—figure supplement 210088011.4−3010120−30202.5100
Figure 5—figure supplement 1A,B10088011.4−301010–240−30202.5100
Figure 5—figure supplement 1C,D10088011.4−3010120−30200.5–10100
Figure 6—figure supplement 1 (independent channels, red)90800−1015100−10301100
Figure 6—figure supplement 1 (cooperative channels, black)908101.514.5−1015100−10301100

Simulating a simple voltage clamp experiment in a membrane patch containing a cluster of six channels demonstrates fundamental differences between the gating of independent and cooperative clusters: The conductance dynamics of a cluster of independent channels is identical to the sum of its individual, independent components. Therefore, the single channel opening probability dictates the most probable cluster state (defined by the channel activation curve and the voltage) and the cluster conductance fluctuates around this state following a binomial distribution (Figure 1C, top). At the half-activation voltage, the cluster average is at half of its total conductance. Fluctuations in the total cluster conductance are relatively fast (due to channel gating time constants of ≈ 1 ms).

The picture changes when ion channel cooperativity is introduced. Assuming strong interactions, all channels in a cluster gate open and close quasi-simultaneously (at a voltage level within the regime where the activation curve ‘bends over’). Intuitively, this coordination within a cluster can be easily understood: the opening of ion channels strongly enhances the probability of their neighbors to open. Similarly, the closing of channels reduces the probability of their neighbors to remain open, so that gating dynamics largely synchronize. The distribution of total cluster conductance now deviates from the unimodal, binomial distribution of clusters with independent channels: it becomes bimodal reflecting the alternation between a fully closed and a fully open cluster state (Figure 1C, bottom). Note that the bistability is accompanied by a novel slow time scale; the cluster switches between the all-open or all-closed states at a drastically reduced rate (6 Hz) compared to its constituting channels when independent (1000 Hz). Effectively, the bistability and the slow switching let the cluster of cooperative channels appear as one slow macrochannel (Blunck et al., 1998).

Hysteresis and bistable gating

The new, slow timescale also changes the response of the cluster to a pulsatile stimulation protocol. During a sufficiently strong depolarizing pulse, all channels in both the cooperative and independent cluster open. Once the pulse is over and voltage returns to baseline, the cluster of independent channels follows and channels swiftly return to the closed state. In the case of the cooperative cluster, however, channels in the cluster remain open, despite the return of the clamped voltage value to its original value before the pulse onset. The strong channel interactions introduce what we would like to term ‘stickiness’: the cluster ‘sticks’ to its new open state (even when voltage levels are back to ‘normal’) and can only be released from this state by a sufficiently strong hyperpolarization (Figure 2A). In other words, cooperative gating induces hysteresis and hence allows channels in a cluster, which are fast switching in isolation, to represent (i.e. ‘remember’) voltage values that occurred multiple hundred milliseconds ago.

Hysteresis and bistability of a cluster of cooperative channels.

(A) Probing an initially closed cluster (time point 0), both independent (blue) and cooperative (red) channels open in response to a depolarizing voltage pulse (1). However, returning to the base line voltage (2), independent channels rapidly close again, whereas a cluster of cooperative channels ‘sticks’ to the open configuration. Only a strong hyperpolarizing pulse (3) can close the cooperative channels and restore the initial condition (4). Tracking membrane potential versus conductance during the pulse protocol (right panels), the hysteresis in gating for the cooperative channels becomes apparent; whether the cluster is open or closed at baseline voltage depends on the previous voltage pulse. The transparent lines correspond to 20 repetitions of the pulse protocol and the solid lines are the respective means. (B) Strong coupling among channels changes the activation in a cluster (red) with respect to an isolated channel (black) and induces a voltage regime of bistability (red shaded). (C) Clusters only become bistable beyond a critical coupling strength (black dashed) and for increasing coupling, the bistable range extends to lower voltages. Cluster and channel parameters are summarized in Table 1.

The hysteresis in gating arises from the bistability of the cluster. Before the depolarizing pulse, the channels in the cluster are closed, which is a stable cluster state, because the single channel is not activated at this membrane potential and there is no facilitation between closed channels. During the pulse, channels open, and this open cluster configuration is stable at the baseline voltage, because the strong coupling allows the channels to secure each other in the open state. Only when the voltage is lowered further, it overcomes the mutual facilitation and the channels close again.

The bistability reflects those effective cluster activation curves that ‘bend over’ (Figure 2B). Consequently, cluster size, channel coupling and the shape of the isolated channel’s activation curve determine the voltage range of bistability. Both an increase in channel coupling j and cluster size S can lead to a stronger overall cooperativity strength J, which is the product of the two. Bistability arises from a critical overall cooperativity strength Jcrit, which depends on the shape of the single channel activation. With increasing J, the bistable range broadens and extends to more hyperpolarized voltages (Figure 2C). Correspondingly, for a cluster with more channels and/or strong coupling, which both result in a larger J, the voltage pulse to close the cluster has to be stronger in order to leave the bistable range.

Due to their self-excitatory effect, channels in a cooperative cluster can be compared to neurons in a recurrently connected network: if the external input is low, the neurons are silent, but during a short stimulation they switch to a very active state, which they sustain after the stimulation by recurrent self-excitation. Such networks have been proposed to underlie the short-term storage of stimuli in working memory (Durstewitz et al., 2000). This poses the question, whether in a similar way, a cluster of cooperative channels can act as a memory unit for a single neuron. In order to address this question, we first need to better understand the emerging slow timescale of switches between cluster states in the cooperative regime.

Prolonged lifetimes of the open and closed cluster state

As we have seen, hysteresis allows clusters to act as a voltage memory. After high voltages, they remain open and after low voltages, they remain closed. Eventually, however, a cluster spontaneously switches its state - in a quasi-synchronous fashion all channels open or close (Figure 3A). These switches originate in channel noise, stochastic gating of individual channels. If multiple channels coincide in their spontaneous gating, their neighbors quickly follow. Both opening and closing of channels spread by cooperative coupling; in a spiral of facilitation build-up or, reversely, in a spiral of facilitation loss (Figure 3—figure supplement 1). These gating avalanches, however, cannot be triggered by single or few channels, so that the life times of a cluster can be orders of magnitude longer than those of the single channel.

Figure 3 shows the lifetimes of the open and closed states - exemplary in a cluster with five channels - as a function of voltage. These lifetimes are defined as the time τOC it takes, once all channels are open until all channels are closed again, and vice versa for τCO. In the center of the bistable range, both open and closed cluster state are very robust against channel fluctuations with mean residence times in the range of seconds (Figure 3A). However, at the verges of the bistable range, channel noise corrupts the stability. At low voltages, an open cluster switches back to the closed state within ≈ 10 ms, whereas at high voltages, within in the same time, a closed cluster spontaneously opens. Thus, for a persistent representation of recent voltage history, the middle voltage range offers the longest memory lifetimes: here both the probability to lose an open cluster state and to lose a closed cluster state are low - a prerequisite for encoding of previous voltages in a persisting cluster state (Figure 3A). The lifetimes in this voltage range define the memory timescale τmax of the cluster, that is the time during which a cluster of cooperative ion channels can act as a reliable memory device. For completeness, Figure 3B shows the numerically derived center of the bistable range Vmax and the memory timescale τmax for clusters of different sizes and with different channel-channel coupling j. We find that both, larger cluster size and stronger channel coupling, favor bistability. Although quantitative data on cooperative interactions is scarce, we note that for a cluster of eight channels coupled with strength j = 17 mV, the mid-voltage of the bistable range decreases to −60 mV and the average lifetime of the open and closed state exceeds hundreds of seconds, rending this paradigm suitable for persistent activation in the context of short-term memory.

In general, our analysis shows that, as expected, larger and more strongly coupled clusters have a more hyperpolarized operation voltage and longer life times. Even for clusters with the same maximal shift J (that is hyperboles in Figure 3B), we observe a pronounced increase in lifetime with cluster size. In other words, large clusters of weakly coupled channels tend to be more stable than small clusters of strongly coupled ones. This observation is explained by a reduction of effective channel noise when more channels are added to a cluster.

Figure 3 with 1 supplement see all
Prolonged life times of cluster states in the bistable regime.

(A) Top: Voltage-clamped in the bistable regime, a cluster of five cooperative channels can stay in the open or closed state for multiple seconds until it spontaneously switches. The noise source are spontaneous single channel gating events, visible as the fast fluctuations around the stable states. Bottom: Entering the bistable regime, the residence times in the open and closed state increase over multiple orders of magnitude. The point of maximal stability (Vmax,τmax) is located at the center; here, a cluster is expected to stay open or closed for seconds, exceeding by far the single channel time constant τ(V) (black line). Mean life times are derived from a first passage analysis, see Materials and methods. (B) Stability of the open and closed cluster configuration depend on the cluster size and the intra channel coupling j. Larger clusters with stronger intra-channel coupling can be stable for multiple seconds and their bistable range moves to more hyperpolarized voltages. Note that small clusters with weak coupling are not bistable (white region). Cluster and channel parameters are summarized in Table 1.

Clusters as cellular memory units can mediate persistent spiking

The analysis so far showed that clusters of cooperative ion channels can implement a memory of recent voltage levels in their opening state. Next, we demonstrate how this property - when exhibited by small clusters of channels not even directly involved in the generation of spikes - can lead to persistent neuronal activity that does not require further network input. Interestingly, we find that clusters of cooperative channels can even solve the computationally harder problem required for a graded form of persistent activity. In this graded form, the neuron can signal previous input strength and input duration in the frequency of long-term stable firing. Hence, we propose that ion channel cooperativity offers a generic mechanism for cellular memory.

Neural signaling events lead to controllable, persistent cluster switches

To this end, we extend a Hodgkin-Huxley type neuron model including (independent) sodium and potassium channels by a set of 100 small clusters of cooperative channels with the generic activation kinetics described above (for details on the neuron model, see Materials and methods). On the ionic nature of the latter, we assume for the moment that they conduct a depolarizing current like for example calcium channels, which are a potential candidate (Moreno et al., 2016). For the small number of cooperative channels, we simulate their stochastic gating, whereas for the other ionic currents, we use a deterministic simulation.

Stimulating the model cell with transient input pulses elicits spiking at a rate that increases with the stimulus amplitude (Figure 4). Correspondingly, the stimulation amplitude determines the response of the cooperative channels: At low input amplitudes (and consequently low ensuing firing rates), few channels in the clusters open transiently, but close again rapidly (Figure 4A, top trace). When the pulse is over, the cell falls silent and does not generate spikes. In contrast, if the firing rate during the pulse is large enough, in a subset of clusters all channels open. Moreover, in agreement with our above observations, these clusters now remain open when the input pulse is over (Figure 4A, middle and bottom traces). This additional depolarizing conductance suffices to keep the cell in a firing modus and the model neuron exhibits firing activity that persists beyond the presentation of the stimulus pulse.

Figure 4 with 2 supplements see all
Clusters of cooperative channels mediate persistent activity at different levels in a model neuron.

(A) Current stimulation with different amplitudes to test the response of clusters to different firing rates. Left: Voltage traces demonstrating stimulated firing at increasing frequencies. Whereas the neuron returns to rest after a low amplitude stimulation, after a stronger drive the neuron continues to spike at a stable frequency fPA, which increases with the frequency of the stimulated firing fdrive. Right: Persistent activity is mediated by the conductance of the clusters, which builds up during high-frequency spiking and remains stable during low-frequency spiking and at rest. (B) A strong hyperpolarizing step closes the clusters and stops persistent firing. (C) Left: The cooperative clusters track the input strength; more clusters open when the cell fires at higher frequencies. Firing below 20 Hz, however, does not open clusters. The number of open clusters can differ across trials (black dots) because of stochastic channel gating. Right: Persistent activity increases with the number of open clusters, thereby allowing the neuron to represent the input strength after the stimulus has ended, for example 3 Hz after 29 Hz firing and 9 Hz after 51 Hz firing. Big dots denote the mean for each stimulation strength, small dots the individual trials. The simulation procedure is described in Materials and methods and parameters are summarized in Table 1.

Furthermore, the larger the amplitude of the stimulating pulse, the larger the persistent firing rate following the pulse (owing to the persistent opening of more clusters). As every cluster of cooperative channels has a stable open and closed configuration, the population conductance across all clusters allows for a quasi-continuous range of persistent currents. When all cluster are open, the maximal current is reached and the persistent firing rate saturates (≈ 15 Hz). According to the gating properties of cooperative channel clusters, sufficiently hyperpolarizing pulses should be able to close clusters and hence ‘reset’ the memory of the previous pulse by annihilation of the persistent activity. Indeed, Figure 4B shows that a hyperpolarizing pulse cancels persistent activity.

In order to ensure maximal stability when there are no inputs to the neuron, the center of the cluster bistability (−70 mV) coincides with the resting potential (−67 mV). The limited extent of the bistable range from around −90 mV to −50 mV allows the clusters to react to spikes or strong inhibition. For the occurrence of graded levels of persistent activity, we choose channels that do not open instantaneously during a spike. Otherwise, a single spike would open all clusters at once. In contrast, for slower channels with a time constant exceeding the width of a single action potential (>10 ms), only a small ratio of the closed channels opens per spike, so that the conductance gradually increases. The channels considered here are even more inert (≈ 100 ms), comparable to slow adaptation currents that require multiple spikes to fully develop (Benda and Herz, 2003). Similarly, clusters of inert channels require multiple spikes in fast succession to switch to the open state (minimum of ≈ 20 Hz, see Figure 4C). In this way, persistent activity at low frequencies cannot open further clusters and remains stable. Similarly, we find that the slow channels make the clusters robust against fluctuation-driven firing (Schreiber et al., 2009) - action potentials triggered by noise in the membrane potential (see Figure 4—figure supplement 2). Finally, the cell is slightly depolarized, so that it spikes in the absence of a stimulus when several clusters are open (see Figure 4C). Even when there are too few open clusters to induce spiking, they still increase the membrane potential and leave a form of ‘silent’ memory (not shown).

Interestingly, the bistable clusters also support an ‘inverse’ form of neuronal memory - persistent activity activated by hyperpolarization and silenced by strong depolarization. For this other form, the cooperative channel dynamics are the same, but they are required to conduct a hyperpolarizing current like for example potassium channels. Then, open clusters act as an additional standing leak current and prevent the cell from firing. Permanently closing the clusters by a hyperpolarizing pulse relieves the cell from this additional leak and allows persistent firing. However, when the cell is stimulated to spike strongly, the clusters return to the open state, restore the standing leak and again silence the cell (see Figure 4—figure supplement 1). Similarly to depolarization-activated persistent activity, the memory relies on the cooperative dynamics of the channels, but how the neuron makes use of this memory changes with their ionic nature.

Cooperative channels mediate graded persistent activity

Graded persistent activity is an intriguing, cell-intrinsic feature of some neurons that has been reported in vitro under carbachol in the entorhinal cortex, the perirhinal cortex and the amygdala (Egorov et al., 2002; Navaroli et al., 2012; Egorov et al., 2006). These neurons fire at high frequency in response to depolarizing current pulses and continue to fire stably when the stimulus is over. The persistent firing rate is lower than the one during the pulse and increases with each new pulse presented. Hyperpolarizing pulses lead to graded decreases in firing rate and finally end the persistent activity, seemingly resetting the cell. While other mechanisms reproducing this firing mode have been suggested in the literature (Loewenstein and Sompolinsky, 2003; Fransén et al., 2006), we here add a novel candidate and show that cooperative clusters allow to capture essential properties of graded persistent activity.

In analogy to the experimental paradigm used by Egorov et al. (2002), our neuron model was subjected to a series of depolarizing current pulses. Figure 5A shows the persistent firing after the first pulse as well as the increase of the frequency after each pulse. The persistent frequency was stable for at least one minute (data not shown). The role of the clusters as the underlying memory variable becomes apparent in the step-like evolution of their conductance, which increases during the pulses and stays constant in between. In this simulation, we observe four distinct stable persistent frequencies (from about 3 Hz to 10 Hz). At 10 Hz, further depolarizing pulses fail to increase the frequency of persistent firing. This saturation is reached, when all clusters have been switched to the open state. In addition, we confirmed that persistent firing can be gradually turned off by hyperpolarizing pulses (Figure 5B).

Figure 5 with 1 supplement see all
Clusters of strongly cooperative ion channels mediate graded persistent activity in a simple neuron model.

(A) Repeated stimulation drives fast spiking in the neuron, which is followed by self-sustained, stable low-frequency activity at increasing rates (top). With each pulse, the driven, fast spiking switches more clusters to the open state and adds further stable conductances (bottom). The built-up conductance allows a persistent current to flow and sustains activity beyond stimulation. (B) In the same manner, strong hyperpolarization closes the clusters and allows to reduce the frequency of persistent activity until finally the cells stops firing. The simulation procedure is described in Materials and methods and parameters are summarized in Table 1.

We further investigated which other channel properties - next to strong cooperativity - are essential for the graded nature of persistent activity (Figure 5—figure supplement 1). We concentrated on the time constant and the conductance of the channels. Fast channels, as discussed in the previous section, are less suited than slow channels. The slow channels ensure that clusters gradually open during stimulated fast spiking and halt during slow spiking, so that the persistent firing remains stable. The conductance of the channels has to be at an appropriate intermediate level: sufficient to drive the neuron over threshold, but limited to provide persistent firing frequencies of maximal 10–15 Hz. Otherwise, faster spiking can open further clusters and the persistent activity increases, so that it is no longer stable.

Concluding, a small set of cooperative clusters of ion channels - when combined with ‘normal’, non-cooperative potassium and sodium channels (responsible for the generation of action potentials) - can reproduce the computationally relevant feature of graded persistent activity.

Cooperative channels induce a switch to mnenomic firing in perirhinal cortex neurons

Last but not least, we test the cooperativity mechanism experimentally. Our simulations indicate that clusters of cooperative ion channels can enable the mnemonic firing mode of graded persistent activity in a simple point neuron. In biological neurons, however, other ionic currents, morphological effects or noise could interfere with the action of the cooperative channels and undermine their role as memory units. Thus, we employed the dynamic clamp technique (Sharp et al., 1993; Prinz et al., 2004) to equip perirhinal cortex (PR) neurons in-vitro with the same cooperative clusters used in our models, acting as ‘virtual channels’ in real cells (Figure 6). PR neurons have been shown to exhibit graded persistent activity under activation of muscarinic acetylcholine receptors (Navaroli et al., 2012). Consequently, we applied no carbachol, so that any form of persistent activity observed in the recorded PR neurons stems from the clusters of cooperative channels introduced via dynamic clamp. Accordingly, all recordings were performed in the presence of a synaptic blocker to exclude network contributions. For details of experiments, see Materials and methods.

Figure 6 with 1 supplement see all
Dynamic clamp experiment.

(A) Intracellular recording of a perirhinal cortex neuron in dynamic clamp mode. The computer simulates the state of clusters of cooperative channels given the measured membrane potentials and emulates the effect of their conductance to the neuron’s activity by inserting a corresponding current. (B) Neurolucida drawing of the recorded neuron. (C) Graded persistent activity of the recorded neuron mediated by the multi stable conductance of the clusters of cooperative channels. (D) Hyperpolarization brings the cell back to rest via intermediate levels of persistent firing. In control recordings with identical, yet- independent channels, no persistent firing was observed (Figure 6—figure supplement 1). Slice preparation, electrophysiological setup and the dynamic clamp emulation of the clusters are described in Material and methods. Cluster parameters are summarized in Table 1.

The dynamic clamp technique allowed us to control the cluster size, channel interactions and channel kinetics of the artificial conductances during the experiment. Like in the model neuron simulations, we adapted the cluster parameters so that the bistability range fitted to the properties of the recorded cell.The resting potential determined the center of the bistable range, whereas the spike shape set boundaries for the bistable range. Voltage values during afterhyperpolarization lied above the lower boundary (to prevent reset by a spike’s afterhyperpolarization), but the upstroke of the action potential exceeded the upper boundary of the bistable range.

Stable recordings with cooperative channel dynamics added via dynamic clamp were obtained for three PR neurons. Under cooperative conditions, all cells exhibited persistent activity in response to current stimulation with pulses (Figure 6C). Persisting activity following pulses was graded, that is the persistent firing rate increased with each presented pulse. As in the mathematical model, hyperpolarizing pulses lowered the rate of persistent activity and when applied multiple times, silenced it (Figure 6D). In the control conditions (identical, yet non-cooperative channel clusters), no persistent firing could be evoked (Figure 6—figure supplement 1).

Periods of stimulation were followed by a slow afterhyperpolarization (sAHP) in both control and cooperative conditions, which for paradigms with cooperative gating generated a 3–4 s lasting refractory period between the stimulation and the self-sustained firing. Interestingly (and in agreement with mathematical modelling), the persistent activity ‘unfolded’ despite these intrinsically generated dynamics. We note that such sAHP were not expected in the original recordings of carbachol-induced persistent activity (Egorov et al., 2002), as carbachol blocks sAHPs.

Taken together, the experimental data support the hypothesis that cooperative ion channels can mediate cellular persistent activity. The persistent dynamics are very robust and could be easily evoked without further knowledge of the intrinsic properties of the recorded cells, such as the multitude of other ionic currents, complex cell morphology, and the presence of channel noise.

Discussion

Despite a large body of evidence for cooperative interactions between channels pivotal to the nervous system (Iwasa et al., 1986; Dekker and Yellen, 2006; Kim et al., 2014; Choi, 2014; Dixon et al., 2015; Moreno et al., 2016; Clatot et al., 2017), it is unknown which function coupled channels have for neural dynamics and computation. Based on a mathematical analysis, we show in this study that a cluster of cooperative channels can gate with hysteresis and find, in both simulations and experiments, that multiple such clusters embedded in the membrane of a neuron mediate mnemonic firing like persistent activity. Therefore, we suggest that ion channel cooperativity might serve as a cell-intrinsic memory mechanism at the voltage level.

From synaptic learning, adaption of the immune system to epigenetics, many biological memory systems are founded on strong positive autofeedback (Lisman, 1985; Burrill and Silver, 2010). Feedback loops are ubiquitous in biology (Thomas and D’Ari, 1990); they amplify signals as in calcium-induced calcium release and enhance sensitivity in controlling protein function by multisite phosphorylation. In particular, strong feedback enables bistability and memory as in networks of recurrently connected neurons (Amit, 1990; Durstewitz et al., 2000). Here, we exploit the same principles for voltage-gated ion channels, where cooperativity can act as a strong auto-feedback and enables clusters of channels to act as bistable macrochannels with a memory of previous voltage levels.

Bistability and memory of a cooperative channel cluster

Our gating analysis shows that clusters only become bistable, when channels cooperate strongly. We suggest two criteria to detect this strong coupling regime: a bimodal distribution of cluster conductance states and a hysteresis in gating. Along these lines, experimental evidence for the strong coupling regime comes from cooperative calcium channels in the heart (Navedo et al., 2010). Although these channels couple weakly for the majority of clusters, a few clusters are bistable, switching between the all open and the all closed state, which implies the presence of strong coupling.

Our study of channel noise confirms that also small clusters exhibit a robust hysteresis and can be used for neuronal memory. Cooperativity is inherently local - only close-by channels can interact - and is therefore restricted to small assemblies of channels (Gutkin and Ermentrout, 2006). In these small clusters, spontaneous channel fluctuations become a prominent source of noise, a phenomenon occluded by the assumption of all-to-all coupling in previous studies (Naundorf et al., 2006; Zarubin et al., 2012). In contrast, the clustered model here captures channel noise and therefore predicts that the lifetime of the open and closed cluster state decreases in smaller clusters. However, we find that experimentally observed cluster sizes of about 10 channels (Moreno et al., 2016) - each gating in the millisecond range - suffice for memory on the timescale of seconds. Consequently, also small clusters offer a robust memory and, in contrast to one large cluster with only two states, the ensemble of small clusters yields a more versatile response and has a larger memory capacity.

Ion channel cooperativity as a mechanism for graded persistent activity

Neurons can exploit these hysteretic cluster macrochannels as a short-term, analogue memory of recent voltage history - spikes open clusters and strong hyperpolarization closes them. In particular, the ensemble of clusters can differentiate stimulation strength - very strong or prolonged firing opens more clusters - so that the level of persistent firing depends on the stimulation history. In this way, the collection of cooperative clusters offers the memory capacity required for forms of graded persistent activity previously observed in neurons from entorhinal (Egorov et al., 2002), perirhinal (Navaroli et al., 2012) and prefrontal cortices (Winograd et al., 2008), as well as the amygdala (Egorov et al., 2006).

Our analysis predicts the following prerequisites for persistent activity. First, the cooperative channels should have a bistable range which comprises the resting membrane potential and whose borders set the voltage activity needed to switch the persistent activity on and off (see Figure 3 and Figure 4). In this way, as our dynamic clamp experiments in the perirhinal cortex showed, the memory is robust for multiple seconds and can be reliably controlled via neural activity. Second, for graded persistent activity, the channels should have a slow time constant, on the order 50–100 ms matching the maximal persistent firing frequency. Slow cooperative channels enable a gradual opening of individual clusters during stimulation, prevent progressing excitation during the persistent period and make the clusters more robust against fluctuation-driven spikes. If the time constant was too fast, all clusters would open with the voltage elevation of the first action potential and, consequently, only one persistent response state (i.e. all clusters open) would be possible.

The type of ion conducted by cooperative channels, however, is less constrained. In principle, persistent activity can be mediated by both depolarizing and hyperpolarizing cooperative channels like calcium or potassium channels, respectively (see Figure 4—figure supplement 1). If clusters consist of depolarizing calcium channels, persistent activity starts after strong spiking as observed by Egorov in the entorhinal cortex (Egorov et al., 2002). If on the other hand, clusters consist of hyperpolarizing potassium channels, persistent activity starts after strong hyperpolarization - similar to hyperpolarization-activated GPA as reported by Winograd in prefrontal cortices (Winograd et al., 2008). Independent of the ion type, we expect that the cooperative channels are separate from the action potential-generating sodium and potassium channels, or at least form only a small subgroup therein. This separation would prevent memory from interfering with the action potential, the first based on persistent currents and the latter on regenerative, memoryless currents.

Experimental evidence and comparison to other model of persistent activity

To date, it is not known, whether neurons exhibiting graded persistent activity express cooperative ion channels. However, in mammalian cells from the hippocampus, a cooperative variant of CaV1.3, a calcium channel wide spread in the brain, has been demonstrated to provide persistent depolarizing currents and increase the firing rates of hippocampal neurons (Moreno et al., 2016). Other candidates implicated in persistent activity, the transient receptor potential cation (TRPC) channels (Zhang et al., 2011), are known to cluster and therefore provide the spatial proximity required for cooperativity (Nilius and Owsianik, 2011). However, a recent study showed that graded persistent activity in the entorhinal cortex of mice does not require TRPC channels (Egorov et al., 2019). Generally, the channels required for persistent activity are still under debate. As our work demonstrates, only a small fraction of channels has to be cooperative to produce persistent activity, so that they might be easily overlooked, especially if cooperativity depends on other intracellular regulators like calmodulin (Moreno et al., 2016).

As a model of persistent activity, cooperativity shares positive feedback as the core principle with other hypothesis like calcium modulated conductances, but differs in several aspects. The central difference is that coupled channels have an inherent, direct feedback mechanism - one channel opening facilitates another - whereas independent channels require an indirect interaction. Accordingly, previous studies suggested a feedback cycle of spiking, calcium inflow and a calcium modulated conductance as the basis of persistent spiking (Rodriguez et al., 2018). A direct consequence is that memory in the case of calcium modulated channels requires spiking, whereas cooperative clusters could also implement a silent memory - a long lasting change of excitability after stimulation. However, silent memories are also possible in extended models of calcium modulation, where a rise in calcium triggers persistent conductance changes through intracellular signaling (Fransén et al., 2006; Winograd et al., 2008). Another difference is that models based on calcium modulation naturally address the finding that blocking calcium channels or strong calcium buffering prevents persistent activity (Egorov et al., 2002; Cui and Strowbridge, 2018). For the cooperativity hypothesis, this central role of calcium could imply a calcium regulated coupling, like in the case of calcified calmodulin mediated channel interactions (Moreno et al., 2016).

As mediators of graded persistent activity, clusters of cooperative channels could be the computational substrate of important cognitive processes such as short term memory and evidence accumulation (Zylberberg and Strowbridge, 2017). Beyond graded persistent activity, the cluster bistability may serve other plasticity mechanism of intrinsic excitability such as dendritic attenuation or boosting of synaptic inputs (Marder et al., 1996; Debanne et al., 2019).

Cooperativity dynamically regulates the computational repertoire of ion channels

An attractive feature of cooperative channels is the computational flexibility. For a cell-intrinsic memory on demand, biophysical modulators could turn on and off the cluster bistability. Bistability only emerges, when the channels are coupled strongly and the cluster consists of a sufficient number of channels. Thus, one regulatory mechanism could be to change the coupling as in the case of calcium regulated cooperativity (Moreno et al., 2016). Another regulatory mechanism could be to control clustering of channels, which is subject to factors like neural activity (Misonou et al., 2004), extracellular pH (Sumino et al., 2014), ionic concentrations (Eisenach et al., 2014) and lipid signaling (Hilgemann et al., 2018). As a composite of multiple channels, a cluster macrochannel is a more flexible conductance than its ‘hard-coded’ parts.

The emergence of memory in a cluster of cooperative, but memoryless channels suggests a general role of ion channel cooperativity. Cooperative interactions guide the formation of novel macrochannels with a gating repertoire absent at the single channel level. At the level of the neuron, the common membrane potential orchestrates sodium and potassium channel gating to generate the spike. In the same way, cooperativity could orchestrate gating in small channel assemblies to enrich neural dynamics. Correspondingly, it has been suggested that different TRP channels from heteromultimeres to create a wide variety of functions (Nilius and Owsianik, 2011). Another example of a heterogenous cluster is the assembly of BK and CaV, 1.3 channels (Berkefeld et al., 2006; Vivas et al., 2017). A quantitative understanding of these channel complexes requires detailed experimental characterization of the channel couplings like in Sato et al. (2018). Still, simple coupling models like the one presented here already reveal the potential of cooperativity to provide emergent gating functions.

In summary, clusters of cooperative channels broaden the computational repertoire of neurons. Extrapolating from the current study on cell-intrinsic memory, ion channel cooperativity can mediate direct feedback loops between channels and therefore could allow to form macrochannels with novel gating dynamics. If, additionally, neuromodulation can control cooperative interactions, ion channel cooperativity provides an extremely versatile cell-intrinsic mechanism to enrich and regulate neural activity.

Materials and methods

Neuron and channel model

Isolated ion channel model

Request a detailed protocol

We assume that the cooperative channels have a single activation gate and model their gating dynamics in isolation according to the calcium channel dynamics from the Morris-Lecar model (Morris and Lecar, 1981).

The activation kinetics are

(1) dmdt=m(V)-mτ(V)

where the steady state activation is

(2) m(V)=12(1+tanh(V-V1/2k)),

and the activation time constant reads

(3) τ(V)=τcosh-1(V-Vmσ).

Correspondingly, the channel has two states, open (O) and closed (C) and their kinetics reads

Cβ(v)α(v)O,

with opening rate

α(V)=m(V)τ(V),

and closing rate

β(V)=1-m(V)τ(V).

The channels have a single channel conductance gcoop and reversal potential Ecoop. Original parameters are Ecoop=100mV, V1/2=-1mV, k=15mV, τ=0.05ms, Vm=-1mV and σ=30mV, see Zarubin et al. (2012). Unless reported otherwise, we choose a single channel conductance of gcoop=2.5pS as reported for calcium channels (Church and Stanley, 1996). Modifications of these parameters are summarized in Table 1.

Cooperativity model

Request a detailed protocol

In order to capture cooperative interactions among channels, we model the activation of a channel as dependent on both the membrane potential and the state of near-by channels (Naundorf et al., 2006). Specifically, we assume that the activation m~ of a channel among o open neighbours has the form

(4) m~(V,o)=m(V+oj).

where j is the coupling strength between two channels. Therefore, if the coupling is positive j>0, opening of a neighbouring channel shifts the activation curve towards lower membrane potentials. As a result, a channel with open neighbours has itself an increased open probability, see Figure 1A.

For the time constant, we posit that the same shift applies,

(5) τ~(V,o)=τ(V+oj).

Cluster model

Request a detailed protocol

Next, we introduce the experimentally observed clustering of the channels. We presume that only channels in the same cluster are close enough to cooperate, so that channels in different clusters gate independently. For a cluster, a characteristic measure of cooperativity is the maximal shift J, which corresponds to the shift of the activation curve when all neighbours of a channel are open. Under the simplifying assumption of constant coupling strength j among all channels, the maximal shift in a cluster of size S amounts to J=(S-1)j.

A note on comparing coupling strengths with previous studies: We choose J for the maximal shift in accordance with the notation of Zarubin et al. (2012). In contrast, Naundorf et al used J to denote the coupling strength between two channels, which in the work presented here is j (Naundorf et al., 2006).

In a cluster of size S with o open channels, all S-o closed channels open with rate

α~(V,o)=α(V+oj)

and all o open channels, having only o-1 open neighbours, close with rate

β~(V,o-1)=β(V+(o-1)j).

Instead of tracking the cluster state in terms of each constituting channel, we can also view the cluster as a macrochannel. For S two-state channels, this macrochannel has S+1 conductance states, from all channels closed to all channels open. In a small time interval, at most one channel opens or closes, so that transitions are restricted to adjacent cluster states with a difference of one open channel,

oβo+1,o(v)αo,o+1(v)o+1.

The macrochannel rates result from the number of channels that are possible candidates for the transition and the single channel transition rate. Therefore, in a cluster with o open channels, the opening of one of the other S-o closed channels happens with rate

(6) αo,o+1(V)=(S-o)α(V+oj).

Correspondingly, from a state with o+1 open channels, the closing of one of the open channels each having o open neighbours reads

(7) βo+1,o(V)=(o+1)β(V+oj).

Figure 1C shows simulated traces of clusters and demonstrates their behaviour as macrochannels.

Neuron model

Request a detailed protocol

We used a conductance-based neuron model with a single isopotential compartment,

(8) CdVdt=Iapp-Icluster-IV,

where Icluster is the current through the clusters of cooperative channels, IV summarizes the other channel currents, Iapp is a stimulus current and C is the capacitance of the membrane.

Cluster current

Request a detailed protocol

We model the cooperative channels arranged in N clusters, each composed of S identical channels. Then, the current through all clusters is determined by the total number of open channels Ocoop among all clusters,

(9) Icluster=gcoopOcoop(V-Ecoop),

where Ecoop denotes the reversal potential of the considered ions. For most of the article, we assume that the cooperative channels conduct a depolarising current with Ecoop=100mV (e.g. Ca2+ or Na+). In Figure 4—figure supplement 1, we then consider the case of a hyperpolarising current with Ecoop=-100mV (e.g. K+).

In a description, where each channels is tracked, Ocoop simply counts the number of channels in the open state. In the alternative macrochannel description of the clusters, the total number of open channels is obtained from the number of macrochannels in the different conductance states, so

Ocoop=o=0Sγoo,

where γo is the number of clusters with o open channels. This occupancy vector γ can also be used to capture the distribution of conductance states of a single cluster over time (Figure 1C).

As opposed to the other ionic conductances, the dynamics of the clusters is modeled on the level of the underlying jump process generated by the single-channel gating events. Such a detailed description is necessary to account for the fact that cooperative interactions reduce the number of independent stochastic units and therefore increase fluctuations (1/N as opposed to 1/SN White et al., 2000). Effectively, the N clusters represent the independent units, so that with around 100 clusters, a jump process description is adequate to account for the discrete nature of the fluctuations.

Leak and action potential mediating currents

We choose a type 1 neuron model to account for the continuous low-frequency firing range as observed in graded persistent activity (Egorov et al., 2002). The model presented is a version of the Traub-Miles model (Benda, 2002), which comprises action potential generating potassium and sodium currents and a leak current,

IV=A(INa+IK+IL).

The original model is formulated with conductance densities and is independent of the neuron surface area A. In the present model, the cluster conductance is defined in absolute terms, because it stems from a concrete number of channels with a fixed conductance. Correspondingly, we turn the remaining currents and the capacitance into absolute quantities by choosing a surface area, for details see Choosing a neuron surface. The capacitance density is C¯=1μF/cm2.

In the following, all gating variables have first order kinetics of the form

(10) x˙=αx(V)(1-x)-βx(V)x.

Sodium current

Request a detailed protocol
INa=g¯Nam3h(V-ENa)

with g¯Na=100mScm2, ENa=48mV and gating variables m and h with rates

αm(V)=0.32kHz(V/mV+54)/(1-exp(-0.25(V/mV+54))),
βm(V)=0.28kHz(V/mV+27)/(exp(0.2(V/mV+27))-1),
αh(V)=0.128kHzexp(-(V/mV+50)/18),
βh(V)=4.0kHz/(exp(-0.2(V/mV+27))+1).

Potassium current

Request a detailed protocol
IK=g¯Kn4(V-EK)

with g¯K=200mScm2, EK=-82mV and gating variable n with rates

αn(V)=0.032kHz(V/mV+52)/(1-exp(-0.2(V/mV+52))),
βn(V)=0.5kHzexp(-(V/mV+57)/40).

Leak current

Request a detailed protocol
IL=g¯L(V-EL)

with g¯L=0.1mScm2 and EL=-67mV.

Choosing a neuron surface

Request a detailed protocol

When we fix the peak conductances of the other ionic currents, the contribution of the cluster current for a fixed number of cooperative channels depend on the neuron area A. For example, for a small neuron area, the cluster conductance becomes relatively larger and could drive faster persistent spiking than for a large neuron area. In choosing a neuron area, we try to meet two characteristics of graded persistent activity, namely its low-frequency range and its quasi-continuous nature.

Specifically, the frequency of persistent firing is usually located in the low-frequency range below 15 Hz, which sets a limit to the conductance of the clusters. In terms of conductance densities, we find that close to rheobase Iapp=0.105μA/cm2 and with a reversal potential Ecoop=100mV common for Ca2+ channels, a cooperative channel conductance density of about g¯coop0.0004mScm2 is sufficient to drive spiking at about 10 Hz.

Furthermore, the levels of persistent firing are assumed to be quasi-continuous, which requires a large number of clusters with a small conductance. In all simulations, we assumed about 100 clusters, which in principle allow 100 levels of persistent firing. In a frequency range of up to 10 Hz, this would correspond roughly to a frequency resolution of 0.1 Hz, which lies below the experimentally observed grading of about 0.5 Hz (Fransén et al., 2006). Furthermore, the clusters have to be of a certain size to be bistable on the timescales of seconds, as shown in Figure 3. Hence, for a high number of clusters N = 100 with size S = 8 and a typical single channel conductance g=2.5pS, the total cooperative channel conductance is gcoop=2μS.

Therefore, we arrive at a required surface area of,

A=gcoopg¯coop=0.005cm2,

which is a large neuron area compared to typical neuron surface areas as reported in literature (e.g. Ambros-Ingerson and Holmes, 2005). This large membrane area also explains the rather high currents that are needed to excite and hyperpolarise the neuron. In principle, a smaller neuron area could be reached by smaller cluster numbers, smaller cluster sizes, a lower reversal potential or an increase of the firing range.

Stimulation with a white noise current

Request a detailed protocol

We investigate how robust the clusters are against noise in the membrane potential (Schreiber et al., 2009). In particular, we mimic fluctuating synaptic input to the neuron, a major noise source, by injecting a white noise current with standard deviation σI and time resolution dtnoise,

Inoise(t)=kik(H(t-kdtnoise)-H(t-(k+1)dtnoise));ik from 𝒩(0,σI2),

where H(x) is the Heaviside step function and 𝒩(0,σI2) denotes a normal distribution with zero mean and variance σI2.

We choose a time resolution of dtnoise=0.5ms, oriented at the temporal width of synaptic currents, and then vary the noise intensity (see Figure 4—figure supplement 2). Note that the simulation step size is much smaller than the temporal resolution of the noise (see Simulation of the neuron model).

Overview of cluster and cooperative channel parameters.

Simulation and analysis of cluster dynamics

Mean channel activation in a cluster

Request a detailed protocol

In a large cluster of cooperative channels, the number of open channels is expected to coincide with its mean o¯=mcoop(V)S, where mcoop(V) is the average activation of the channels in the cooperative ensemble. Moreover, the average activation reflects how much the open neighbours shift the activation curve,

mcoop(V)=m(V+o¯j),

where we use Equation 4 for the cooperativity altered activation and m(V) is the single channel activation function. A rewrite of the shift term o¯j=mcoop(V)Sjmcoop(V)J results in a self-consistency relation (Naundorf et al., 2006, Supplementary Notes 2),

(11) mcoop(V)=m(V+mcoop(V)J).

We numerically solve Equation 11 for a range of membrane potentials to obtain the activation curves mcoop(V) for different couplings J. As discussed in the results on cluster dynamics, the activation becomes bistable for sufficiently strong coupling J>Jcrit (Figure 1).

The critical coupling strength depends on the form of the single channel activation m(V). For the activation curve in Equation 4, it reads (derivation below)

Jcrit=2k.

Hence, the critical coupling Jcrit coincides with the width of the activation curve. Put differently, bistability emerges, when cooperative facilitation by the neighbouring open channels can keep a channel open despite a low membrane potential at which the channel would usually be closed. Thus, when all neighbours are open, the shift has to exceed the width of a channel’s activation curve (Figure 1A).

Derivation of the critical coupling strength

Request a detailed protocol

For the derivation of the critical coupling strength, we follow the argument of Huang et al. (2012), (Equations 4-10). They consider cooperative sodium channels and are interested in the coupling strength, where the channel activation becomes a step function. This coupling strength in principle also induces bistability; the steps in the activation curve occur at the edges of the bistable range. However, sodium channels transit to the inactivated state when they are open and thus they cannot be bistable.

If a bistable range exists, the number of solutions of the self-consistency Equation 11 has to change from one to three (bistable) and back to one again, when the voltage is increased. From this observation, we can deduce a condition in the coupling strength J. Namely, each solution mc corresponds to an intersection of the left and right hand of Equation 11,

mc=m(V+mcJ)

and for the voltages where the number of solutions change, these intersections have to be tangential, that is

ddmcmc=ddmcm(V+mcJ)
1=J2k(1(2m(V+mcJ)1)2).

Here, we use that dm(V)dV=12k(1-(2m(V)-1)2). Again using the self-consistency relation mc=m(V+mcJ), we observe that the solution has to obey

Jkmc2-Jkmc+12=0.

Finally, as the solution has to be real, bistability sets the following condition on J:(Jk)22JK0J2k=Jcrit.

Mean first passage times between cluster states

In order to study the stability of the open and closed cluster state, we calculate the mean first passage times between these two states. We employ the macrochannel description, where for a cluster of size S the dynamics form a continuous-time Markov process with transition matrix

Qij(V)=δi,j+1αi-1,i(V)+δi,j-1βi+1,i(V); 0i,jS.

Qij(V) denotes the voltage-dependent transition rate from cluster state j to i and summarizes the opening and closing rates from Equation 6 and 7.

Over a small time interval, the transition rates become transition probabilities and the continuous Markov chain can be discretized. That is, choosing a small time interval Δt, such that for all i,j it holds QijΔt1, the transition probabilities are Pij=QijΔt for ij and Pii=1-jiPji. For a discrete Markov chain, the mean first passage steps K (with Kij denoting the mean number of steps required to reach state i from state j) ) can be obtained by solving the system of equations given by Allen (2010), see p.69)

K=E+(K-diag(K))P,

where E is a (S+1)×(S+1) matrix of ones. Multiplication with the time interval recovers the mean first passage times M=KΔt.

For the stability of a bistable cluster, the mean first passage times between the open and closed state are of particular interest: the average time τ¯OC it takes until a cluster spontaneously switches from the state with all channels open to the one with all channels closed or vice versa τ¯CO. In terms of M, they read τ¯OC=M0,S and τ¯CO=MS,0. Figure 3A depicts the mean residence times as a function of the membrane voltage.

Numerical simulation of cluster dynamics

Request a detailed protocol

For the simulation of the cluster dynamics at a constant membrane potential, we use both a fixed time step method and the Gillespie algorithm (Gillespie, 1977).

In the fixed time step method, we track the state of each channel in the cluster. For each time step, we look up the number of open channels, evaluate the corresponding opening and closing rates and obtain the transition probabilites by multiplying the rates with the time step. Finally, a random number generator is used to update the state of the channels according to the transition probabilities. This method is used for the cluster simulations in Figure 1 and Figure 2.

For the application of the Gillespie algorithm (Gillespie, 1977), we switch to the macrochannel description, where the states only specify the number of open channels, but ignore the configuration of the individual channels. Specifically, for clusters of size S, the S+1 conductance states represent the species in the Gillespie algorithm, whereas the 2S transitions between adjacent cluster states with their respective rates form the reactions.

Simulation of the neuron model

In the neuron model, the coupled dynamics of the continuous variables, voltage (Equation 8) and gating variables (Equation 10), and the jump process of the cooperative channels form a hybrid stochastic system (Anderson et al., 2015). For the simulation of spiking activity, we use brute, a custom written fixed time step algorithm.

BRUTE algorithm

Request a detailed protocol

The brute algorithm tracks the state of each channel and uses a fixed time method to update their states. In order to take into account the cooperative interactions among channels in one cluster, it considers each cluster separately and calculates the opening and closing rate of the channels therein. For N clusters of size S, the channel states can be cast into a matrix C{0,1}N×S, which the algorithm updates in each time step.

In the following, the membrane potential and the continuous gating variables are summarized in 𝐱 and the current time is t. Channel states are encoded as 0 (closed) or 1 (open). Then, in pseudo code, the update algorithm with time step Δt is

  1. Initialize: set t=t0, set the initial membrane potential and gating variables 𝐱0 and set the initial channel states C{0,1}N×S

  2. Generate an array of random numbers r[0,1]N×S

  3. Numerically integrate voltage and gating variables dynamics (Equation 8; 10) until t+Δt

  4. Update each cluster, that is each row i of C as follows

    1. calculate the number of open channels oi=jCij

    2. update each channel Cij

      • if Cij=1 and rij<β(V,oi-1)Δt then Cij=0

      • else if Cij=0 and rij<α(V,o)Δt then Cij=1

  5. Increment t by Δt and check if simulation is finished, otherwise continue with step 2

In all simulations, we chose a time step Δt=0.25us. At this small time resolution, we observed a convergence of the cluster dynamics. A python implementation of the algorithm with a simple example is available at https://itbgit.biologie.hu-berlin.de/cooperativity/brutelib.git (copy archived at https://github.com/elifesciences-publications/49974-brutelib).

Dynamic clamp experiment

Preparations of mouse brain slices

Request a detailed protocol

Dynamic clamp experiments were conducted with horizontal mouse brain slices (300–450 µm thick) containing the hippocampus, entorhinal and perirhinal cortices. Brain slices were obtained from male C57BL/6N mice (4–8 weeks old, 20–25 g) using standard proceedings (Roth et al., 2016). Mice were purchased from Charles River Laboratories (Sulzfeld, Germany, Strain Code: 027) and were taken care of in the Interfaculty Biomedical Research Facility in Heidelberg. Housing was provided in Makrolon II cages with a maximum of three animals and tissue nesting material made of cellulose. Animals had ad libitum access to food and water. All experimental protocols were conducted in compliance with German law and with the approval of the state government of Baden-Württemberg (Project T100/15).

In order to minimize the stress of euthanasia, mice were sedated by exposure to CO2 in a rising concentration (20–30 l/h sourced from a compressed gas cylinder) until the animal fell unconscious. Mice were subsequently killed by decapitation and the brain was quickly removed and transferred to 4°C cold, carbogen buffered (95% O2, 5% CO2 at pH 7.4) artificial cerebrospinal fluid (ACSF) containing the following (in mM): 124 NaCl, 3 KCl, 1.8 MgSO4, 1.6 CaCl2, 10 glucose, 1.25 NaH2PO4, 26 NaH2CO3. Brain slices were cut using a vibratome (Leica VT1200S, Nussloch, Germany). Then, slices were transferred to a Haas-type interface chamber (Haas et al., 1979), perfused with ACSF at a rate of 1.2–1.4 ml/min at 34 ±1°C. Slices rested for at least 2 hr before electrophysiological recordings.

Electrophysiological recordings

Request a detailed protocol

Single-cell recordings were obtained from principal neurons of the perirhinal cortex layer II. These neurons show graded persistent firing under activation of muscarinic cholinergic receptors (Navaroli et al., 2012).

The recordings were obtained with sharp microelectrodes (tip resistance 100–130 MΩ) pulled from 1.0 mm borosilicate glass capillaries (Harvard Apparatus, Cambridge, UK, Cat. No. 30–0019) on a DMZ Universal Electrode Puller (Zeitz, Martinsried, Germany) and filled with 2 M K-acetate that in some cases contained 1% biocytin.

Recorded signals were low-pass filtered at 10 kHz, amplified x10 using an SEC-05X amplifier (npi electronic, Tamm, Germany) and digitized at 20 kHz with an analog to digital converter (MICRO 1401 mkII ADC, CED, Cambridge, UK). Signals were visualized and saved using Spike2 software (CED, Cambridge, UK). All intracellular recordings were conducted in bridge mode and bridge balance was monitored and adjusted during experiments.

Positive and negative current pulses (duration and amplitudes controlled by Spike2 Software and SEC-05X amplifier) were applied via the recording electrode to determine input resistance and firing properties. Input resistance was estimated during the experiment and later calculated off-line using Matlab software (The Mathworks, Natick, MA). Resting membrane potential was calculated by subtracting the potential offset after withdrawal from the cell at the end of the recording.

All cells selected for measurements had an input resistance above 20 MΩ, a stable membrane potential throughout the recording and exhibited firing properties of principal neurons. Membrane potential of neurons was controlled by manually adjustable DC current injection through the recording electrode to determine the current needed for intended firing frequencies and depolarisation near threshold.

To prevent spontaneous network activity in slices, all experiments were performed in the presence of ionotropic glutamate receptor blockers 6-cyano-7-nitroquinoxaline-2,3-dione (CNQX, 30 µM, Cat. No. 1045) and DL-2-amino-5-phosphonovaleric acid (APV, 10 µM, Cat. No. 0105) obtained from Tocris (Bristol, UK). Drugs were bath applied by continuous perfusion. Biocytin (Cat. No. B4261) was obtained from Sigma-Aldrich (Taufkirchen, Germany).

Biocytin staining

After recordings the slice was fixed in 4% PFA in phosphate buffer (PB) at 4°C for 90 min and then stored in phosphate buffered saline (PBS) at 4°C until further processing. It was washed 3 × 15 min to remove excess PFA and incubated with Avidin - Alexa Fluor 488 conjugate (1:1000; Life Technologies, Carlsbad, CA) diluted in PBS containing 5% normal goat serum and 0.2% Triton X at room temperature under light protection for 2 hr. The slice was then washed in PBS for 15 min, incubated with 4,6-diamidino-2-phenylindole (DAPI; 1:10,000; Carl Roth, Germany) in H2O for 3 min, washed again in PBS (15 min) all at room temperature and afterwards embedded in Mowiol 4–88 (Sigma-Aldrich, Taufkichen, Germany).

Fluorescence images were acquired with a Nikon A1+ Confocal Microscope (Nikon, Düsseldorf, Germany) and reconstruction of the neuron was obtained from z-stack confocal images using Neurolucida tracing software (MBF Bioscience, Williston, VT).

Dynamic clamp

Request a detailed protocol

Data acquisition and dynamic clamp loop were controlled by RELACS, V0.9.8, RRID:SCR_017280. The feedback loop run at a frequency of 20 kHz and consisted of sampling the membrane potential, updating the state of the clusters of cooperative channels and injection of the corresponding current (Equation 9). The software allowed online adjustments of the following model parameters:

  • the channel to channel coupling strength j

  • cluster number N and cluster size S

  • single channel conductance gcoop and reversal potential Ecoop

  • parameters of the channel kinetics, namely maximum of the time constant τ and the half-activation voltage V1/2

Real time cluster update in the dynamic clamp loop

In order to meet the real time requirement imposed by dynamic clamp, we employ the macrochannel description of the clusters (see Cluster model). In each time step, the macrochannel description updates the cluster population in the S+1 conductance states of a cluster by evaluating the 2S possible transitions. This provides a huge reduction of transitions, when compared to the NS state update required to each channel.

In pseudo code, the update algorithm works as follows:

  1. Initialize: set initial cluster population γ0{0,,N}S+1

  2. Get currently measured membrane potential V

  3. Generate two arrays of random numbers ropen[0,1]S and rclose[0,1]S for the 2S opening and closing reactions

  4. Calculate propensities of the 2S reactions via the macrochannel transition rates in Equation 6 and Equation 7.

    • λoopen=γoαo,o+1(V) for o{0,S-1}

    • λoclose=γo+1βo+1,o(V) for o{0,S-1}

  5. Update cluster populations

    • S opening reactions: if roopen<λoopenΔt then increase γo+1 by one and decrease γo by 1

    • S closing reactions: if roclose<λocloseΔt then increase γo by one and decrease γo+1 by 1

  6. Wait until the time interval Δt has passed and continue at step 2

The successive update of the cluster populations can lead to a problem for states, which have a population of one and where both the opening and the closing reaction happen. In our solution, we chose to assure a positive population for all states and reject the closing reaction in such a case. In an offline comparison with the exact Gillespie algorithm, we tested that our update algorithm accurately captured the evolution of the clusters despite this slight bias toward opening reactions.

The algorithm with instructions how to use it in the recording software RELACS is available at https://itbgit.biologie.hu-berlin.de/cooperativity/dynamic_clamp_model (copy archived at https://github.com/elifesciences-publications/49974-dynamic_clamp_model).

Experimental protocol

Request a detailed protocol

In all recordings, we first measured the f-I curve of the neuron and determined the resting potential and the form of the action potentials. With these neuron properties at hand, we chose the model parameters for the dynamic clamp experiment. First, we aimed at a total conductance of the cooperative channels, that allowed a current to flow sufficient for about 10 Hz persistent firing. Correspondingly, we chose the number of cooperative channels, given that channel conductance and reversal potential were orientated at values from calcium channels. Second, we split the channels into clusters such that we obtained both a large number of clusters and a sufficient cluster size for long-term stability. Third, we adjusted the bistable regime of the clusters such that its center coincided with the resting potential and its borders remain far below the spike peak. To this end, we selected the appropriate coupling strength and half-activation voltage of the channels.

We used the stimulation protocol of consecutive pulses, either depolarising or hyperpolarising, to test for graded persistent activity. Pulses had a length of 4–6 s with a period of about 25 s. Often, we saw channels opening persistently during the depolarising pulses, but no persistent activity. In these cases, we would increase the baseline current to bring the neurons closer to threshold. Additionally, we run control protocols with zero coupling to test whether the cell displayed no persistent activity in the absence of cooperativity (Figure 6—figure supplement 1).

Analysis

First, we categorised the recordings according to the amplitude and the number of pulses, for example UP-4 for a series of four depolarising pulses or DOWN-1 for one hyperpolarising pulse. Then, the inter spike intervals were analysed during and after the pulses as well as before and after the pulse protocol. For the inter spike intervals, we excluded the first 1000 ms of each interval to account for transient currents. Next, we defined a simple success criteria for a protocol: during an UP (DOWN) protocol, the frequency of persistent firing increases (decreases) with each pulse. During increase or decrease successive periods of silence were allowed to include runs where the persistent firing would only start after for example the second pulse. From the three recorded cells, we got successful runs of every protocol: DOWN-1 (two successful/3 total), DOWN-3 (2/2), DOWN-4 (6/8). UP-1 (7/15), UP-3 (1/1), UP-4 (4/7). An overview of all protocols is available at https://gin.g-node.org/ppfeiffer/cooperative_channels_in_biological_neurons_via_dynamic_clamp, together with the raw data and python scripts for analysis.

For the presented recordings, the model parameters are summarized in Table 1.

Data availability

The model for the cooperative ion channels and the neuron are described in the manuscript. Software for simulation of the cooperative channels is provided in git repositories linked in the manuscript (https://itbgit.biologie.hu-berlin.de/cooperativity/brutelib; copy archived athttps://github.com/elifesciences-publications/49974-brutelib; and https://itbgit.biologie.hu-berlin.de/cooperativity/dynamic_clamp_model; copy archived at https://github.com/elifesciences-publications/49974-dynamic_clamp_model). The recordings and analysis obtained in the dynamic clamp experiments are explained in the manuscript and are publicly available at https://gin.g-node.org/doi/cooperative_channels_in_biological_neurons_via_dynamic_clamp.

The following data sets were generated
    1. Pfeiffer P
    2. Egorov A
    3. Lorenz F
    4. Schleimer J
    5. Draguhn A
    6. Schreiber S
    (2019) G-Node
    Clusters of cooperative ion channels enable a membrane potential-based mechanism for short-term memory.
    https://doi.org/10.12751/g-node.08853e

References

  1. Book
    1. Allen LJS
    (2010)
    An Introduction to Stochastic Processes with Applications to Biology
    Chapman & Hall/CRC.
  2. Book
    1. Benda J
    (2002)
    Single Neuron Dynamics - Models Linking Theory and Experiment
    Humboldt-Universität zu Berlin.
  3. Book
    1. Thomas R
    2. D’Ari R
    (1990)
    Biological Feedback
    CRC Press.

Article and author information

Author details

  1. Paul Pfeiffer

    1. Institute for Theoretical Biology, Humboldt-Universität zu Berlin, Berlin, Germany
    2. Bernstein Center for Computational Neuroscience, Humboldt-Universität zu Berlin, Berlin, Germany
    Contribution
    Conceptualization, Data curation, Software, Formal analysis, Funding acquisition, Validation, Investigation, Visualization, Methodology
    For correspondence
    pfeiffpa@hu-berlin.de
    Competing interests
    No competing interests declared
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0001-5324-5886
  2. Alexei V Egorov

    Institute of Physiology and Pathophysiology, Heidelberg University, Heidelberg, Germany
    Contribution
    Supervision, Investigation, Methodology
    Competing interests
    No competing interests declared
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0003-4899-8407
  3. Franziska Lorenz

    Institute of Physiology and Pathophysiology, Heidelberg University, Heidelberg, Germany
    Contribution
    Data curation, Software, Investigation, Methodology
    Competing interests
    No competing interests declared
  4. Jan-Hendrik Schleimer

    1. Institute for Theoretical Biology, Humboldt-Universität zu Berlin, Berlin, Germany
    2. Bernstein Center for Computational Neuroscience, Humboldt-Universität zu Berlin, Berlin, Germany
    Contribution
    Conceptualization, Software, Formal analysis, Supervision, Validation, Methodology
    Competing interests
    No competing interests declared
  5. Andreas Draguhn

    Institute of Physiology and Pathophysiology, Heidelberg University, Heidelberg, Germany
    Contribution
    Resources, Supervision, Funding acquisition, Methodology, Project administration
    Competing interests
    No competing interests declared
  6. Susanne Schreiber

    1. Institute for Theoretical Biology, Humboldt-Universität zu Berlin, Berlin, Germany
    2. Bernstein Center for Computational Neuroscience, Humboldt-Universität zu Berlin, Berlin, Germany
    Contribution
    Conceptualization, Resources, Formal analysis, Supervision, Funding acquisition, Investigation, Visualization, Methodology, Project administration
    For correspondence
    s.schreiber@hu-berlin.de
    Competing interests
    No competing interests declared
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0003-3913-5650

Funding

Bundesministerium für Bildung und Forschung (01GQ0901)

  • Susanne Schreiber

Bundesministerium für Bildung und Forschung (01GQ1403)

  • Susanne Schreiber
  • Jan-Hendrik Schleimer

Deutsche Forschungsgemeinschaft (GRK 1589/2)

  • Paul Pfeiffer
  • Susanne Schreiber

Volkswagen Foundation (SmartStart2)

  • Paul Pfeiffer

Deutsche Forschungsgemeinschaft (SFB1134,AO1)

  • Alexei V Egorov
  • Franziska Lorenz
  • Andreas Draguhn

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 by BMBF (01GQ0901, 01GQ1403), DFG (GRK 1589/2), SFB1134, A01, and Volkswagen Foundation (SmartStart). We are grateful to Jan Benda for his dedicated support with dynamic clamp. We thank Fabian Santi and Dominik Spicher for fruitful contributions during an earlier stage of the project. Finally, we thank Nelson Niemeyer for valuable feedback on the final form of the manuscript.

We acknowledge support by the German Research Foundation (DFG) and the Open Access Publication Fund of Humboldt-Universität zu Berlin.

Ethics

Animal experimentation: All experimental protocols were conducted in compliance with German law and with the approval of the state government of Baden-Württemberg (Project T100/15).

Copyright

© 2020, Pfeiffer 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

  • 3,208
    views
  • 411
    downloads
  • 23
    citations

Views, downloads and citations are aggregated across all versions of this paper published by eLife.

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

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)

  1. Paul Pfeiffer
  2. Alexei V Egorov
  3. Franziska Lorenz
  4. Jan-Hendrik Schleimer
  5. Andreas Draguhn
  6. Susanne Schreiber
(2020)
Clusters of cooperative ion channels enable a membrane-potential-based mechanism for short-term memory
eLife 9:e49974.
https://doi.org/10.7554/eLife.49974

Share this article

https://doi.org/10.7554/eLife.49974

Further reading

    1. Neuroscience
    Shan Jiang, Lu Chen ... Chang-Rui Chen
    Research Article

    General anesthesia (GA) is an indispensable procedure necessary for safely and compassionately administering a significant number of surgical procedures and invasive diagnostic tests. However, the undesired stress response associated with GA causes delayed recovery and even increased morbidity in the clinic. Here, a core hypothalamic ensemble, corticotropin-releasing hormone neurons in the paraventricular nucleus of the hypothalamus (PVHCRH neurons), is discovered to play a role in regulating sevoflurane GA. Chemogenetic activation of these neurons delay the induction of and accelerated emergence from sevoflurane GA, whereas chemogenetic inhibition of PVHCRH neurons accelerates induction and delays awakening. Moreover, optogenetic stimulation of PVHCRH neurons induce rapid cortical activation during both the steady and deep sevoflurane GA state with burst-suppression oscillations. Interestingly, chemogenetic inhibition of PVHCRH neurons relieve the sevoflurane GA-elicited stress response (e.g., excessive self-grooming and elevated corticosterone level). These findings identify PVHCRH neurons modulate states of anesthesia in sevoflurane GA, being a part of anesthesia regulatory network of sevoflurane.

    1. Neuroscience
    Marie-Eve Picard, Miriam Kunz ... Pierre Rainville
    Short Report

    Pain is a private experience observable through various verbal and non-verbal behavioural manifestations, each of which may relate to different pain-related functions. Despite the importance of understanding the cerebral mechanisms underlying those manifestations, there is currently limited knowledge of the neural correlates of the facial expression of pain. In this functional magnetic resonance imaging (fMRI) study, noxious heat stimulation was applied in healthy volunteers and we tested if previously published brain signatures of pain were sensitive to pain expression. We then applied a multivariate pattern analysis to the fMRI data to predict the facial expression of pain. Results revealed the inability of previously developed pain neurosignatures to predict the facial expression of pain. We thus propose a facial expression of pain signature (FEPS) conveying distinctive information about the brain response to nociceptive stimulations with minimal or no overlap with other pain-relevant brain signatures associated with nociception, pain ratings, thermal pain aversiveness, or pain valuation. The FEPS may provide a distinctive functional characterization of the distributed cerebral response to nociceptive pain associated with the socio-communicative role of non-verbal pain expression. This underscores the complexity of pain phenomenology by reinforcing the view that neurosignatures conceived as biomarkers must be interpreted in relation to the specific pain manifestation(s) predicted and their underlying function(s). Future studies should explore other pain-relevant manifestations and assess the specificity of the FEPS against simulated pain expressions and other types of aversive or emotional states.