On the Firing Rate Dependency of the Phase Response Curve of Rat Purkinje Neurons In Vitro

The Phase Response Curve (PRC), a simple input-output characterization of single cells, can provide insights into individual and collective properties of neurons and networks, by quantifying the impact of an infinitesimal depolarizing current pulse on the time of occurrence of subsequent action potentials, while a neuron is firing tonically. Recently, the PRC theory applied to cerebellar Purkinje cells revealed that these behave as phase-independent integrators at low firing rates, and switch to a phase-dependent mode at high rates. Given the implications for computation and information processing in the cerebellum and the possible role of synchrony in the communication with its post-synaptic targets, we further explored the firing rate dependency of the PRC in Purkinje cells. We isolated key factors for the experimental estimation of the PRC and developed a closed-loop approach to reliably compute the PRC across diverse firing rates in the same cell. Our results show unambiguously that the PRC of individual Purkinje cells is firing rate dependent and that it smoothly transitions from phase independent integrator to a phase dependent mode. Using computational models we show that neither channel noise nor a realistic cell morphology are responsible for the rate dependent shift in the phase response curve.

In the particular case of neurons, PRCs can be employed to infer several network properties that are influenced by intrinsic membrane mechanisms. It has been shown that the PRC of tonically firing Purkinje Cells is flat at low firing rates, which has profound implications for information processing in the cerebellum. Here, we propose a novel method to estimate the PRC of single Purkinje cells at various firing rates and use it to unveil the smooth transition between flat and phasic PRC. Furthermore, we address potential explanations for the observed transition using computational modeling.

In addition, the extended range of PCs firing rates during behavior suggests that the rate of APs, its sudden transitions, its coherence across PCs, and the AP timing synchronization may contribute to information representation, processing, and downstream relaying. Thus, investigating how distinct firing regimes affect spontaneous and evoked response properties is imperative for dissecting cerebellar computation. Recently, key results from the mathematical theory of coupled oscillators sparked a lot of interest: a simple input-output characterization of the units composing a network, known as their phase response (or phase resetting) carve (PRC), is sufficient to classify and predict individual and collective properties. In the context of tonically firing neurons, the PRC quantifies the impact of an infinitesimal depolarizing current pulse on the time of occurrence of subsequent APs [5—10]. As the cell oscillates regularly, the pulse advances or delays the time of the next AP, depending on the oscillation phase (p corresponding to the time of pulse delivery. The resulting change of the time of the next AP can also be quantified in terms of the cell’s firing period and thus expressed as a phase shift Ag). By capturing the relationship between the evoked phase shift Ag) and the phase (p at which the input pulse occurred, the PRC predicts how, upon receiving weak synaptic inputs, neurons transiently delay or accelerate AP firing, contribute to network-wide AP synchrony, integrate external inputs or detect their temporal coincidences. So far, not only has the PRC been considered in theoretical and computational studies, but it has also been computed in experimental works (see [11] for a review), where different methods have been devised for its estimation [11—13]. Recently, Phoka et al. 2010 [14] proposed a correction to a traditional estimation method and tested it in PCs of juvenile mice. Unexpectedly, they reported that the PC’s intrinsic firing rate has a profound effect on the response properties: the PRC of PCs firing at low rates displays a flat profile, suggesting that neurons behave like phase-independent inputs integrators; on the other hand, the PRC of PCs firing at high firing rates has a prominent peak, indicating a phase preference similar to coincidence detectors. While it was not the first time that PRCs were shown to undergo changes over a range of AP frequencies [15] , the wide physiological range of PCs spontaneous firing rates and their ease of experimental access in in vitro preparations, made the report on the rate dependence of the PRC relevant. Furthermore, intrinsic membrane properties might promote synchrony in a way that is relevant to information processing [16] , particularly in the cerebellum [17, 18].

We aimed at its systematic exploration, both at the single-cell and at the population levels, which may be directly relevant for modeling studies. In particular, in the light of the known bistable behavior of PCs, and their ability to abruptly toggle between distinct AP firing rates, we found it urgent to clarify whether the changes in PRC occur abruptly or smoothly, for increasing AP frequencies. In addition, we tested the effect of the current pulse amplitude, verifying that the PRC prominent phase-depen-dency of PCs firing at high rates is an intrinsic property and not an artifact of the stimulation protocol. The key contribution of this work is twofold: (i) we developed a novel ad hoc closed-loop electrophysiological protocol to regulate PCs slow scale adaptation and achieve highly significant PRC estimates at fixed firing rates in a relatively short experimental time. By such an approach, (ii) we confirmed and extended the observations of [14], considerably improving the earlier observation statistics, and demonstrating conclusively and unambiguously that no abrupt sWitch in PRC occurs. Instead, PCs smoothly shift from integrators to coincidence detectors, as their AP frequency increases. Finally, verifying that these observations are not affected by the particular PRC estimation method, we tested our conclusions employing both the corrected direct method, as in [14], and an indirect method [11—13].

Ethics statement All procedures were performed according to institutional and national ethical guidelines (license no. LA1100469 from the Belgian Federal Public Service Health, Food Chain Safety and Environment).

Briefly, after isolating the cerebellar vermis, the tissue was glued with cyanoacrylate glue to a flat metal platform and surrounded by agar blocks to improve stability during slicing [3, 20]; the tissue was then cut in 250 um slices using a Vibratome (VTlOOOS, Leica Microsystems, Wetzlar, Germany) in ice-cold artificial cerebrospinal fluid (ACSF), containing (in mM): 125 NaCl, 2.5 KCl, 1.25 NaHZPO4, 26 NaHCO3, 25 glucose, 2 CaClz, and 1 MgClz, balanced with 95% 02 and 5% C02. The slices were incubated for 30 — 45 min at 32°C and then stored at room temperature, until they were transferred to the recording chamber of a fixed-stage upright microscope (DMLFS, Leica Microsystems, Wetzlar, Germany). The microscope was equipped with differential interference contrast (DIC) Video-microscopy and mounted a 63X water immersion objective.

Some PCs were filled with Lucifer yellow and imaged by epifluorescence microscopy, confirming that the entire dendrite was always in the plane of the slice. Whole-cell patch-clamp recordings were performed at 33 i 1°C, employing an EPC10 amplifier (HEKA, Lambrecht/Pfalz, Germany) or an Axon Multiclamp 700B amplifier (Molecular Devices, USA), both used in current-clamp mode. Patch electrodes were pulled from thick-walled borosilicate glass capillaries (1BF150, World Precision Instruments, Hitchin, UK) with a horizontal puller (P97, Sutter, Novato, USA) to a resistance of 3 — 6 MO. Electrodes were filled with an intracellular solution containing (in mM): 130 methanesulfonic acid, 10 HEPES, 7 KCl, 0.05 EGTA, 2 NaZATP, 2 MgATP, 0.5 NaZGTP, and pH adjusted to 7.3 with KOH. All recordings were obtained employing ACSF as the extracellular solution, balanced with 95% 02 and 5% C02, and routinely supplemented with 10 uM SR95331 (a selective antagonist of GABAA receptors) to abolish incoming spontaneous synaptic potentials. Amplified analog signals were low-pass filtered at 10 kHz, sampled at a rate of 30 kHz and digitized at 16 bits with a DAQ board (PC1-6229, National Instruments, USA). The same board was used to generate the amplifier control commands waveforms, synthesized at the same rate and resolution of the data acquisition. Stimulation and response data were generated and collected by using the public domain software LCG [21] , and analyzed by custom scripts written in MATLAB (The Mathworks, Natick, MA). The amplifier builtin online capacitance compensation circuitry was always applied, while the online bridge balancing circuitry was employed alternatively to the offline (software) active electrode compensation (ABC) [22], whose implementation is builtin in LCG [21]. Liquid junction potentials were left uncorrected and all the chemicals and drugs were obtained from Sigma-Aldrich (Diegem, Belgium). Analysis scripts, LCG configuration files and LCG command line strings to precisely replicate our eXperimental protocol and modeling are available from ModelDB [23] at http://senselab.med. Phase response curve estimates

Applying direct methods [24] in tonically firing cells, such as the PCs, required the repeated injection of very brief square pulses of current (i.e., Ipulse = 50 — 150 pA, Tpulse = 0.5 — 1 ms, at least 1400 repetitions), each timed at a different phases (p of the cell f1ring cycle (e.g., Fig. 1, panels E and A). The phase-shift Ag) of the next AP induced by each pulse (Fig. 1, panels E, B and C), was first quantified and then normalized by the total injected charge Q = Ipulse - Tpulse, [8, 25] allowing comparison across stimulation conditions. Briefly, upon (online or offline, see below) digital detection of the timing tk of individual AP peaks, the mean (181) of the distribution of inter-spike intervals 181k 2 (tk+1 — tk) was computed, and taken as an estimate of the (regular) f1ring period. The occurrence of each pulse was expressed as the corresponding phase (p = T/ (181), by relating its absolute time of occurrence tpulse to the AP immediately before (i.e., say tj),

Note that due to jitter in the next AP, the value of (p may slightly exceed its upper theoretical limit (p = 1 (i.e., (p E [0; 1 + 8] ). Because with no pulse the next AP would have occurred at tj + (181), the actual phase-shift induced by the external perturbation was determined as Ag) 2 ((ISI) — ISIPermbedV (181), with ISIperturbed = tj+1 — tj and where tj+1 is the (perturbed) time of the AP immediately following the pulse. By this convention, positive (negative) values of Ag) represent phase advances (delays). Finally, normalizing Ag) to the charge Q of each pulse, the traditional direct estimate of the PRC can be expressed as: estimator

Experimental setup. (A) A real-time closed-loop in vitro setup for experimentally estimating the PRC at a fixed firing rate, based on a Proportional-lntegral-Derivative (PID) controller and on a reactive-clamp paradigm, was employed to rapidly and optimally explore the dependency of the PRC of PCs on the cell’s firing rate. To this aim, a brief external current-pulse (B-D at different scales) was repeatedly delivered at different phases (p of the cell cycle (E), and the resulting impact on the time of the next AP was quantified as a

Despite its potential bias due to the non-station-ary firing regimes, we used it here solely as a control method and for confirming the firing rate dependency of the PRC in PCs. In tonically firing cells, such as the PCs, WSTA required repeatedly recording the times {t}k of APs elicited by weak-amplitude fluctuating currents I (1‘), generated as exponentially filtered white-noise [26] f(t), lasting for ~ 30 s:

Instead of binning and constructing a histogram of the sampled PRC data-points, a standard nonparametric smoothing technique based on Gaussian kernel convolution [28] was employed as in [14], aimed at increasing the signal-to-noise ratio of each PRC estimate: In discrete coordinates, the convolution integral became a sum over each of the N data points available, With the kernel K(x) centered over each available phase. The optimal kernel Throughout the text and in the figures, the smoothed PRC estimates have been indicated for simplicity as Z((p).

Due to the intrinsic variability of each cell’s inter-spike intervals, repeating the pulse injections over and over in time conveniently allowed us to sample uniformly the range of (p, while stimulating the cell at a frequency much lower than its firing rate (i.e., 2—6 pulse/sec). As a consequence of the need to study the firing rate dependency, large parts of the recording were often discarded when the cell was not firing at a fixed rate, further increasing the time needed to obtain a PRC estimate (often greater than 30 min for a single firing rate).

In a first series of experiments, a constant holding current was therefore applied (i.e., on the top of Ipulse or of I (1‘)), adapting its value manually from —0.2 to 1 nA, to alter the firing rate of PCs by depolarizing or hyperpolarizing their membranes. Depolarizing or hyperpolarizing a cell instantaneously alters its firing rate, although several minutes are typically required for the cell to reach a (new) steady firing rate. This caused long waiting intervals before estimating the PRC at a given firing rate. In addition, occasional slow drifts of the mean inter-spike intervals occurred, over a window of several seconds, thus altering or biasing the shape of the PRC. To address these limitations, in an additional set of experiments, we made use of a spike rate controller [29], using a closed-loop paradigm similar to the one employed in [30] and inspired by the response-clamp paradigm [31] was adopted: the frequency-clamp. In short, an iterative estimate 13” of the cell’s instantaneous firing rate was updated online after each AP, detected in real-time (i.e., as a positive crossing of a voltage threshold), using the following formula: where T = 1s acts as the time scale over which the instantaneous firing rate is estimated, weighing each new AP and the previous firing history [21, 31]. The running value of 13k was com-— 13k. This was fed into a Proportional-Integral-Derivative controller (PID), realized via software in LCG, and employed to automatically update in closed-loop the value of the constant holding current pared to a target frequency Ptarget and employed to define an error signal ek 2 Raw where gp, g1, gD are the proportional, integral, and derivative gains, respectively (i.e., gp = 0.001 pA/HZ, gl 2 0.1 pA/Hz, gD = 0 pA/HZ). The value of Ptarget was also used as the initial value for the estimator P0, in order to reduce undesired transients. While the output of the PID controller was updated every time an AP was detected, it was held constant during the time interval starting from the AP preceding the perturbation pulse to the second spike following it, in contrast to what was done in [30]. In other words, the PID controller was temporarily disconnected, while holding its most recent output, before delivering the external current perturbation required for estimating the PRC, in order to minimize any artifact. A detailed description of the experimental setup is given elsewhere [21].

In addition, as in closed-loop the timing of the external current pulse required for the PRC estimation could be precisely chosen in reaction to an AP and after a certain time delay, we optimally synthesized the values of these delays in order to sample the range [0; 1] of (p with maximal efficiency (i.e., more uniformly than pseudo-random number generation). We employed a Sobol sequence, first used in [32] and described in [33]: briefly, after online detection of an AP, the external pulse was delivered in a reactive-clamp fashion [34] after a delay Ti, generated as the i-th element of a Sobol sequence (Grey code variant [35] , after discarding the initial points). This was repeated at least 1400 times at a rate of one perturbation every 6 APs, independently of the cell’s firing rate.

Furthermore the PRCs estimated are in good agreement with those obtained for PCs using a similar method [14] , which provides an additional control to our methods. In order to study the influence of the PID controller in our PRC estimates we used the Khaliq and Raman [36] model (see SI Fig, Discussion and Modeling methods).

In order to obtain a concise description of a PRC and to compare those obtained at different AP firing rates, for the same cell and across cells, we adopted the peak-to-baseline ratio r as in [14]: Where ml and me are the values of the late and of the early local extrema (i.e., largest peaks in absolute values) for each of the two halves of Z((p) (i.e., in (p E [0; 0.5] and in (p E [O.5;1]). Note that When ml and me have opposite signs (e.g., as in type II PRCs, [37]), r = 1. This also allowed us to concisely quantify the dependency of the PRC shape on the firing rate F, by fitting to r(F) the parameters of a sigmoidal function

All simulations were performed using the NEURON simulation environment [38]. The simulation code is available on ModelDB or by request from the corresponding author.

In a first set of simulations, we used the single-compartment, conductance-based PC model described in [36] and available on ModelDB at the URL https:// senselab.med.yale.edu/modeldb, accession number 48332. To test the hypothesis that voltage fluctuations endogenously generated by the random opening and closing of ion channels might influence the shape of the PRC, we incorporated channel noise into the model using the method described in [39]. Briefly, the fluctuations induced by channel noise can be accounted for by extending the dynamics of the ionic conductances present in the model according to the equation Where g is the maximal conductance, p0(t) is the fraction of open channels (for a deterministic model, g(t) = g 190(0), N is the number of states of the equivalent kinetic scheme and each 17,-(t) is the solution to a stochastic differential equation of the form where 5,-(1‘) is a delta-correlated Gaussian process with zero mean and unitary variance. For the general case of arbitrary kinetic schemes, the N — 1 time constants T,- and standard deviations a,are obtained numerically from the N x N transition matrix of the system that contains the transition rates between all possible states in the kinetic scheme. In the case of the Khaliq-Raman model, this approach was employed only for the resurgent sodium current, which is described by a kinetic scheme that cannot be mapped into the composition of multiple two-state subunits. For all other ionic conductances, the coefficients T,- and a,- were analytically calculated using the procedure detailed in [39].

We chose two values of CV, low (around 5%) and high (around 10%), corresponding to values of both length and diameter of 160 and 80 pm, respectively. In the deterministic model, the length and diameter were set to 80 um and an additional noisy current, modeled as delta-correlated Gaussian white noise, was injected to obtain comparable values of CV.

A constant current offset was injected to vary the baseline firing frequency of the model. In some simulations, this offset current was computed by a PID controller, to replicate and validate in silico the closed-loop technique employed in the experiments. Multi—compartment model. To elucidate whether the presence of an extensive dendritic tree might influence the firing rate dependency of the PRC, we used the De Schutter-Bower model [40, 41].

In both cases, the algorithm employed to compute the PRC differed from that used in the single compartment model and resembled the one adopted in [42, 43]: briefly, after evolving the model until it reached a steady state, we identified two spike times to and t1 such that the 181 1‘1 — to was of the appropriate duration. Then, we simulated the model again until to — 5 ms and saved the full state of the model at this point in time. Finally, starting from to — 5 ms, we evolved the model until t1 + 10 ms for N = 50 trials: current perturbations (0.5 ms duration and 0.2 nA amplitude) were applied at times given by where t; is the perturbation time in the i-th trial. In the case of the model Without synaptic inputs we used the PMlO model [40], set the temperature of the simulation to 28°C and injected a somatic current of varying amplitude to span a Wide range of firing rates. In the case of the model With synaptic inputs, we used the PM9 model With synapses distributed on the dendritic tree as described in [41], set the temperature to 37°C and fixed the presynaptic excitatory and inhibitory firing rates to 35 and 2 HZ, respectively.

The passive membrane properties of the cells were measured in terms of the input resistance (21.8 i 4.7 M52) and membrane time constant (40.7 i 19.8 ms). We employed the current-clamp configuration and focused on the input-output relationship between the phase shift Ag), induced in the cell firing cycle by an external current pulse, and the phase (p at which the pulse was repeatedly delivered (Fig. 1; see Materials and Methods). Known as the phase response curve (PRC), this characterization has been increasingly employed to study and classify cellular excitability [11], accompanying conventional descriptions such as, e.g., the frequency-current curve. Previous studies suggested that at least 500 repeated stimulations (i.e., trials) are required to accurately estimate PRCs in tonically firing neurons using direct estimation methods [24]. Here we used 946 to 17278 trials (mean = 2650) to compute 112 PRCs from 42 PCs using the direct method. While firing spontaneously and regularly in vitro over long periods, and responding to DC depolarizing or hyperpolarizing holding currents by increasing or decreasing their firing rate, PCs eXperience substantial firing rate adaptation over very long time scales. Methods to estimate the PRC under transient conditions have been recently developed: [44] however, in order to facilitate the comparison of our results to those of [14] and because the firing rate range that we want to probe is large, we choose to acquire PRCs in a regime of stationary firing. This often leads to discarding a series of stimulation trials in order to avoid an artificially skewed distribution of inter-spike intervals. Similarly, it is often necessary to wait until PCs reach a steady state firing rate before initiating the repeated stimulation protocol. This makes computing several PRCs in the same cell impractical, for instance when one is interested in investigating the effects attributable to the cell’s firing regime, pharmacological manipulation, or recording conditions. We therefore designed and implemented via software a real-time closed-loop control system (Fig. 1A), employed to speedup convergence to the firing rate steady state and reduce very slow fluctuations. In addition, following [45] we employed a quasi-random

Using the Khaliq-Raman model (see Methods), we investigated how the closed loop system affects the PRC estimate. Phase response curves estimated with open or closed-loop methods were strikingly similar (81 Fig). This system allowed us to both maintain the cell at the desired target frequency and significantly reduce the time required for a reliable estimation of the PRC. We note that the approach of [44] could equally be used to study the firing rate dependency of the PRCs.

Fig. 2A displays 12 PRCs, estimated under stable recording conditions, while repeatedly altering the PC firing rate in a shuffled order. While at low firing rates the profile of the PRC appears relatively flat and independent off the input phase (p, at higher rates the profile changes: a late-phase peak (i.e., in the range [0.5; 1]) becomes sharper and shifts to the left, while the average amplitudes in [0; 0.5] decrease. This was quantified for this PC, and four other individual cells, by defining the peak-to-baseline ratio (Fig. 2B; see Methods). This is a measure of the absolute difference between absolute peak values in the ranges [0; 0.5] and [0.5; 1], and it is maximally 1 when these two peaks have different signs. For all individual neurons, the peak-to-baseline ratio increased smoothly in the range of physiological firing rates under consideration, and it could be best fit by a sigmoidal function.

Since we observed a rate-dependent shift in the PRC late-peak (Figs. 2A and 3A) we asked whether this reflects a time-to-spike preference. To test this possibility we applied the following change of variables (p = 1 + tAp/ (181), where the tAp is the relative time to the AP following the stimulus, and plotted the PRC as a function of time, for distinct firing rates (Fig. 3D). The location of the maxima of these plots (Figs. insets in 2B and 3B), displayed a marked dependence on the firing rate in the range [0; 100] Hz while, at very high firing rates, i.e. above 100 HZ, it became rate-independent and equal to 2 ms. This suggests that the rate-dependence of the PRC late-peak (Figs. 2A and 3A) does not result from a rate-independent time-to-spike preference. Our observations are not changed when using the truncated Gaussian method [45] as illustrated in S2 Fig. In an initial set of experiments, we wondered whether the rather flat profile of the PRC observed at low firing rates was an intrinsic property of the cells. Very weak input stimuli may be in fact effectively ignored by the cell and result in a phase-independent PRC profile. The benefits of systematically acquiring PRCs, across distinct conditions along stable recording sessions, were again exploited: we injected in the same PCs an external input with different amplitudes, over distinct firing rates. An example of such experiments is reported in panels AB of Fig. 4, where the impact of the pulse amplitude is apparent: the signal-to-noise ratio increases for stronger pulses, while the PRC is indeed phase-independent. Fig. 4C-D further visualize graphically the 68% confidence intervals of the PRC estimates, revealing an almost two fold reduction when using doubled pulse amplitudes.

When expressed in terms of a population summary, the PRCs confirmed our previous observation (Fig. 3A). When pooling the PRCs obtained across cells, in

Direct method for the estimation of the PRC. (A) A direct estimate of the PRC, obtained for the same PC, is plotted after data smoothing (see Methods) while altering the cell’s firing rate in the range 20 —160 Hz): a strong dependency on the firing rate is apparent. The transition from an approximately flat to a phase-dependent PRC profile does not occur abruptly, but smoothly: in each subplot, the horizontal gray dashed lines represent Z(cp) = 0, while the continuous black thick traces are the PRCs, estimated at distinct firing rates (Fig. 1). Black circles indicate the location of the extrema for each of the two halves of the curves (i.e., in cp e [0; 0.5] or [0.5; 1], emphasized by the vertical thin black line), used to concisely characterize the PRC shape according to its peak-to-baseline ratio (see Methods). The graded PRC shape dependency on the firing rate is confirmed in three other PCs (B, markers) and quantified by their peak-to-baseline ratio. The black cun/e represents the function (1 + e‘(F‘a)/b)‘1, with best-fit parameters a = 44.1 , b = 20.5. The inset further displays the location of the PRC peak, relative to the time of the AP following the stimulus (i.e., rpeak = (rpeak — 1) - (ISI)), for the same five cells. 127 Hz 82 Hz 112 Hz 97 Hz Peak to baseline ratio 0 io

Population summary using the direct method. (A) Population summary obtained over distinct firing rates, by averaging PRCs across individual cells. (B) Individual responses are plotted in gray and pooled according to the corresponding firing rate. Cell numbers are further indicated in parentheses. This summary, quantified by the peak-to-baseline ratio as in Fig. 20 (individual cells, n = 42: gray markers; averages from A: red markers), confirms (C) our obsen/ations in single PCs (Fig. 2). The black cun/e represents the function 15 HZ-wide bins, according to the PC firing rate at which each curve was measured, a moderate amount of variability was observed (Fig. 3B). The overall quantification in terms of the peak-to-baseline ratio, already discussed in Fig. 2B, is comparable to that measured for single cells, revealing and confirming the same marked smooth dependency on the firing rate (Fig. 3C). As in Fig. 2B, the inset of Fig. 3C displays the location in time of the peak of the PRC against the firing rate, and the average PRCs were also plotted as a function of time to the next AP

3D). While slightly noisier than the data acquired Within the same PCs, this evidence prompts us to exclude that the late-peak in the PRC is an artifact of the phase normalization.

We computed 95 PRCs in 16 PCs, with firing rates below 100 HZ. The average PRC profiles, obtained from pooled data as in Fig. 3B across several firing rates, are shown in Fig. 5A-B were the PRCs are plotted both as a function of the phase of the time to the next AP. As for Figs. 2 and 3, performing a quantification based on the peak-to-baseline ratio revealed a qualitatively similar frequency dependence of the PRC profile on the PC firing rate.

To some extent, we attribute this inaccuracy to the WSTA method: while for direct methods the stationarity of the firing rate could be precisely monitored and controlled in closed-loop, indirect methods do not allow the same precision, as they require the injection of a noisy current waveform that elicits a train of APs with some variability. It should be noted that a method has been developed recently that allows a more accurate estimation of PRCs while using fewer spikes [44, 46]. We did not employ this method since we were mostly concerned with validating the results obtained with the direct method. Nonetheless, the presence of a late peak in the PRCs at high firing rates confirms that our observations do not depend on the PRC estimation method in use. When referred to the time of the next AP, the location of the PRC peak averaged to a value of —1.7 i 0.3 ms (Fig. 5C, inset), with the exception of some (n = 4) cases where the PC fired at low rates and the PRC did not exhibit a peak at late phases.

Active conductances are thought to modulate the shape of the PRC, and therefore computational modeling [7, 47, 48] could be a powerful tool to dissect the ionic bases of the PRC. The 0.6 Peak time conductance-based models developed in [36, 49] to recapitulate several experimental observations in Purkinje cells were used in [14] to attempt to reproduce the PRC shape and in particular its dependency on the firing rate. Such attempts have proved largely unsuccessful and up to today a model able to reproduce the rate dependency observed experimentally remains elusive. Having observed a rate dependency in the CV of PCs in our experimental data (81 Fig, panel C) we asked Whether spiking variability could account for the flat PRC profile observed at low firing rates. To test this hypothesis, we employed the Khaliq-Raman model [36] and computed PRCs at low and high firing rates. We compared the results obtained When spiking variability was introduced via additive noise fluctuations (i.e., by injecting a noisy current into the model

Stochastic channel noise is not responsible for the firing rate dependency of the PRC. The PRCs of the model incorporating channel noise (gray curves) are strikingly similar to those of the deterministic model (black curves). The variability of the spike trains does not affect the curves either at low (A and C) or high (B and D) firing rates (left and right columns, respectively). The membrane area of the model with channel noise was chosen to obtain the desired value of CV at low firing rate (around 10% for A, and 5% for C) and that same area was used for the simulations at high firing rate (panels B and D). Comparable values of CV were obtained in the deterministic model by changing the magnitude of the fluctuations of the injected current. For a given perturbation size, the model is less sensitive at high firing rates (B and D).

The results are shown in Fig. 6: solid lines denote the stochastic model, Whereas dashed lines represent the deterministic model With external noisy current. In both conditions, the model is more sensitive to perturbations When firing at low rates (left panels (A and C)), Which translates both to a larger amplitude of the PRC and to a higher CV, for the same number of channels (and therefore magnitude of the internally generated fluctuations) or variance of the noisy current. This result holds true both in the case of low and high variability, as exemplified by the top and bottom panels, respectively. Importantly, in all four conditions tested (low and high firing rate, low and high CV), the PRCs of the stochastic and deterministic models are strikingly similar and, for a given firing rate, do not depend on the amount of endogenous or exogenous noise. Furthermore, the models fail at reproducing the flat PRCs observed in vitro at low firing rates and the greater sensitivity to perturbation amplitude at high firing rates typical of the eXperiments.

We compute PRCs spanning firing rates in the range 33 to 77 HZ (Fig. S3 Fig) and observed a decrease in Z((p) for perturbations delivered in the early phase of the spiking cycle. Additionally, in our simulations early-phase Z((p) became negative for firing rates above approximately 50 HZ. The shape at late phases remained unchanged even for high firing rates, Which is in contrast With the experiments (compare S3 Fig With Fig. 2), in Which the late peak drifts towards the middle of the

The peak-to-baseline ratio of the multi-compartmental model seems to capture some of the experimentally observed features of the PRC, although interpretation of these results requires great care, see the Discussion. Purkinje cells are endowed with very large dendritic trees with thousands of contacts from parallel fibers, which might provide a dendritic load—even in the relatively quiet conditions of the slice preparation—by means of spontaneous synaptic release. To test whether distributed synaptic inputs might influence the shape of the PRC, we delivered random synaptic inputs to the dendrites of the model neuron. We found that, even when the instantaneous rates of inhibitory and excitatory inputs were relatively stable, synaptic inputs greatly affected the shape of the PRC (S3 Fig, panel A). These simulations suggest that a large portion of the variability observed in the PRC shape in Purkinje cells, even for the same firing rate (Fig. 3B), could be attributable to spontaneous activation of dendritic conductances, given that even similar presynaptic patterns can result in very different PRCs (S3 Fig, panels B and C). We conclude that realistic dendritic activation interferes with the PRC measurement, and, even though for some selected trials (see Methods) the model seems able to partly reproduce the experimental observations, the PRCs computed in these conditions are extremely variable.

However, while their main focus was to introduce and test a novel technique for the PRC estimation, their experimental data set consisted of a small population of 16 cells. In addition, only for a very small number of PCs, could the response of the same cell be investigated under stable recording conditions, while studying two firing regimes (i.e., low and high firing rates). In fact, despite the increased accuracy, the PRC estimation technique they introduced requires a high number of stimulation trials for each recorded cell: for instance, at least 10000 trials were used in Fig. 5 of [14], which correspond to ~ 1 h when sampling low firing rates. Aiming at replicating their intriguing discovery, we immediately realized that this made rather impractical a more systematic exploration of a wider set of firing rates in the same cell: stable recording conditions could not be maintained except under extraordinary circumstances. Our efforts confirm and extend the previous report, aiming at a consistent, accurate, and conclusive evaluation of the PRC in PCs. We specifically implemented a real-time system composed of a closed-loop controller of the cell instantaneous firing rate. This allowed us to compensate for the intrinsic variability in the firing rate, the extremely long transients before steady-state, and the very slow drifts. Ultimately, we could obtain accurate PRC estimates with fewer trials and over reduced intervals of time and thus explore a broad range of firing frequencies (i.e., 20 — 150 HZ). Of course, the controller was disabled each time a stimulation was applied, to decrease the chance of introducing artifacts. Our data clearly and unambiguously demonstrate a strong modulation in the shape of the PRC by the firing rate, in all PCs we recorded and irrespective of the estimation method and the amplitude of the external perturbation pulse (Figs. 2—5). In particular, PCs smoothly alter their dynamical response properties, from non-ideal integrators at low firing rate (i.e., ~ 20 HZ) to phase-dependent integrators at higher firing rates (i.e., ~ 80 HZ) (Figs. 2, 3 and 5 and also SS Fig).

In [50] and [8], less than 180 trials were used to compute the PRC of individual cortical neurons by direct methods, although current pulses were stronger in amplitude and longer in duration than those used in PCs in [14]. In addition, cortical neurons fire APs only upon external holding current injection. In comparison, PRCs of PCs are generally very noisy and require a substantially higher number of trials, however employing larger pulse amplitude reduces PRC raw data variability and increases confidence of its estimate (Fig. 4). Along the lines of improving the confidence of the PRC estimates, novel (e.g., indirect) methods and closed-loop eXperimental protocol recently received increasing attention [13, 29—32, 51]. Our work has been greatly inspired by those previous proposals as we combined some of the existing paradigms together, benefiting from advantages of the (corrected) direct method. This was an important design element in our experiments, enabling to study PCs across a variety of firing rates and ultimately helped to unambiguously conclude on a smooth transition of the PRC shape, and improved the firing rate resolution. While PRCs estimated by the indirect method did not perfectly match those obtained by the (corrected) direct method, smooth firing rate dependency of PRCs was nonetheless confirmed.

To clarify this point, we compared open with closed-loop estimates in a single-compartment conductance-based model of a Purkinje cell [36] (81 Fig, panel A). The usage of a PID controller resulted in a slight increase in the coefficient of variation of the unperturbed trials (approximately 0.07% increase) while leaving the curves almost unchanged. We observed no increase in the variability of inter-spike intervals in the experiments in which the PID controller was employed (filled markers in 81 Fig, panel B).

For PCs, mathematical modeling might help identifying candidate ionic bases for the observed firing rate dependency of PRCs. In [14], several attempts to reproduce the firing rate dependency of the PRC were described, in particular employing the conduc-tance-based models proposed by Khaliq and Raman [36] and by Akemann and Knopfel [49]: none reproduced the phase-independent PRC profile at low firing rates and the overall firing rate modulation. In our experimental data set, the inter-spike intervals distributions displayed a significant inverse correlation of the coefficient of variation (CV) with the firing rate (i.e., Pearson’s r = —0.4, p < 10—6, and a slope of —0.25/kHz—Sl Fig, panel B) regardless of whether open or closed-loop methods used. We first simulated the deterministic KR model, in which we injected a noisy current given

The value of IS was chosen in order to obtain the desired variability of the 181s (corresponding to a CV between 5 and 10%). The results are shown in Fig. 6 (black traces), for high and low firing rates (left and right columns, respectively). Note how the amount of variability in the spiking response primarily influences the dispersion of the points in the PRC but not its shape at either high or low firing rates. We then speculated that stochastic fluctuations of voltage-gated channels could play a role in disrupting phase-preference at low firing rates, given the (inverse) rate-dependent AP jitter variance. An augmented version of the Khaliq-Raman PC model, with stochastic channel gating [39] and geometry adjusted to reproduce realistic inter-spike interval distributions was used to investigate whether stochastic fluctuations of ionic channels could be responsible for the observed phenomena. This was of particular interest since the Khaliq-Raman model describes also potassium channels with large single-channel conductance (i.e., > 200 nS) [36] where channel noise could play a larger role. We observed no difference between the stochastic and the deterministic models, even though the CV was reduced at high firing rates in the stochastic model (for the same cell size and channel density). Nonetheless, biophysically accurate (i.e., voltage-dependent) channel flickering had a similar effect on the PRC as what was obtained by additive Gaussian noise current injection, in [14] and our simulations.

Given the known relationship between excitability type and PRC profile [5] , it is surprising that the model displays a fully positive PRC with a peak around (p = 0.5 at low firing rates, which then turns into a slightly more bimodal profile as the firing rate is increased, in a way that is reminiscent of a type II PRC, with a peak around (p = 0.7 and a slightly negative component at early phases—see Fig. 6 and S3 Fig of [14] (the model is available from ModelDB [23] at http://senselab.med.yale.edu/modeldb (accession number 155735).

This has to our knowledge not been previously highlighted and suggests that the Khaliq-Raman (as well as the Akemann and Knopfel) model may provide grounds for further investigations on the modulation of the PRC by the firing rate. Another overlooked point in previous studies is that as the firing rate is increased the sensitivity of the neuron to the perturbation decreases dramatically (Fig. 6—low versus high firing rate). This is in disagreement with what was observed in the experiments and illustrates how a biophysically realistic model can fail to describe simple neuronal behaviors while taking in consideration complex channel kinetics.

Having shown that the shape of the PRC is not affected by the CV of the 181, we first simulated the model with the injection of a DC current (i.e., Is = O in Eq. 14, leading to a CV = 0) of varying amplitude to span several firing rates. Additionally, we set the simulation temperature to 28°C, in order to span a broader range of steady state firing frequencies, particularly in the low end of the spectrum. While the model failed to reproduce the shape of experimental PRCs at low firing rates, it produced curves characterized by a peak at early phases that becomes negative as the firing rate is increased (S3 Fig, panel A). This behavior results in lower peak to baseline values at low firing rates (S3 Fig, panel B) even though the shape of the curves does not resemble that of the experimental PRCs. Finally, similarly to the KR model, but in contrast with the experiments, the model is more sensitive to perturbations at low firing rates.

However this study’s focus is on whether the flat PRC at low rates observed in vitro could be explained by the presence of the dendritic tree, and not in the effects of complex dendritic morphologies on somatic PRCs. It is then valid to conclude that the PRCs firing rate dependency in PCs can not be explained by a complex model with realistic morphology.

Furthermore, we wanted to conduct simulations with the correct temperature factor since temperature severely affects channel dynamics. In order to simultaneously add dendritic load and increase the temperature of the model to 37°C, we introduced synaptic conductances as described in [41], see Methods. Under these conditions, the model fires irregularly so we focused in computing PRCs for selected inter-spike intervals, for which the expected 181 was known (see Methods). Despite our focus on inter-spike intervals that did not have abrupt conductance fluctuations (S3 Fig, panel B), even a slight fluctuation in the rate of dendritic events severely affected the shape of the somatic PRCs (S3 Fig, panel C). Nevertheless for selected inter-spike intervals, the model with dendritic load could approximate the strong difference in PRC shape at low versus higher rates (S3 Fig, panel A). However, given the high sensitivity to dendritic conductances, care should be taken when interpreting this curves as discussed above. Finally, for questions related to network dynamics, network simulations using the PRCs measured in this study could result in more accurate predictions than simulations with conductance based models that do not accurately reproduce the rate dependency of the PRC in the PC.

Considering a framework of elementary computation, we discuss the most immediate consequence of a firing rate-dependent PRC.

The cell would then favour a rate code at low firing rates and a time based code at high firing rates which is a profound implication for information processing in the cerebellum. It has been suggested that neuronal coding depends both on the operating mode and the properties of the stimulus [52] , therefore the fact that intracellular properties favor different coding schemes is relevant only if the inputs to the neurons are adjuvant to this hypothesis. It turns out that, a recent study that looked at the relation between population coding by PCs and saccadic eye movements found that the same PC can exhibit multiplexed spike coding, suggesting that despite the reduced level of overall correlation between PC firing rates a small fraction of the spikes can couple to the network activity and report different behavioral information [53]. Our study suggests that the intrinsic properties of PCs favor multiplexed information coding and that there is a smooth dependency on the firing rate of how sensitive PCs are to synchronous inputs.

Furthermore, in the cerebellar microcircuitry, it was shown that the interplay between the activation of parallel fibers and of interneurons affects PCs firing, and was associated with a variety of motor and sensory tasks [54, 55]. Moreover, the incoming information is likely to be processed and conveyed downstream by means of the synchronization of several PCs [18] and high frequency synchrony has been reported in the Purkinje cell layer [17] and has been associated with PC collaterals. We provide a thorough characterization of single cell properties that might be involved in synchronization, hence contributing for the observed phenomena. It is in this perspective that the study of PRC of PCs is of great relevance, as it provides insights in the synchronization mechanisms of neuronal populations. The fact that Purkinje cell’s phase response curve smoothly transitions from phase independent at low firing rates to highly bimodal at high firing rates could then be used as an intrinsic mechanism to modulate the contribution of individual inputs to spiking synchrony across Purkinje cells.

(A) The Khaliq-Raman model neuron was simulated With or Without a closed-loop (PID) regulation of its firing rate, as in the experiments. The model PRC was estimated as in Fig. 3 and studied in four conditions: lOW or high firing rate and With or Without PID. The PRC estimates using the PID are strikingly similar to those estimated using the open-loop method. (B) The CV of the eXperiments With PID (filled markers) are comparable to that of the experiments Without PID (open 82 Fig. PRCs estimated with the truncated Gaussian method. Same as Fig. 3 for the truncated Gaussian method. (A) Population summary obtained over distinct firing rates, by averaging PRCs across individual cells. (B) Individual responses are plotted in gray and pooled according to the corresponding firing rate. Cell numbers are further indicated in parentheses. This summary, quantified by the peak-to-baseline ratio as in Fig. 2C (individual cells, n = 42: gray markers; averages from A: red markers), confirms (C) our observations in single PCs (Fig. 2). The black curve represents the function (1 + e_(F_“)/b)'1, With best-fit parameters a = 50.5, b = 24.2 optimized over the set of 42 PCs. The inset further displays the location of the summary PRCs peak, relative to the time of the AP following the stimulus (i.e., Tpeak = ((ppeak — 1) - (181)), as in Fig. 2B. (D) The average PRCs are equivalently represented as a function of time (i.e., T = ((p — 1) - (181)) for the 20 ms preceding the perturbed AP.

(A) Direct estimate of the PRC, obtained using a step of current to control the firing rate of the cell (data analysed as in Fig. 2). Unlike what was seen in the experiments, the amplitude of Z((p) is bigger for lower rates suggesting that the same perturbation pulse advances a bigger percentage of the 181 at low firing rates. A small peak is observed at early (B) The peak at early phases induces a decrease in the peak-to-baseline ratio (the black curve represents the function (1 + e_(F ‘W b)_1, with best-fit parameters a = 19.0, b = 9.5) that however is not indicative of a flat PRC. Differently from the experiments, Z(phi) is larger at low firing rates, implying that the phase advance, in response to the same perturbation, is greater at lower firing rates. (EPS)

(A) Direct estimate of the PRC, obtained using a dendritic activation to control the firing rate of the cell. (B) Instantaneous rate of presynaptic firing rate during two inter-spike intervals of comparable duration (black and gray). Excitatory firing rate are represented by full lines and inhibitory firing rate by dashed lines. (C) Phase responses obtained from delivering perturbations during two (black and gray) inter-spike intervals of comparable duration. Dendritic conductances can greatly affect the shape of the PRC. (EPS)

The value of Z((p = 0.3) as a function of firing rate, both across PCs (open black circles) and for the same cell (filled gray squares). A gradual transition from positive to more negative values is apparent suggesting a smooth transition. The red markers represent the mean of all experiments. Black and gray curves the linear fit for all experiments and for the single cell respectively. (EPS)

H. Cuntz, A. Roth, S. Hong and T. Netoff for helpful discussions. We are also grateful to Q. Robberechts and

K. Veys and Y. Le Franc for help setting up the experiments and to Drs. M. Negrello and R. Man for comments on an earlier version of the manuscript.

Performed the experiments: IC. Analyzed the data: IC. Contributed reagents/materials/ analysis tools: IC DL EDS MG. Wrote the paper: IC DL MG. Performed computer simulations: IC DL.

Appears in 115 sentences as: firing rate (77) firing rates (62)

In *On the Firing Rate Dependency of the Phase Response Curve of Rat Purkinje Neurons In Vitro*

- Recently, the PRC theory applied to cerebellar Purkinje cells revealed that these behave as phase-independent integrators at low firing rates , and switch to a phase-dependent mode at high rates.Page 1, “Abstract”
- Given the implications for computation and information processing in the cerebellum and the possible role of synchrony in the communication with its post-synaptic targets, we further explored the firing rate dependency of the PRC in Purkinje cells.Page 1, “Abstract”
- We isolated key factors for the experimental estimation of the PRC and developed a closed-loop approach to reliably compute the PRC across diverse firing rates in the same cell.Page 1, “Abstract”
- Our results show unambiguously that the PRC of individual Purkinje cells is firing rate dependent and that it smoothly transitions from phase independent integrator to a phase dependent mode.Page 1, “Abstract”
- It has been shown that the PRC of tonically firing Purkinje Cells is flat at low firing rates , which has profound implications for information processing in the cerebellum.Page 1, “Author Summary”
- Here, we propose a novel method to estimate the PRC of single Purkinje cells at various firing rates and use it to unveil the smooth transition between flat and phasic PRC.Page 1, “Author Summary”
- The intrinsic electrical activity of Purkinje cells (PCs) exhibits a large repertoire of dynamical behaviors, including spontaneous firing of simple action potentials (APs), bistability of the firing rate , and hysteresis [1—4].Page 2, “Introduction”
- In addition, the extended range of PCs firing rates during behavior suggests that the rate of APs, its sudden transitions, its coherence across PCs, and the AP timing synchronization may contribute to information representation, processing, and downstream relaying.Page 2, “Introduction”
- Unexpectedly, they reported that the PC’s intrinsic firing rate has a profound effect on the response properties: the PRC of PCs firing at low rates displays a flat profile, suggesting that neurons behave like phase-independent inputs integrators; on the other hand, the PRC of PCs firing at high firing rates has a prominent peak, indicating a phase preference similar to coincidence detectors.Page 2, “Introduction”
- While it was not the first time that PRCs were shown to undergo changes over a range of AP frequencies [15] , the wide physiological range of PCs spontaneous firing rates and their ease of experimental access in in vitro preparations, made the report on the rate dependence of the PRC relevant.Page 2, “Introduction”
- In particular, in the light of the known bistable behavior of PCs, and their ability to abruptly toggle between distinct AP firing rates , we found it urgent to clarify whether the changes in PRC occur abruptly or smoothly, for increasing AP frequencies.Page 2, “Introduction”

See all papers in *March 2015* that mention firing rate.

See all papers in *PLOS Comp. Biol.* that mention firing rate.

Back to top.

Appears in 5 sentences as: computational model (1) Computational modeling (1) computational modeling (2) computational models (1)

In *On the Firing Rate Dependency of the Phase Response Curve of Rat Purkinje Neurons In Vitro*

- Using computational models we show that neither channel noise nor a realistic cell morphology are responsible for the rate dependent shift in the phase response curve.Page 1, “Abstract”
- Furthermore, we address potential explanations for the observed transition using computational modeling .Page 2, “Author Summary”
- Computational modelingPage 13, “Computational modeling”
- Active conductances are thought to modulate the shape of the PRC, and therefore computational modeling [7, 47, 48] could be a powerful tool to dissect the ionic bases of the PRC.Page 13, “Computational modeling”
- PRCs computed for a highly detailed computational model of a Purkinje cell with synaptic activation in the dendritic tree.Page 20, “Supporting Information”

See all papers in *March 2015* that mention computational modeling.

See all papers in *PLOS Comp. Biol.* that mention computational modeling.

Back to top.

Appears in 5 sentences as: single cell (3) single cells (2)

In *On the Firing Rate Dependency of the Phase Response Curve of Rat Purkinje Neurons In Vitro*

- The Phase Response Curve (PRC), a simple input-output characterization of single cells , can provide insights into individual and collective properties of neurons and networks, by quantifying the impact of an infinitesimal depolarizing current pulse on the time of occurrence of subsequent action potentials, while a neuron is firing tonically.Page 1, “Abstract”
- 2B, is comparable to that measured for single cells , revealing and confirming the same marked smooth dependency on the firing rate (Fig.Page 12, “Population summary”
- The phase response curve allows us to discuss the impact of single cell properties on net-work-level phenomena, as a putative way to relay strong inhibition downstream (by synchroni-sation of PCs).Page 19, “Implications for information processing”
- We provide a thorough characterization of single cell properties that might be involved in synchronization, hence contributing for the observed phenomena.Page 19, “Implications for information processing”
- Black and gray curves the linear fit for all experiments and for the single cell respectively.Page 20, “Supporting Information”

See all papers in *March 2015* that mention single cell.

See all papers in *PLOS Comp. Biol.* that mention single cell.

Back to top.

Appears in 5 sentences as: steady state (5)

In *On the Firing Rate Dependency of the Phase Response Curve of Rat Purkinje Neurons In Vitro*

- For the computation of the PRC, we evolved the model until it reached a steady state and then applied pulses of current (0.5 ms duration and 0.5 nA amplitude) at random times with a mean period between perturbations of 2.5 HZ.Page 8, “Modeling”
- In both cases, the algorithm employed to compute the PRC differed from that used in the single compartment model and resembled the one adopted in [42, 43]: briefly, after evolving the model until it reached a steady state , we identified two spike times to and t1 such that the 181 1‘1 — to was of the appropriate duration.Page 9, “Modeling”
- Similarly, it is often necessary to wait until PCs reach a steady state firing rate before initiating the repeated stimulation protocol.Page 9, “Results”
- 1A), employed to speedup convergence to the firing rate steady state and reduce very slow fluctuations.Page 9, “Results”
- Additionally, we set the simulation temperature to 28°C, in order to span a broader range of steady state firing frequencies, particularly in the low end of the spectrum.Page 18, “by”

See all papers in *March 2015* that mention steady state.

See all papers in *PLOS Comp. Biol.* that mention steady state.

Back to top.

Appears in 5 sentences as: synaptic inputs (6)

In *On the Firing Rate Dependency of the Phase Response Curve of Rat Purkinje Neurons In Vitro*

- By capturing the relationship between the evoked phase shift Ag) and the phase (p at which the input pulse occurred, the PRC predicts how, upon receiving weak synaptic inputs , neurons transiently delay or accelerate AP firing, contribute to network-wide AP synchrony, integrate external inputs or detect their temporal coincidences.Page 2, “Introduction”
- In the case of the model Without synaptic inputs we used the PMlO model [40], set the temperature of the simulation to 28°C and injected a somatic current of varying amplitude to span a Wide range of firing rates.Page 9, “Modeling”
- In the case of the model With synaptic inputs , we used the PM9 model With synapses distributed on the dendritic tree as described in [41], set the temperature to 37°C and fixed the presynaptic excitatory and inhibitory firing rates to 35 and 2 HZ, respectively.Page 9, “Modeling”
- To test whether distributed synaptic inputs might influence the shape of the PRC, we delivered random synaptic inputs to the dendrites of the model neuron.Page 16, “Computational modeling”
- We found that, even when the instantaneous rates of inhibitory and excitatory inputs were relatively stable, synaptic inputs greatly affected the shape of the PRC (S3 Fig, panel A).Page 16, “Computational modeling”

See all papers in *March 2015* that mention synaptic inputs.

See all papers in *PLOS Comp. Biol.* that mention synaptic inputs.

Back to top.

Appears in 4 sentences as: real-time (4)

In *On the Firing Rate Dependency of the Phase Response Curve of Rat Purkinje Neurons In Vitro*

- (A) A real-time closed-loop in vitro setup for experimentally estimating the PRC at a fixed firing rate, based on a Proportional-lntegral-Derivative (PID) controller and on a reactive-clamp paradigm, was employed to rapidly and optimally explore the dependency of the PRC of PCs on the cell’s firing rate.Page 5, “Whole-cell patch clamp electrophysiology and pharmacology”
- In short, an iterative estimate 13” of the cell’s instantaneous firing rate was updated online after each AP, detected in real-time (i.e., as a positive crossing of a voltage threshold), using the following formula: where T = 1s acts as the time scale over which the instantaneous firing rate is estimated, weighing each new AP and the previous firing history [21, 31].Page 7, “Closed-loop stimulation to achieve accurate frequency-clamp”
- We therefore designed and implemented via software a real-time closed-loop control system (Fig.Page 9, “Results”
- We specifically implemented a real-time system composed of a closed-loop controller of the cell instantaneous firing rate.Page 16, “Discussion”

See all papers in *March 2015* that mention real-time.

See all papers in *PLOS Comp. Biol.* that mention real-time.

Back to top.

Appears in 3 sentences as: current injection (3)

In *On the Firing Rate Dependency of the Phase Response Curve of Rat Purkinje Neurons In Vitro*

- In addition, cortical neurons fire APs only upon external holding current injection .Page 16, “PRC estimation methods”
- Nonetheless, biophysically accurate (i.e., voltage-dependent) channel flickering had a similar effect on the PRC as what was obtained by additive Gaussian noise current injection , in [14] and our simulations.Page 18, “by”
- We note that, in response to external current injection , the Khaliq-Raman model fires at arbitrarily low firing rates, implying a type I excitability.Page 18, “by”

See all papers in *March 2015* that mention current injection.

See all papers in *PLOS Comp. Biol.* that mention current injection.

Back to top.

Appears in 3 sentences as: experimental data (3)

In *On the Firing Rate Dependency of the Phase Response Curve of Rat Purkinje Neurons In Vitro*

- Having observed a rate dependency in the CV of PCs in our experimental data (81 Fig, panel C) we asked Whether spiking variability could account for the flat PRC profile observed at low firing rates.Page 14, “Computational modeling”
- However, while their main focus was to introduce and test a novel technique for the PRC estimation, their experimental data set consisted of a small population of 16 cells.Page 16, “Discussion”
- In our experimental data set, the inter-spike intervals distributions displayed a significant inverse correlation of the coefficient of variation (CV) with the firing rate (i.e., Pearson’s r = —0.4, p < 10—6, and a slope of —0.25/kHz—Sl Fig, panel B) regardless of whether open or closed-loop methods used.Page 17, “Conductance-based neuron modeling”

See all papers in *March 2015* that mention experimental data.

See all papers in *PLOS Comp. Biol.* that mention experimental data.

Back to top.

Appears in 3 sentences as: input-output (3)

In *On the Firing Rate Dependency of the Phase Response Curve of Rat Purkinje Neurons In Vitro*

- The Phase Response Curve (PRC), a simple input-output characterization of single cells, can provide insights into individual and collective properties of neurons and networks, by quantifying the impact of an infinitesimal depolarizing current pulse on the time of occurrence of subsequent action potentials, while a neuron is firing tonically.Page 1, “Abstract”
- Recently, key results from the mathematical theory of coupled oscillators sparked a lot of interest: a simple input-output characterization of the units composing a network, known as their phase response (or phase resetting) carve (PRC), is sufficient to classify and predict individual and collective properties.Page 2, “Introduction”
- We employed the current-clamp configuration and focused on the input-output relationship between the phase shift Ag), induced in the cell firing cycle by an external current pulse, and the phase (p at which the pulse was repeatedly delivered (Fig.Page 9, “Results”

See all papers in *March 2015* that mention input-output.

See all papers in *PLOS Comp. Biol.* that mention input-output.

Back to top.