JN Add DOIs to your references at manuscript stage!
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
 QUICK SEARCH:   [advanced]


     


J Neurophysiol 92: 977-996, 2004. First published March 24, 2004; doi:10.1152/jn.00067.2004
0022-3077/04 $5.00
This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow All Versions of this Article:
92/2/977    most recent
00067.2004v1
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Right arrow Citation Map
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in Web of Science
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Web of Science (24)
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Giugliano, M.
Right arrow Articles by Streit, J.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Giugliano, M.
Right arrow Articles by Streit, J.

Single-Neuron Discharge Properties and Network Activity in Dissociated Cultures of Neocortex

M. Giugliano1, P. Darbon1, M. Arsiero1,2, H.-R. Lüscher1 and J. Streit1

1Institute of Physiology, University of Bern, 3012 Bern, Switzerland; and 2Department of Biophysical and Electronic Engineering, University of Genoa, 16145 Genoa, Italy

Submitted 20 January 2004; accepted in final form 17 March 2004


 ABSTRACT
 
 TOP
 ABSTRACT
 INTRODUCTION
 METHODS
 RESULTS
 DISCUSSION
 GRANTS
 ACKNOWLEDGMENTS
 REFERENCES
 
Cultures of neurons from rat neocortex exhibit spontaneous, temporally patterned, network activity. Such a distributed activity in vitro constitutes a possible framework for combining theoretical and experimental approaches, linking the single-neuron discharge properties to network phenomena. In this work, we addressed the issue of closing the loop, from the identification of the single-cell discharge properties to the prediction of collective network phenomena. Thus, we compared these predictions with the spontaneously emerging network activity in vitro, detected by substrate arrays of microelectrodes. Therefore, we characterized the single-cell discharge properties to Gauss-distributed noisy currents, under pharmacological blockade of the synaptic transmission. Such stochastic currents emulate a realistic input from the network. The mean (m) and variance (s2) of the injected current were varied independently, reminiscent of the extended mean-field description of a variety of possible presynaptic network organizations and mean activity levels, and the neuronal response was evaluated in terms of the steady-state mean firing rate (f). Experimental current-to-spike–rate responses f(m, s2) were similar to those of neurons in brain slices, and could be quantitatively described by leaky integrate-and-fire (IF) point neurons. The identified model parameters were then used in numerical simulations of a network of IF neurons. Such a network reproduced a collective activity, matching the spontaneous irregular population bursting, observed in cultured networks. We finally interpret such a collective activity and its link with model details by the mean-field theory. We conclude that the IF model is an adequate minimal description of synaptic integration and neuronal excitability, when collective network activities are considered in vitro.


 INTRODUCTION
 
 TOP
 ABSTRACT
 INTRODUCTION
 METHODS
 RESULTS
 DISCUSSION
 GRANTS
 ACKNOWLEDGMENTS
 REFERENCES
 
In vivo neocortical neurons receive a continuous barrage of excitatory and inhibitory postsynaptic potentials (EPSPs/IPSPs) (Destexhe and Paré 1999Go). Such an intense background activity arises from the intracortical presynaptic neurons, spontaneously spiking at low rates, and it induces the postsynaptic membrane voltage to fluctuate as in a random walk, resulting in an irregular spike emission (Destexhe et al. 2003Go). Under such conditions, the biophysical properties of the neurons are substantially altered compared to those of a neuron at rest (see Bernander et al. 1991Go and references therein; Destexhe et al. 2003Go). Such considerations define a different scenario with respect to traditional noise-free in vitro experimental para-digms, used either in the characterization of the integrative properties of single neurons or in the investigation of emerging network phenomena. As a consequence, observations carried out in vitro may not directly transfer to the in vivo situations (Steriade 2001Go).

Even more important, the characterization of the single-neuron response properties, performed under appropriate and realistic conditions, is a key element to predict and understand how a population of neurons collectively interacts and processes information in the intact brain. This view is supported by several theoretical studies (Amit and Tsodyks 1992Go; Brunel 2000Go; Mattia and Del Giudice 2002Go; Salinas 2003Go), where the single-neuron response properties were used to make predictions about the collective phenomena, such as the global spontaneous irregular activity (Amit and Brunel, 1997bGo), the emergence of network-driven oscillations (Brunel and Wang 2003Go; Fuhrmann et al. 2002Go), and of selective delay-activity states (Amit and Brunel 1997bGo; Wang 2001Go; Yakovlev et al. 1998Go).

Under these perspectives, novel electrophysiological paradigms were recently proposed by several experimenters, focusing on the integrative (Destexhe and Paré 1999Go, 2000Go; Poliakov et al. 1997Go; Rauch et al. 2003Go), computational (Chance et al. 2002Go; Mainen and Sejnowski 1995Go; Protopapas and Bower 2001Go), and adaptive properties of single neurons (Fuhrmann et al. 2002Go; Paninski et al. 2003Go) and synapses (Froemke and Dan 2002Go; Sjöström et al. 2001Go). In particular, Rauch et al. (2003)Go proposed to mimic the conditions that neurons experience in the recurrent networks of the intact cortex. They reproduced a realistic background synaptic activity as a computer-synthesized noisy current that was injected into the soma, as shown in Fig. 1, A and B (see also Destexhe and Paré 2000Go). Under such conditions, layer V pyramidal neurons in acute neocortical slices respond as integrate-and-fire (IF) point neurons (Rauch et al. 2003Go). This suggests that an extremely simplified model neuron may be used to accurately describe the discharge properties of a single cell under simulated in vivo conditions and to predict and interpret collective emergent phenomena. Unfortunately, because it is not straightforward to observe, induce, or re-create coordinated network activity in neocortical slices (but see Giugliano and Lüscher 2003Go; Giugliano et al. 2003Go; and Shu et al. 2003Go), there was no attempt at validating any network-level prediction in vitro, resulting from the proposed minimal characterization of the single cell. In addition, a direct confirmation of the same predictions in awake andbehaving animals is at present impossible or still requires extremely challenging technical issues to be overcome.



View larger version (69K):
[in this window]
[in a new window]
 
FIG. 1. As opposed to conventional DC stimulation (A, top trace), the electrical response evoked by a nondeterministic input current (B), mimicking a realistic input drive of cortical networks, is determined by either the steady current or the random fluctuations. Each stimulus waveform was generated as an independent realization of the stochastic Ornstein–Uhlenbeck process (C), fully specified by the steady-state mean m, variance s2, and autocorrelation time length {tau}I, and characterized by a stationary Gaussian amplitude distribution and an exponentially decaying autocorrelation function Rii (see insets). Evoked time course of the membrane voltage, corresponding to the current plotted in C, resembled a random walk with irregular spike emission (D, top trace), similar to the intracellular activity recorded in vivo. For each pair (m, s2), the neuronal response was quantified in terms of the mean firing rate (D, bottom panel) estimated at the steady state (the reported instantaneous firing rate reaches its steady state after 2–5 s as it adapts exponentially, with a time constant of 2 s. Experiments were performed under current-clamp conditions. Current was injected into the soma of E cultured dissociated neocortical neurons (scale bar: 50 µm).

 
For such reasons, we decided to focus on the quantitative characterization of the discharge properties of neocortical cultured neurons, dissociated from the neocortex. Such a preparation has been widely used for investigating single-neuron properties as well as the network activity, which spontaneously emerges in vitro (see for a review Marom and Shahaf 2002Go). In particular, developing networks of synaptically connected cortical neurons can be grown and cultured on arrays of substrate microelectrodes (MEAs), so that the distributed electrical activity can be recorded for a long time, with a high spatial resolution (Potter 2001Go).

In this paper, we show that mature cultured neurons retain the features that characterize neurons in acute brain slices, investigated under the same experimental paradigm. Under these conditions, they respond as IF point neurons, not only with respect to the output mean firing rate but also to the subthreshold membrane voltage statistics, confirming and extending the work of Rauch et al. (2003)Go. As a step toward the quantitative description of the mechanisms underlying population phenomena in cultured networks (Jimbo et al. 1999Go; Kamioka et al. 1996Go; Maeda et al. 1995Go; Marom and Shahaf, 2002Go), we computer-simulated the electrical activity of a network of model neurons, incorporating the single-neuron details identified in the experiments. We further interpret the results of our simulations by the extended mean-field theory (Renart et al. 2003Go), and we compare them with in vitro experimental recordings, using the MEAs.


 METHODS
 
 TOP
 ABSTRACT
 INTRODUCTION
 METHODS
 RESULTS
 DISCUSSION
 GRANTS
 ACKNOWLEDGMENTS
 REFERENCES
 
Cultures

Cultures of neocortical neurons were obtained from the somatosensory/motor cortex of newborn Wistar rats (P1-2), following standard procedures. Rats were anesthetized (0.4 ml Vetanarcol) and killed by decapitation. Transverse brain slices, 225 µm thick, were exposed to a 0.3% trypsin solution for 3 min at 37°C for enzymatic digestion. Cells were mechanically dissociated and plated on substrate arrays of planar microelectrodes (MEAs) or glass coverslips, at a density of 150,000 or 75,000/150 µl, respectively. MEAs were microfabricated as described previously (Tscherter et al. 2001Go) and coated for 1 h with diluted (1:50) Matrigel (Falcon/Biocoat, Becton Dickinson AG, Basel, Switzerland). The glass coverslips were coated with polylysine (1 mg/ml overnight at 37°C or Matrigel (1:50). Cells were restricted to a small area ({approx}50 mm2) using cloning glass cylinders attached to the MEAs or coverslips by silicone sealant. Cultures were incubated at 36.5°C in a 5% CO2 air atmosphere and maintained in 150 µl nutrient medium. The medium contained MEM Eagle (Sigma) supplemented with fetal bovine serum 10%, glucose 0.2%, B27, and Glutamax (Gibco BRL, Life Technologies AG, Basel, Switzerland) for the cultures on glass coverslips, and serum-free Neurobasal medium supplemented with B27 and Glutamax (Gibco BRL) for the cultures on the MEAs. Half of the medium was changed weekly.

Recordings were made in a chamber mounted on the stage of an upright microscope (Nikon, Tokyo, Japan). Patch-clamp experiments were carried out on cells in which a few processes could be visually identified (e.g., Fig. 1E), from cultures of 4–7 wk of in vitro age. Within this period, no systematic change of single-neuron and network properties was detected, suggesting a complete maturation of the neurons and indicating that the cultures' development had reached a steady state (Kamioka et al. 1996Go; Marom and Shahaf 2002Go).

Pharmacology

Several minutes before starting each recording session, the culture medium was replaced by an extracellular solution containing (in mM): NaCl 145, KCl 4, MgCl2 1, CaCl2 2, HEPES 5, Na-pyruvate 2, glucose 5, at pH 7.4. Recordings were performed either in the absence of a continuous flow of solution, with solution changes every 10–15 min, or in the presence of continuous superfusion at 1 ml/min. No differences were detected between these protocols. All recordings were made at room temperature (23–24°C).

As extensively described in the literature, mature cultures of dissociated neocortical neurons exhibit a spontaneous electrical activity, which results in simultaneous bursting network activity and mainly depends on the synaptic development and coupling between the neurons (Jimbo et al. 1999Go; Kamioka et al. 1996Go; Maeda et al. 1995Go; Marom and Shahaf 2002Go). In the patch-clamp experiments, given that our goal was to characterize the single-neuron response properties to a class of computer-synthesized stimulus waveforms, single-neuron recordings were performed while blocking glutamatergic synaptic transmission with D-APV (D-2-amino-5-phosphonovalerate 50 µM), a competitive antagonist of the NMDA (N-methyl-D-aspartate) receptors, and CNQX (6-cyano-7-nitroquinoxaline-2-3-dione, 10 µM) (both Sigma, Buchs, Switzerland), a competitive antagonist of non-NMDA receptors. These substances were bath applied and completely suppressed incoming synaptic activity and spontaneous spiking in the neurons. In the experiments involving multisite extracellular recordings of the neuronal electrical activity by the MEAs, an antagonist of GABAA receptors (bicuculline methochloride, 10 µM) (Tocris, Anawa Trading SA, Wangen, Switzerland) was bath applied to block synaptic inhibition, with the aim of focusing on the disinhibited pattern generation as well as establishing a more direct comparison to the computer-simulated networks, which included only excitatory model neurons.

Whole cell patch-clamp recording from single neurons

The patch-clamp technique was used, in the whole cell configuration (Hamill et al. 1981Go), by using an Axoclamp-2B amplifier (Axon Instruments, Union City, CA). Signals were recorded in current clamp, filtered at 1 kHz, sampled at 5 kHz, and digitized by a 12-bit A/D converter (Digidata 1200) and pClamp 8 software (Axon Instruments). Electrodes were pulled from filamented borosilicate glass capillaries (GC150F, Harvard Apparatus GmbH, March-Hugstetten, Germany) on a horizontal puller (DMZ, Zeitz Instrumente GmbH, Munich, Germany) and their resistance was 5.5 ± 0.6 M{Omega}. Electrodes were filled with a solution containing (in mM: K-gluconate 100, KCl 20, HEPES 10, Mg-ATP 4, Na2-GTP 0.3, Na2 phosphocreatine 10; pH 7.3, 290 mOsm. Other (standard) pipette solutions were reported not to alter significantly the response properties of the cells, under the very same experimental protocol (Rauch et al. 2003Go). High-resistance seals (2–4 G{Omega}) were formed and the whole cell configuration was achieved by the application of a negative pressure pulse. The bath application of D-APV and CNQX always followed the establishment of the whole cell patch configuration.

Conventional procedures, consisting of repetitive hyperpolarizing step currents, were used to obtain a direct estimate of the passive (cable) properties of cultured neurons (Iansek and Redman 1973Go). Such estimates are related to the time constant {tau}m and capacitance Cm of an equivalent lumped RC compartment (Abbott and Dayan 2001Go).

Emulating a stationary realistic input from the network

Independent realizations of the stochastic Ornstein–Uhlenbeck process were computer-synthesized and injected under current clamp (see Fig. 1C). (Interactive stimuli-synthesis and analysis software tools were developed and are available on request.) Such nondeterministic current stimuli mimic a realistic barrage of excitatory/inhibitory postsynaptic currents (EPSCs/IPSCs) for a cell embedded in a large in vivo network, spontaneously and randomly active at low rates (Destexhe et al. 2003Go), as well as in a cultured network where spontaneous neurotransmitter release at individual synapses, as well as other sources of inhomogeneity and randomness determine an irregular background synaptic noise in vitro (Maeda et al. 1995Go). Although a conductance injection would have appeared more appropriate (e.g., by means of the dynamic-clamp technique; see Destexhe and Pare 2000), it has been proven that noisy conductance-driven and noisy current-driven stimulations are equivalent with respect to the evoked neuronal steady-state mean firing rates, apart from 2 stimulus-independent scaling factors, for the input mean and variance (La Camera et al. 2003Go; Rauch et al. 2003Go). More important, the comparison of the experimental data to the predictions of mathematical model neurons has been carried out under the same current-driven conditions. Under these hypotheses and because the impact of a single EPSC/IPSC on the postsynaptic membrane voltage is weak in evoking suprathreshold responses in vitro (Nakanishi and Kukita 1998Go; Nakanishi et al. 1999Go), the overall current experienced by a generic postsynaptic neuron can be approximated as a diffusion stochastic process (Destexhe and Paré 2000Go; Fourcaud and Brunel 2002Go). Therefore, under extended mean-field hypotheses on the interactions between neurons in a population (Amit and Brunel 1997bGo), such noisy stimuli may statistically account for a wide range of feedforward/recurrent network architectures and regimes. For instance, by indicating with Ne/i the number of excitatory/inhibitory afferents with stationary mean activation rates fe/i, under the hypothesis that synaptic inputs are approximately independent, the distribution of the resulting postsynaptic somatic current amplitude becomes Gaussian (i.e., by the central-limit theorem), with steady-state mean m and variance s2 given by the expressions reported below (Amit and Brunel 1997bGo)

where Ie/i and {tau}e/i are the effective peak-amplitude and decay time constant at the soma, for individual excitatory and inhibitory postsynaptic currents, respectively. Throughout this work, we set {tau}e = {tau}i = {tau}I {1; 5} ms, thereby mimicking (AMPA- and GABAA-mediated) fast synaptic currents (Destexhe et al. 1994Go).

We note that the features of the incoming synaptic current arising from the particular spike-timing precision and reliability of individual presynaptic neurons were not explicitly investigated (Mainen and Sejnowski 1995Go; Jolivet et al. 2004). Consistently, the neuronal response to such noisy current input was routinely analyzed by characterizing the output mean spike rate at the steady state. However, provided that the hypothesis on the statistical independence of the presynaptic activity holds, the impact of neuronal precision and reliability is expected to weakly contribute to the network activity we discuss in the present work. Spontaneous synaptic release and other sources of network randomness will in fact introduce uncorrelated stochastic components to the overall synaptic inputs to any neuron. Therefore, under such an hypothesis, the overall resulting current incoming to a generic postsynaptic neuron of the network is still Gauss-distributed and completely characterized by mean m and variance s2.

Noisy current-clamp protocol

The stimulation protocol consisted of the repeated somatic injection of independent current realizations I(t), each lasting 20 s and interleaved by 30–60 s of recovery time. For any pair (m, s2), the following iterative expression was used to synthesize a realization I(t) of the process (Cox and Miller 1965Go)

(1)
where {xi}t is a unitary Gauss-distributed random variable (Press et al. 1992Go), updated at every time step at a rate of 5 kHz (i.e., dt = 0.2 ms) (Fig. 1C). The exploration of the plane (m, s2) typically requested 30–45 min to be completed and it was carried out in a shuffled order, randomly repeating each pair twice. During such a time interval, the lack of spontaneous synaptic events, the stability of the resting membrane potential, and the cell input resistance were continuously monitored. The entire procedure was stopped in case of drifts in any of these observables.

Analysis of the single-neuron response properties

Collected data consisted of the voltage responses to each noisy stimulus, lasting 20 s and associated with the pair (m, s2) (see Fig. 1D). To account for current-clamp offsets, the actual injected current was also monitored and the actual value of (m, s2) directly estimated from it. Raw voltage traces were processed in Matlab (The MathWorks, Natick, MA) by a peak-detection algorithm to extract individual spike times and shape. When no substantial change in the shape of the individual action potentials, or in the instantaneous firing rate at the beginning and at the end of the elicited spike train occurred (representative of possible nonstationarities), the trial was accepted. The steady-state mean spiking frequency f was estimated, discarding an initial transient (i.e., 2–5 s) and averaging across the 2 available repetitions. After a successful completion of the stimulation protocol, the experimental curve {fh = fh(mh, sh), h = 1, 2,... , M} (for at least 3 distinct values of s and typically M ~ 20) was plotted and compared to the theoretical responses of 2 model neurons, driven by the same current statistics (Fig. 2). The recordings were further compared to the data available for pyramidal layer V neurons in acute slices, collected under the same experimental conditions (Rauch et al. 2003Go).



View larger version (54K):
[in this window]
[in a new window]
 
FIG. 2. A: Sketch of the electrical equivalent circuits of the single-compartment models considered for the fit of the experimental responses: for the sake of simplicity neither the absolute refractoriness nor the spike-frequency–dependent adaptation was represented. The condition V(t) = {theta} corresponds to the emission of an action potential at time t. For both models, under the current-driven stimulation protocol (see METHODS), a description of the steady-state discharge properties can be analytically derived. B: the response functions (Eqs. 2 and 3) were plotted for increasing fluctuations amplitude s, including (i.e., aLIF, aCLIFF) or not including (i.e., LIF, CLIFF) the spike-frequency–dependent adaptation. Adaptation generally decreases the spike rates, while preserving the sensitivity to fluctuations (s2), as discussed by La Camera et al. (2002)Go. The plots reveal a different dependency on s in the two models, evident at small and large s (50–800 pA, step 50 pA).

 
Mathematical model of the single-neuron response

The mathematical modeling and the procedures for the model fit to the collected response data are standard and closely followed Rauch et al. (2003)Go. Briefly, 2 single-compartment standard mathematical descriptions of neuronal excitability were considered: the Lapicque's or leaky integrate-and-fire neuron LIF (Abbott and Dayan 2001Go; Tuckwell 1988Go) and the constant leakage integrate-and-fire neuron with a floor (CLIFF) (Fusi and Mattia 1999Go; Mongillo and Amit 2001Go; Rauch et al. 2003Go) (see Fig. 2). These models have been studied in depth (Fourcaud and Brunel 2002Go; Mattia and Del Giudice 2002Go) and widely used in simulations of large-scale networks (Mattia and Del Giudice 2000Go; Reutimann et al. 2003Go) and hardware implementations (Chicca et al. 2003Go). As opposed to the biophysically realistic conductance-based models (Abbott and Dayan 2001Go), these descriptions are characterized by a single state variable V (i.e., the membrane potential) and by a reduced set of effective constant parameters: the absolute refractory period {tau}arp and reset voltage H, the membrane capacitance C, the voltage threshold for spike emission {theta}, and the subthreshold voltage decay rate –C(EV)/{tau} for the LIF ({tau} being the membrane effective time constant and E the resting potential) or –{lambda} for the CLIFF model. Finally, both models incorporated a simplified spike-frequency–dependent adaptation, modeling the contribution of intracellular calcium- and/or sodium-activated outward currents to the net membrane current, modulated by a constant factor {alpha} and implemented as described in Liu and Wang (2001)Go and van Vreeswijk and Hansel (2001)Go. The stationary effect of such an adaptation is to reduce the gain of the frequency–current curve by a factor that does not depend on the adaptation time constant {tau}a (see Eqs. 2 and 3). We refer as aLIF and aCLIFF, to the LIF and CLIFF models incorporating the spike-frequency–dependent adaptation, respectively (see Fig. 2B).

Model parameters fit to the experimental data

The full analysis of the response of the spiking neuron models to the input current, specified by Eq. 1 (i.e., a colored noise), can be performed only under approximate treatments (Fourcaud and Brunel 2002Go) but it was not considered in the present work. Instead, it is considerably easier to analyze the models' response to an idealized "equivalent" white current (i.e., {delta}-correlated), characterized by the infinitesimal mean µ and variance {sigma}2 determining an asymptotically equivalent effect on the subthreshold membrane voltage V (Rauch et al. 2003Go). The last is a satisfactory approximation, when the evoked mean interspike interval is much larger than {tau}I, as verified by computer simulations. Under these hypotheses, the steady-state current-to-rate response function {Phi}(m, s2) can be analytically determined (Fusi and Mattia 1999Go; Tuckwell 1988Go) (Fig. 2B) and the corresponding parameter space may be searched for the best fit to the stationary data points {fh(mh, sh), h = 1, 2,... , M}, collected for each cell. The implicit formulas corresponding to the mean firing rate {Phi}aLIF/aCLIFF under noisy current stimulation are

(2)

(3)
indicating for the aLIF model


with {tau} = RC the membrane time constant (Fig. 2A), whereas for the aCLIFF model

 Without loss of generality, all the voltages were referred to 0 mV and the spike-emission threshold {theta} was kept fixed at 20 mV, with respect to the resting membrane potential (i.e., E = 0 mV and {theta} = E + 20 mV).

For each cell, simulated-annealing optimization techniques (Press et al. 1992Go) were used to fit the above reported theoretical response functions to the data, by minimizing the following mean-square error with respect to the model parameters

Each term of the sum is weighted by the inverse of an accuracy interval, estimated by using a simplified phenomenological binomial model of the spike-emission processes (Rauch et al. 2003Go). Briefly, by specifying a desired confidence level q on the spike rate estimate {fh, h = 1, 2,... , M} and indicating with Nsph the number of spikes emitted in the time interval T as a realization of a binomial random variable, the mean rate f and its experimental estimate fh = Nsph/T satisfy the following relationship

with

The last relationships imply that, by setting K ≥ 1, the resulting accuracy interval {delta}h± for fh corresponds approximately to a confidence of at least 68%.

Finally, {chi}2 and its minimum {chi}min2 are random variables, known to be approximately distributed according to a {chi}2 distribution (Press et al. 1992Go). This makes it possible to refer to the probability Prob ({chi}2 > {chi}min2) as a standard indication of the goodness of the fit. Model fits were accepted when such a probability was >0.1. The minimal value of K of successful {chi}2-test was chosen as a quantitative analogue measure of the quality of the fit, comparing the performances of the aLIF, LIF, aCLIFF, and CLIFF models over the entire data set. Finally, the parameters search was repeated twice for each experiment: under free-search conditions and introducing additive quadratic cost-penalty constraints to the simulated annealing energy (Press et al. 1992Go). This aimed at discouraging the exploration of the parameter space in a limit-valued region (i.e., H -> {theta} and {tau}arp -> +{infty}). We refer to these conditions as free search, or no-penalty best fit, and as penalty best fit.

MEA recording and analysis of the network activity

In some experiments, MEAs were used as a substrate for the cultured networks, so that neuronal somata as well as axons were allowed to develop close to the individual MEA microelectrodes. Such a proximity makes possible the extracellular detection of the emission of action potentials by one or more neighboring cells (Streit et al. 2001Go). MEAs contained 68 platinum planar electrodes, spaced at 200-µm intervals and laying out in the form of a rectangle. Recording channels showing activity were selected by eye and their recordings digitized at 6 kHz per channel and stored on a hard disk. The detection of extracellularly recorded action potentials (i.e., fast voltage transients) and further analysis were performed off-line in IGOR. (WaveMetrics, Lake Oswego, OR), as described previously (Streit et al. 2001Go; Tscherter et al. 2001Go). No attempt was made to sort spikes detected by the same MEA electrode. The electrical noise of individual channels was very stable. The selectivity of event detection was routinely checked by using recordings obtained in the presence of tetrodotoxin (TTX, 1.5 µM, Sigma), as a "zero" reference.

The spontaneous network activity, detected by the MEA, consisted of asynchronous activity and population bursts (PBs). The asynchronous activity was defined as one or a few spikes detected by a single or several electrodes, whereas the PBs, consisting of episodes of activity occurring simultaneously across several recording channels, spread over the network. The processed network activity was visualized in the form of event raster plots and of the instantaneous population mean firing rate (see Fig. 5A). The last was computed by counting the total number of detected events from all recording channels, within a sliding time window of 10 ms. With the aim of summarizing the features of the detected network activity, the mean and the coefficient of variation (CV) of the interburst intervals (IBI) distribution and of the burst durations (PBd) distribution were estimated over at least 10 min of continuous recording (15 PBs minimum).



View larger version (49K):
[in this window]
[in a new window]
 
FIG. 5. The spontaneous emergence of patterned network activity: microelectrode array (MEA) recordings and computer simulations. A: disinhibited network activity, emerging in dissociated cultures of neocortical neurons, was detected by the MEAs. The raster plot indicates the time of occurrence of the events extracellularly detected by 7 substrate electrodes, and the resulting population mean firing rate (bottom panel) was estimated as described in METHODS (scale bars: 60 s, 10 Hz). Both B and C contain the raster plot (top panels) of the spikes emitted by a subset of the neurons of a simulated IF network, together with the population mean firing rate (bottom panels). The gray shading reports the instantaneous adaptation level, averaged over all the model neurons. The dotted horizontal gray lines represent the prediction from the mean-field theory for the intraburst firing rate fburst and the resting activity frest (see Fig. 8B) (scale bars: 60 s, 30 Hz). Depending on the synaptic coupling between neurons, the simulated network bursting occurs at a frequency of about 0.016 Hz (B), determined by fluctuations and by the probability of spontaneous network transitions (see Fig. 7 and the RESULTS), or at a frequency of 0.05 Hz (C), imposed by the relaxation of the adaptation mechanisms ({tau}a). Spontaneous synaptic release or other sources of randomness, together with the finite-size fluctuations in the network collective activity, represent an effective mechanism to considerably slow down oscillation frequency to a time scale that is much longer than the recovery time constants {tau}a. D: to compare quantitatively the simulation results to real network data, the results from 4 MEA experiments under bicuculline were plotted, by simultaneously studying the CV and the mean of the interburst intervals (IBIs) distribution and of the duration of the population bursts (PBd). Numbers associated with each marker help to identify the same experiment in both plots and indicate the number of PBs, collected in 10 min of continuous recording and considered for the statistics. Computer simulations of the network of aLIF model neurons reproduced the same features of the spontaneous patterned activity, simultaneously matching the IBIs and the PBd, while sweeping Je from 12 to 18 pA (indicated by the arrows). Interestingly, the statistics of the simulated IBIs resembled those of a Poisson point process with dead time, whose functional relationship relating coefficient of variation and mean of the distribution, is reported as CVIBIs = (meanIBIsTdead)/meanIBIs. Parameters used were: Ne = 100, Cee = 0.38, {tau}e = 5 ms, {tau}a = 2.1 s, {alpha} = 0.75 pA·s, m0 = 25.1 pA, s0 = 92 pA, Tdead = 6.1 s; B, Je = 13 pA; and C, Je = 16 pA.

 
Mathematical model of a network of IF neurons

We studied and computer-simulated the collective electrical activity of small networks of Ne = 100–1,000 interacting excitatory aLIF identical neurons (van Vreeswijk and Hansel 2001Go), using the single-neuron effective parameters identified in the experiments (Table 1). We chose to neglect network inhomogeneities to keep the interpretation of the results as simple as possible. Anyway, our interpretative framework can be extended to cover a similar situation, following the approach described in Amit and Brunel (1997a)Go.


View this table:
[in this window]
[in a new window]
 
TABLE 1. Numerical values of passive membrane properties {tau}m and Cm, aLIF effective parameters, corresponding fit performances, and number of cells n that passed the fit test, reported for penalty-fit and free-search cases

 
In the model, each excitatory neuron is characterized at time t by its membrane potential Vi(t) and its adaptation current Xi(t). These quantities follow the dynamical equations


Indicating with tki the times of emission of a spike by the ith neuron [i.e., V(tki) = {theta}], the previous equations are complemented by the conditions


whereas, indicating with Cij the connectivity matrix and by {Theta}(t) the unitary step function (i.e., {Theta}(t) = 0, t < 0 and {Theta}(t) = 1, t > 0), the total synaptic current Ii(t) into the neuron ith is given by

Focusing on cultured neocortical networks, we chose an unstructured connection topology (Marom and Shahaf 2002Go) with a probability Cee of synaptic connection between any 2 neurons of 0.3–0.4 (i.e. Prob {Cij = 1|i != j} = Cee) (Nakanishi and Kukita 1998Go; Nakanishi et al. 2001Go). In agreement with the hypotheses underlying the stimuli injected in the experiments, the synaptic interactions were described by currents rather than conductance changes. The adaptation currents were also described as current changes, as well. Describing these variables by conductance changes does not qualitatively affect the results reported. Synaptic interactions between 2 connected neurons were triggered by the presynaptic emission of action potentials, after an effective delay {delta} of 1.5 ms, which included the axonal propagation delay and synaptic release latency (Nakanishi and Kukita 1998Go). The resulting individual postsynaptic currents were characterized by an instantaneous rise to Je and by an exponential decay with a time constant {tau}e = 5 ms, and no activity-dependent short- or long-term change. The effect of spontaneous synaptic release and of other sources of randomness (Maeda et al. 1995Go) was incorporated into an activity-independent additional random synaptic drive. The effective value of {alpha} unavoidably included the superimposed stationary contribution of both fast and slower adaptation mechanisms, estimated in the patch-clamp experiments over 20 s of stimulation time. With the aim of focusing on the time scales characterizing the network bursting (i.e., ~0.1–1 s), sometimes we decreased the value of {alpha} that affected the PBd and not the IBIs statistics. For the sake of simplicity, the intrinsic cumulative inactivation experimentally measured was not included in the simulations (but see Giugliano et al. 2002Go).

Under the same extended mean-field hypotheses that underlie the noisy currents that were injected into real neurons, the collective activity of the simulated network can be fully predicted in terms of its stationary mean firing rate f, by the knowledge of {Phi}aLIF(m, s2) and the details of the synaptic connectivity (Amit and Brunel 1997bGo). In the present case, the statistics of the total recurrent synaptic current experienced by a generic neuron of the network is approximately Gauss-distributed and it can be fully described by its mean and variance

(4)
where f is the network mean firing rate and Je is the effective peak-amplitude for the individual postsynaptic currents, uniform across all the synapses of the network. m0 and s02 are 2 fixed parameters, corresponding to the spontaneous neurotransmitter release and other sources of randomness, independent of f.

Statistics

The Pearson's r linear correlation as well as the Kendall's Tau nonparametric (rank-order) test (Press et al. 1992Go) were used to assess statistical correlations. The last provides a correlation measure together with an estimated significance level pK, which corresponds to the probability of obtaining the same correlation from statistically independent samples. Averages are expressed as means ± SE.


 RESULTS
 
 TOP
 ABSTRACT
 INTRODUCTION
 METHODS
 RESULTS
 DISCUSSION
 GRANTS
 ACKNOWLEDGMENTS
 REFERENCES
 
Cultured neurons respond as simple IF neurons

As opposed to a DC stimulation (see Fig. 1A), under noisy current injection, the neuronal membrane voltage evolves in time as in a random walk, leading to irregular spike emission (see Fig. 1, A and B). Its subthreshold amplitude distribution becomes bell-shaped, with mean and SD increasing with the steady-state mean m and variance s2 of the injected current, respectively (see Figs. 1D and 4).



View larger version (40K):
[in this window]
[in a new window]
 
FIG. 4. The steady-state amplitude distribution of the subthreshold membrane voltage (open circles) was in good agreement with the aLIF model prediction (continuous thick line), although its parameters were tuned to match the mean firing rates only. AD: the panels indicate the results from 4 different experiments, and in each subplot the voltage amplitude density distribution is plotted for increasing values of m, from left to right, and of s from top to bottom, within the same panel. For each cell, the scaling factors (i.e., V -> {beta}VE) were fixed across the (m, s2) plane and their values were: A, E = 53.14 mV, {beta} = 0.36; B, E = 49.27 mV, {beta} = 0.51; C, E = 49.02 mV, {beta} = 0.64; and D, E = 55.32 mV, {beta} = 0.67. The values of E weakly correlated to the experimentally measured resting membrane potentials, estimated as: A, E = –59.24 mV; B, E = –56.25 mV; C, E = –62.19 mV; and D, E = –55.9 mV (r = 0.47; Tau = 0.67, pK = 0.17).

 
With respect to the steady-state spiking frequency, cultured neocortical neurons respond as simple IF neurons (Fig. 2), similarly to pyramidal neurons (Rauch et al. 2003Go) and with a similar response gain. Actually, the curves f(m, s2), collected for each cell, were quantitatively fitted with high accuracy by the IF model response functions, in the frequency range where the experimental response was stationary (Table 1). Spike-frequency–dependent adaptation ({alpha}) was indeed required by the models to fit the data points at the steady state (La Camera et al. 2002Go).

No difference was found in the response properties of the neurons to noisy currents characterized by an autocorrelation time length {tau}1 = 1 ms (n = 17) and 5 ms (n = 18). Thus, in both cases the responses could be captured by the white-noise input approximation, mentioned under METHODS.

However, collected data resembled the aLIF and not the aCLIFF responses (Fig. 3), as opposed to pyramidal cells in acute slices, where both models were reported to describe the data set equally well (Rauch et al. 2003Go). The aLIF model was substantially better in fitting the responses of the cells over the entire set of experiments (n = 35) (see Fig. 3 and Table 1).



View larger version (37K):
[in this window]
[in a new window]
 
FIG. 3. The plots report the results from 6 different experiments. Current-to-rate response of cultured neurons, dissociated from the rat neocortex, can be captured remarkably well by the integrate-and-fire (IF) dynamics. The responses were evaluated at the steady state as a function of m (markers), for increasing values of s (i.e., s1 < s2 < s3, estimated from the injected current realizations) and compared to the corresponding best-fit prediction (lines), provided by the aLIF (left panels) and the aCLIFF models (right panels) (see also Fig. 2B). The aLIF model reproduced the responses of the cultured neurons with high accuracy, whereas the aCLIFF model routinely failed in passing the fit-confidence test. In both cases of free search and penalty best-fit search, the aLIF showed superior performances and a lower K (see METHODS) (e.g., in the free-search best fit, K = 1.31 ± 0.09 for the aLIF; K = 2.03 ± 0.12 for the aCLIFF model). We note however that, even when the {chi}2-test gave low performances for both models, the average absolute discrepancy on the predicted frequency was always well below 2 Hz (i.e., 0.75 ± 0.07 Hz for the aLIF and 1.2 ± 0.11 Hz for the aCLIFF).

 
Moreover, cultured neurons showed very similar response properties compared to neurons from acute slices, when repetitively stimulated by the same noise realization. This was shown to yield a much higher precision in the timing of individual spikes, compared to a DC stimulus (Mainen and Sejnowski 1995Go), it was confirmed in a set of experiments [n = 5, 26–43 days in vitro (DIVs), not shown], and it suggests that the somatic spike emission mechanism of cultured neurons is not intrinsically noisy.

Best-fit effective parameters and passive membrane properties

As mentioned in METHODS, the passive membrane properties ({tau}m and Cm) of the patched neurons were routinely measured. As opposed to the previous report by Rauch et al. (2003)Go, the passive properties and the corresponding best-fit effective parameters of the IF models showed a significant cross-correlation, both for the best-fit free search and the penalty search. We report the results of the Pearson's r-test and, between parentheses, those of the Kendall's Tau test. In the best-fit free search, between Cm and C the correlation coefficient was 0.66 (0.54) for the aLIF and 0.73 for the aCLIFF model, whereas between {tau}m and {tau} it was 0.77 (0.57) for the aLIF. In the best-fit search with penalties, between Cm and C the coefficient was 0.56 (0.44) for the aLIF and 0.68 for the aCLIFF model, whereas between {tau}m and {tau} it was 0.81 (0.61) for the aLIF. Such results were validated by a very high significance of the nonparametric test (i.e., pK < 10–5; see METHODS).

Interspike interval variability and subthreshold voltage distribution

As stated in the previous sections, by using the set of best-fit effective parameters, the steady-state response of the aLIF model accurately matches the corresponding experimental mean firing rates, in the plane (m, s2). In layer V pyramidal neurons the best-fit parameters could sometimes account for the CV of the interspike intervals distribution as well (Rauch et al. 2003Go), although the fit criterion involved the mean firing rates only. Carrying out a similar analysis on the aLIF model in cultured cells, it turned out that model responses poorly matched the CV, experimentally estimated at the steady state (not shown).

However, when the steady-state subthreshold distribution of the membrane voltage, recorded in the experiments, was estimated and compared to the aLIF model prediction, a good agreement was observed over the available set of pairs (m, s2) (see Fig. 4). To quantitatively compare the aLIF model behavior to the experimental voltage traces, individual action potentials were clipped and a unique offset and a scaling factor were required for the model internal state-variable V to optimally match the voltage distribution (i.e., V -> {beta}VE'). For each cell, these additional 2 parameters do not affect the current-to-rate response function and they were the same over the pairs (m, s2) (Fig. 4).

Slow/cumulative inactivation and the stationary spike frequency

As indicated in METHODS, special care was taken in assessing the stationarity of the neuronal responses, with the aim of direct comparison with the model predictions, available at the steady state. Nevertheless the temporal dynamics of the output firing rate was characterized by fast (frequency-dependent) adaptation components, occurring over a time scale of several hundreds of milliseconds (see Fig. 1D), and by a slower component, occurring over a time scale of several seconds. Although the steady-state effect of both adaptation processes was captured by the model (i.e., by {alpha}; see METHODS), a cumulative (inactivation-related) component contributed to set an upper limit to the maximal stationary response frequency, sustained by the neuron over the entire duration of the stimulation (Fleidervish et al. 1996Go; Powers et al. 1999Go; Rauch et al. 2003Go; Sanchez-Vives et al. 2000Go; Sawczuk et al. 1997Go). For instance, the injection of a noisy stimulus current for 20 s with m > 200–300 pA or more, depending on the cell input resistance, resulted in a slowly decaying instantaneous firing rate, eventually turning into a cumulative inactivation of the action-potential generation (not shown). This was reminiscent of the slow cumulative sodium-current inactivation characterized by Fleidervish et al. (1996)Go, and made it impossible to quantify the steady-state responses above about 30 Hz, with the current stimulation protocol.

The presence of such a cumulative inactivation was explicitly tested (n = 8, 21–43 DIVs) by extending the protocol described in (Fleidervish et al. 1996Go; Schwindt et al. 1989Go), consisting in a repeated-pulse stimulation, lasting 1 s, with a very short recovery time. Under noisy current injection, by using the same current realization for each repetition, the same phenomenon occurred, although the voltage fluctuations, induced by the nondeterministic stimulus waveform, delayed the onset of the inactivation at parity of m, and sometimes transiently reversed the inactivation for a few tens of milliseconds, compared to DC stimuli. As expected, larger values of {tau}1 induced slow modulations on the membrane voltage trajectory (see Svirskis and Rinzel 2000Go), making the episodes of partial transient recovery more frequent, at parity of m and s (not shown).

The spontaneous emergence of patterned network activity in vitro

Networks of neocortical dissociated neurons show spontaneous collective patterned activity (Fig. 5A). Such an activity starts as asynchronous and spatially uncorrelated firing, toward the end of the first week in culture, and evolves into nonperiodic, synchronized, population bursting after 3–4 wk (Kamioka et al. 1996Go). Such an activity is reported to stay unchanged for more than 8 wk and thus represents the mature state of the network (Marom and Shahaf 2002Go). During such a period, single cells emit rare and irregular spikes or bursts of action potential, superimposed on spontaneous voltage fluctuations around a membrane potential of about –60 mV (Nakanishi and Kukita 1998Go). A similar spontaneous network-driven activity constitutes a typical feature of dissociated neuronal networks (Marom and Shahaf 2002Go) and was also observed in adult neocortex (Sanchez-Vives and McCormick 2000Go).

The recording of the network activity by means of MEAs led to a quantitative characterization of the distributions of the interburst interval (IBI) and of the population burst durations (PBd), in 7 cortical cultures. Under control extracellular medium, the IBIs were characterized by a mean ranging from 4.6 ± 0.4 to 30.3 ± 4.8 s, and by a CV ranging from 44 to 70%. The PBs were characterized by a mean ranging from 54 ± 8 to 146 ± 8.5 ms, and by a CV ranging from 26 to 80%. Under pharmacological disinhibition (see METHODS), the PBs became longer with a mean ranging from 600 ± 70 ms to 1.7 ± 0.053 s, and a CV ranging from 7 to 30%, whereas IBIs were characterized by a mean ranging from 12.7 ± 1.3 to 51.8 ± 11.1 s, and by a CV ranging from 31 to 96%. Figure 5D summarizes graphically the results from 4 experiments, under bicuculline.

Simulations of networks of aLIF model neurons

By using the single-neuron effective parameters, we computer-simulated the collective electrical activity, emerging from a small homogeneous population of excitatory aLIF neurons (see METHODS). We considered estimates for the synaptic connectivity available from the literature (i.e., Cee; see METHODS), and we set (m0, s0) in Eq. 4 to match the spontaneous low-rate asynchronous background activity observed in the MEAs experiments ({approx}1 Hz).

Similarly to previous theoretical reports (Segev et al. 2001Go; Tateno 2002Go; van Vreeswijk and Hansel 2001Go; Wiedemann and Lüthi 2003Go), but with a stronger original motivation for the use of the aLIF model, and a quantitative goal to compare computer simulations to the available MEAs recordings, we found that the simulated network activity consisted of asynchronous activity and of population bursts.

Although the individual neurons of the simulated network were not intrinsic burster cells and no pacemaker mechanisms had been introduced in the simulations, the emerging activity evolved into a collective network-driven bursting, depending on the strength of synaptic coupling Je (see Fig. 8). In agreement with the experiments of Maeda et al. (1995)Go, the spatial origin of PBs varied randomly with each burst, consistent with the conclusions about the lack of a unique pacemaker mechanism, driving the network.



View larger version (16K):
[in this window]
[in a new window]
 
FIG. 8. In the absence of adaptation, as extensively discussed in the literature, network bistability is possible in a homogeneous excitatory population, in a small range of synaptic coupling Je (A, {alpha} = 0), as reported by the phase diagram. This is a direct consequence of the profile of the network response function, as plotted in Fig. 6A. Under the hypotheses of a separation between time scales (see RESULTS), a similar phase diagram can be considered for the full model network, including adaptation (B). By studying the quasi-stationary/instantaneous stable states, it is possible to make quantitative predictions on the mean firing rates during asynchronous regimes as well as during population bursts. Insets: actual traces of the simulated network activity, for different synaptic coupling, and comparison of the actual resting and bursting activity levels to the predictions of the theory (black and gray dotted lines). In both panels, black continuous and dashed traces indicate stationary and quasi-stationary stable states, respectively. Dashed gray lines report the location of unstable stationary and quasi-stationary equilibria, whose distance from frest determines the regular/irregular character of network oscillations. Markers represent the network mean firing rates, measured in the computer simulations under different regimes. Although the simulated network considered was very small, Ne = 100, the agreement with the theory is remarkable. Parameters as in Fig. 5.

 
On a first approximation, the simulated IBIs had a lower bound set by the time scale of the recovery processes of adaptation {tau}a, which was set in the simulations in the range of 0.7–2 s (see Fig. 5, B and C, gray shading). Transiently and following each PB, the network activity was almost completely suppressed by the adaptation currents in each neuron and it later recovered. On the other hand, the synaptic coupling as well as the strength of the spontaneous synaptic release or other sources of randomness (i.e., m0 and s0) were correlated to the bursting frequency.

Three different global regimes were observed. The first one corresponds to a situation in which the excitatory synaptic interactions between neurons are very weak or the network connectivity is very low. Under such conditions, the network of model neurons was exclusively dominated by low-rate asynchronous activity ({approx}1 Hz). No PB occurred either spontaneously or evoked by any brief depolarizing stimulus (as opposed to Fig. 7A). Actually, in such a regime the global dynamics was dominated by a single low-rate stable state. These conditions approximate early developmental stages of cultured networks, described in the literature to display a very similar electrical activity pattern.



View larger version (33K):
[in this window]
[in a new window]
 
FIG. 7. The comparison between the numerical computer simulations and the predictions of the extended mean-field analysis shows a remarkable agreement, in the case of no adaptation (i.e., {alpha} = 0; parameters as in Fig. 5). A: In the top panels, the theoretical analysis of the stationary network activity has been carried out under 3 different conditions: m0 = 25.1 + {Delta}m0 pA, m0 = 25.1 pA, and m0 = 25.1 – {Delta}m0 pA and the stable (unstable) equilibrium states have been evaluated and graphically represented by filled circles (star), as done in Fig. 6A. The remaining traces report (from the top to the bottom) the temporal evolution of m0, the membrane voltage of 3 IF neurons taken by chance, the raster plot of the spikes emitted by a subset of the neurons, and the resulting population mean firing rate, overlapped to the predicted steady-state regimes (gray dotted lines). B: the same simulation was repeated after increasing the excitatory synaptic coupling between neurons, without any depolarizing external stimulation. Network finite-size effects may induce a spontaneous transition from one stable state to the other (A, Je = 10.75 pA; B, Je = 11.5 pA; scale bars; 200 ms and 25 Hz/50 mV).

 
The second regime corresponds to an increased strength of excitatory synaptic interactions between neurons. It has been shown that as a real network matures (from 3 to 40 days), the frequency and propagation velocity of the PBs markedly increases (Maeda et al. 1995Go), suggesting an increase in the number of synaptic contacts as well as in the synaptic efficacy. Therefore, this situation might be regarded as an evolution of the previous regime, determined by neurite outgrowth, synaptic development, and by homeostatic and activity-dependent plasticities (Desai et al. 1999Go; Turrigiano et al. 1998Go). While the low-rate asynchronous activity was still characterizing the spontaneous electrical activity of the network at the same frequency, the occurrence of spontaneous PBs was possible with a probability dependent on the average synaptic strength and, inversely, on the size of the network. The last dependency is determined by the random fluctuations in the population firing rate, which are related to finite-size effects (Mattia and Del Giudice 2002Go). The features of spontaneous bursting were also determined by the same factors, qualitatively resulting in a rare and unpredictable occurrence of short and irregular PBs, or in a frequent and more regular occurrence of longer and more regular PBs. Under such conditions, a brief external triggering stimulus could successfully recruit the neurons to transiently sustain an evoked PBs, whose duration is inversely determined by {alpha}. Actually, in such a regime the global dynamics of the network is transiently bistable and once a PB is started, the adaptation slowly redefines the location and the existence of the stable states for the network dynamics, until there is suddenly only a single stable state at 0 Hz. This somehow corresponds to a reset for the collective network activity. In details, the adaptation hyperpolarizing contribution, which starts to build up in every neuron recruited by the PB, decreases the output firing rate until the recurrent synaptic inputs to any neuron stop.

In such a regime, the statistics of the simulated IBIs and PBd matched the experiments, performed under pharmacological disinhibition (see Fig. 5, BD).

The third and last regime occurs for stronger excitatory synaptic interactions in the network. Under such conditions, the network is characterized by a high-rate asynchronous regime ({approx}55–80 Hz). In this regime, the firing of the individual neurons is more regular than in the previous ones. Moreover, any attempt at transiently silencing the activity of the network, by hyperpolarizing a large fraction of the neurons of the population, would not prevent the network to later recover its global state (as opposed to the simulations reported in Fig. 7). Although, under physiological conditions, it is probably not possible for a cortical network to sustain such a firing regime for a long time, an interesting phenomenon was observed for an intermediate synaptic interactions strength: occasional population breaks. This appears as rare and unpredictable simultaneous short interruptions of the global activity and is determined by the finite-size fluctuations of the activity.

The single-neuron response and the mean-field theory interpretation

In the present section, we show that the quantitative knowledge of the single-neuron response function, identified in the previous experiments, is very relevant to predict and interpret the emergence of the network activity described above. As widely discussed in the literature, a homogeneous network of synaptically interacting excitatory neurons may be regarded as a single dynamical system. Its stationary states can then be predicted and interpreted, in the limit of an infinite number of neurons Ne -> {infty}, by using the mean-field Eq. 4 and studying {Phi}aLIF [m(f), s(f)] as a function of f (Amit and Brunel 1997bGo).

We first consider the collective activity in the absence of spike-frequency–dependent adaptation (i.e., {alpha} = 0). Under such hypotheses, because of the simultaneous dependency of {Phi}aLIF on m and s2, the collective firing rate of the network may be characterized by 2 stable dynamic-equilibrium states, in a small range of average synaptic coupling Je (Figs. 6A and 8, left panel) (Amit and Brunel 1997bGo; Fusi and Mattia 1999Go). Such global activity configurations correspond to the solutions f* of the following self-consistent network equation, further satisfying a stability condition

These solutions have been graphically identified as the intersection points of {Phi}aLIF(f) with the unitary-slope line (see Figs. 6A and 7, insets).



View larger version (18K):
[in this window]
[in a new window]
 
FIG. 6. On the basis of the single-neuron response properties, it is possible to predict the collective activity emerging in a homogeneous population of excitatory neurons. The plot represents the profile of {Phi}aLIF as a function of f, (A) for increasing values of the average synaptic efficacy Je. In a small range of strong synaptic coupling Je, the network is characterized by two distinct dynamic equilibrium states for the collective activity ({alpha} = 0, circles: stable states; star: unstable state) (see Fig. 8A). B: {Phi}aLIF is plotted as a function of f, while comparing the instantaneous network response profile (see RESULTS) immediately before a PB (resting) and after the complete buildup of adaptation currents, at the end of the same PB (fully adapted) (see Fig. 8B). As discussed by many authors, adaptation reduces the gain of the network responsiveness curve. When the network is active at a low frequency (resting) regime frest, the impact of the spike-frequency–dependent adaptation is negligible and an additional dynamic equilibrium stable state, at a high frequency fburst, characterizes the collective dynamics. However, if a transition frest -> fburst occurs, the adaptation will slowly build up, bending the network response profile until no stable dynamic equilibrium state at high frequency can be sustained anymore [f(t) -> fH]. Although a similar graphical analysis does not provide direct information about the temporal evolution of the population activity (but see Fig. 8), it can be inferred that the actual network mean rate will be flipping between two states (i.e., ~50 and ~1 Hz) in an alternating and activity-dependent way, as confirmed by the computer simulations of Fig. 5. Parameters are as in Fig. 5.

 
In such dynamic equilibria, the population activity is the result of reverberating interactions between the neurons, which are sustained by excitatory synapses. The transitions from one state to the other can be triggered by an external brief stimulus affecting each neuron (i.e., m0 -> m0 ± {Delta}m0) and resulting into an intrinsically bistable network activity (Fig. 7A). This has been already described as a possible neuronal correlate of the selective delay-activity states in vivo (Yakovlev et al. 1998Go) and UP-states in vitro (Compte et al. 2003Go). Interestingly, such transitions can even occur spontaneously in small networks (Fig. 6B), as a result of the fluctuations induced by finite-size effects (Mattia and Del Giudice 2002Go).

The mechanisms of a PB

Considering the full network model, where individual neurons keep adapting their output rate as a function of the activity, it is possible to carry out a simple approximate analysis. Provided that the mechanisms responsible for the excitability reduction (e.g., the adaptation) act on a time scale ({tau}a) that is longer compared to the single-neuron dynamics ({tau}), an analysis of the quasi-stationary equilibria may tell us a lot about the collective activity of the network. In other words, by assuming that adaptation is delayed and transiently uncoupled from the neuronal dynamics, we may consider it as frozen and determined by the previous global regime.

Let us consider the situation of weak excitatory synaptic interactions between the model neurons: a low-rate stable regime frest characterizes the network dynamics because it can be immediately determined from {Phi}aLIF. By making explicit the dependency on m and s, as well as on {alpha} (see METHODS and Eqs. 2 and 4), we can write the following self-consistent equation and numerically find the stable solution frest

Such a solution is approximately located at 1 Hz and the instantaneous network firing rate f(t) fluctuates around it, as a result of finite-size network effects. However, every neuron of the network will instantaneously respond to the synaptic input drive according to

with frest a constant. In other words, the network is not immediately experiencing the impact of adaptation. Therefore, assuming that the network was at "rest," we may study and determine any additional stable dynamical attractor f*, which may exist even though transiently. This is done by numerically solving

The last equation is satisfied at least by f* = frest, and for increasing synaptic coupling, f* = fburst is also a solution (fburst > frest; see Figs. 6B and 8, right panel). Because of the profile of the single-neuron response function, under these circumstances there always exists an unstable equilibrium f'{theta} (frest < f'{theta} < fburst), separating the 2 basins of attraction (Fig. 6).

From the last considerations, it follows that when a fluctuations in the global activity is large enough to overcome the distance {Delta}'{theta} = (f'{theta}frest), the stability of frest may be (transiently) lost and the entire network synchronously shifts to a new regime where f(t) = fburst. Therefore, the unstable state acts as a kind of threshold or a no-return point, as the network dynamics is attracted toward a new stable state, reminiscent of what happens to the membrane voltage of a neuron during the generation of an action potential.

Anyway, such a regime cannot be sustained indefinitely because fburst is not a solution of the following self-consistent equation, which now includes the updated effects of adaptation

Instead, such an equation is satisfied by f* = fH {approx} 0 Hz, with fH < fburst (see Figs. 6B and 8, right panel). Therefore, while adaptation progressively builds up in individual neurons, the network activity decreases as the locations of the fburst and f'{theta} tend to become coincident, until their existence is lost (Fig. 6B). As soon as this happens, the steady-state network dynamics suddenly converges to fH and the neurons generally stop firing, reminiscent of the hyperpolarization experienced by the membrane voltage after an action potential. This accounts quantitatively for the generation of PBs, as tested by network simulations (Fig. 8, markers). The amount of time spent in the PB is therefore related to the distance {Delta}' = (fburst f'{theta}) and to the time requested by adaptation for the full buildup (i.e., {alpha} and {tau}a). Qualitatively, it can be concluded that for an increasing Je, {Delta}'{theta} decreases while {Delta}' increases; thus the mean IBI decreases and the PBd increases.

Similarly to the previous situation, such a new regime cannot be sustained indefinitely because fH is not a solution of

Because fH < frest, the network will recover the resting activity level frest, as described at the beginning.

Population breaks

As summarized in Fig. 8, adaptation changes the phase diagram of the network, creating a region of transient network-driven bistability. Interestingly, for stronger synaptic coupling, frest is no longer a stable solution of

Under such conditions, a high-rate equilibrium that we indicate again as fburst, dominates the network dynamics. Such a situation might be transiently disrupted by a fluctuation in the global activity, similarly to what has been discussed for frest. In fact, the numerical solutions of

include f* = fH and an unstable state f''{theta}, with fH < f''{theta} < fburst. As a consequence, if activity fluctuates enough, say to , with and , the entire network would transiently stop firing and synchronously shifts to a new low-rate regime that we call fH. Finally, such a regime cannot be sustained indefinitely, because fH is not a solution of

and the network will recover the activity level fburst, as described at the beginning. In analogy to the bursting events, the time spent in fH is related to the distance {Delta}" = (f''{theta}fH).

The variability of the IBIs distribution

From the previous considerations, it is possible to predict and interpret the regular/irregular character of the IBIs. In the last paragraphs, an analogy between PBs and action potentials occurring in an excitable membrane was proposed. In fact, similarly to the temporal evolution of the membrane potential in the single-neuron IF dynamics discussed previously (see METHODS), the activity of a network randomly fluctuates as a result of synaptic noise. Occasionally these fluctuations may be large enough to overcome an excitability threshold (i.e., {Delta}'{theta}, {Delta}''{theta}). When this happens, a major explosive event occurs (f -> fburst) and later the activity is strongly refractory to any further generation of PBs (f {approx} fH). Our proposal is to regard the generation of PBs as similar to the generation of an action potential, in a model of integration of a noisy input. Although the finite-size fluctuations are quite different from an Ornstein–Uhlenbeck process, because no drift is present and its variance increases with the mean network activity (Mattia and Del Giudice 2002Go) and because the network dynamics is more complicated than a linear integration of inputs (Gerstner 2000Go), a comparison to an IF-like model may be indeed proposed. By mapping the population mean firing rate f into the membrane voltage of an abstract LIF model neuron, we make the previous comparisons explicit, setting the resting membrane voltage to frest, the spike threshold to f'{theta}, the reset potential to fH, and the absolute refractory period proportional to {tau}a.

By increasing the network synaptic coupling Je, and thus decreasing the distance f'{theta}frest (see Fig. 8), the rate of threshold crossing is expected to monotonically increase, while preserving an irregular character. Such predictions were confirmed by the simulated network activity, where IBIs statistics are approximated very well by a Poisson process with a refractory time (Gerstner and Kistler 2002Go) (Fig. 5D, dotted thick line). This is reminiscent of a well-known balanced (i.e., drift-free) integration process (see Shadlen and Newsome 1998Go), where the threshold crossings are determined by subthreshold fluctuations only, in a noise-dominated regime.


 DISCUSSION
 
 TOP
 ABSTRACT
 INTRODUCTION
 METHODS
 RESULTS
 DISCUSSION
 GRANTS
 ACKNOWLEDGMENTS
 REFERENCES
 
In this work, we quantitatively evaluated the single-neuron discharge properties in dissociated cultures of neocortex. We proved that a simplified point-neuron IF model is an adequate description when the network mean firing rates are considered. Under noisy current input, the spike response properties of cultured neocortical neurons qualitatively resemble those of cells in acute slice preparations. This is of interest given that dissociated neurons undergo a different development of intrinsic biophysical properties, compared to the neurons in vivo.

Then, we presented and analyzed the collective activity arising in a simulated network of interacting excitatory model neurons in terms and on the basis of the network-response properties, emerging from f = {Phi}(m, s). The matching with the MEAs recording is satisfactory (Fig. 5, BD), indicating that the discharge response properties to noisy current stimuli and the experimentally characterized spike-frequency adaptation are sufficient to account for the emerging collective activity, observed in the experiments. Although the theoretical approach that we used holds in the limit of an infinite number of neurons Ne -> {infty}, the quantitative agreement between the approximate theoretical predictions and the numerical simulations is remarkable (Figs. 7 and 8), even with small-sized networks (e.g., Ne ~ 100) and incorporating the effect of adaptation. This further supports the validity of the mean-field hypotheses on the neuronal interactions in a recurrent network and the motivation for the noisy current-clamp protocol.

Cultured neurons respond as aLIF neurons at steady state

One of the results of the present work is that cultured neurons, dissociated from rat neocortex, show qualitatively similar response properties as those of pyramidal cells of acute brain slices, under noisy current clamp. Further supporting the observations of Rauch et al. (2003)Go, our results imply that the IF model is an appropriate description of the spike emission process under a realistic network input drive. As a consequence, when neurons are studied under a similar nondeterministic experimental paradigm, the morphological complexity, the biophysical intrinsic details, and the nonlinearities of the active ion currents may collapse to a reduced effective point-neuron dynamics (see also Rudolph and Destexhe 2003Go).

More data and experiments are indeed needed to interpret and understand why the predictive power of the IF model sometimes extend to the interspike intervals' higher-order statistics, such as the coefficient of variation, or to the subthreshold voltage distribution (Jolivet et al. 2004). However, under the assumption of the extended mean-field theory of a network of recurrently interacting neurons, higher-order features of the evoked spike trains and the subthreshold voltage dynamics do not affect the collective population activity, as reflected in the statistics of I(t) (see Eq. 4) (Amit and Brunel 1997bGo). Instead, these might play an important role in the context of synaptic plasticity (Buonomano and Merzenich 1998Go). Thus, our results suggest that the aLIF model, incorporating the effective parameters, may be adequate to simulate not only the collective activity but also long-term plasticities at the network level (Del Giudice and Mattia 2001Go; Fusi 2003Go).

The spontaneous emergence of patterned network activity

Among the global activity regimes, predicted by the theory and tested by numerical simulations, the first two regimes are the most physiological. Importantly, they further quantitatively reproduce the results obtained in our experiments (see Fig. 5D), as well as those corresponding to early developmental stages, in which the effective coupling between neurons is weak (Kamioka et al. 1996Go).

In our simulations, the synchronization across the neurons during PBs arises as a consequence of a global modulation of the population activity. This is a consequence of the alternating and opposing interplay of two mechanisms. The first is represented by the probability of spontaneous transitions, which occasionally recruit enough recurrent connections to trigger and sustain a persistent global spiking regime. The second is represented by the buildup of an activity-dependent mechanism, able to decrease the network excitability. The quasi-stationary analysis presented in the previous sections, to quantitatively interpret the collective activity of the network simulations, relies on the assumption that the refractory mechanisms responsible for an activity-dependent reduction of network excitability act on a time scale that is longer compared to that of the single-neuron dynamics. Because the experimental single-neuron response properties are mainly affected by spike-frequency adaptation, in our analysis we examined such an intrinsic neuronal refractory mechanism, but the proposed procedure might be extended to other (coexisting) mechanisms. However, the hypothesis of separation among dynamical time scales might not be always satisfied, and the consequences include the overestimation of the actual mean frequencies characterizing the PB regimes (i.e., fburst). This might play a minor role in the case of slow/cumulative inactivation, as mathematically described by Fleidervish et al. (1996)Go or by Giugliano et al. (2002)Go, or it could be more severe in the case of short-term synaptic depression (Tabak et al. 2000Go; Tsodyks et al. 2000Go; Wiedemann and Lüthi 2003Go), phenomenologically modeled following Tsodyks et al. (1998)Go because the related time constant changes as the presynaptic frequency changes.

The framework we discussed in our work may explain, in quantitative terms, the spontaneous periodic synchronized bursting activity (Fig. 5) observed in mature neuronal cultures and experimentally correlated to the in vitro synaptogenesis and development (Kamioka et al. 1996Go). Consistent with the in vitro connectivity pattern that is restricted to spatially neighboring sites, a cultured network could be thought as a homogeneous chain of weakly synaptically connected subpopulations of a few hundreds of neurons each, as described here. Each population would then have equal probability to start a PB, spreading to the entire culture by means of the longer range sparse excitatory connectivity, and characterize by the same properties of PBs occurring in an isolated subpopulation. Physical network sectioning experiments of Maeda et al. (1995)Go, and of Nakanishi and Kukita (1998)Go, provide evidence for such a global network architecture, consistent with the hypothesis that synchronized bursting is emerging from synaptic activity rather than by way of gap junctions and/or diffusible factors.

These considerations greatly expand the relevance of the experimental protocol we described and may contribute to a full understanding of the spontaneous dynamics emerging in in vitro preparations. As opposed to previous computational studies, in which extensive numerical computer simulations have been used, a unifying theory was considered. The same theory inspired the single-neuron experimental protocol and it allowed us to gain a deeper insight on the way a particular type of refractory mechanism might affect the global dynamics. Different kinds of refractoriness have a different impact, at the steady state and at quasi-steady states on the statistical properties of input currents: in terms of a very slowly changing hyperpolarizing membrane current, the spike-frequency adaptation affects the mean input current only, whereas synaptic depression induces a change to both the input mean and variance, through a change in Je (see Eq. 4). Intuitively, this would transiently deform the profile of the network response function (Fig. 6B), modifying the location of the intersections with the unitary-slope line, in a peculiar way.

Spontaneous activity and synaptic release are not just triggering such episodes; they also determine the configuration of the collective dynamics, as proven by the dependency on s2 of the single-cell discharge properties.

A limit of our study is that we did not take into explicit consideration the impact of NMDA glutamatergic receptors (Maeda et al. 1995Go). Although a similar mathematical approach could be possible (Brunel and Sergi 1998Go; Fourcaud and Brunel 2002Go), in the present work we ignored the longer time constant associated with the synaptic currents mediated by the NMDA receptors and we implicitly regarded it as contributing to the overall mean synaptic efficacies (i.e., Je). However, Wang (1999)Go already extensively showed that in a similar network architecture, synaptic components mediated by NMDA receptors support and greatly enhance network bistability.

Finally, a major limitation of the LIF model is constituted by the stereotyped description of the spike-initiation mechanisms. Real neurons have a nonlinear voltage–current relationship near threshold, where the activation of the fast sodium currents becomes relevant. This is not the case for the LIF, where the integration of the incoming synaptic inputs is always voltage independent. By comparing the response of biophysically realistic conductance-based and LIF model neurons to weakly modulated oscillatory input currents, it has been recently shown that the oversimplifications of the last affect the response properties to high-frequency input components (i.e., >50–100 Hz) and, in the same frequency range, they lead to an unphysiological dependency on the temporal correlations of the background synaptic noise (i.e., on {tau}I) (Fourcaud-Trocmé et al. 2003Go; H. Köndgen, C. Geisler, S. Fusi, X. J., Wang, H.-R. Lüscher, and M. Giugliano, unpublished observations).

These considerations are unlikely to qualitatively change the results described in the present work, given that the frequency of the PBs is generally well below 10 Hz in cultured networks. However, the prediction of the responsiveness of individual neurons to a sudden increase in their net synaptic input might be slightly affected, resulting in a wrong estimate of the slope of the rise and decay of each simulated PB. However, the work of Fourcaud-Trocmé et al. (2003)Go did not systematically take into consideration the contribution of adaptation, which may behave as a high-pass filter at low frequencies (Fuhrmann et al. 2002Go), thus compensating for the loss in responsiveness.

Finally, by extending the class of IF models to incorporate a nonlinear behavior at the threshold (Fourcaud-Trocmé et al. 2003Go), it is possible to recover the correct low- and high-frequency neuronal response properties as well as to reproduce the timing of individual spikes (Jolivet et al. 2004). Because the steady-state current-to-rate response function {Phi}(m, s2) can be still analytically determined for such extended IF models, an experimental identification of the model parameters may be devised, similar to what has been described here.

Network activity and single-neuron model details

An additional important issue to be discussed is represented by the quantitative link between model details and network activity. However, to determine a strong correlation between the measured single-neuron properties and the emerging network activity in vitro, the availability of additional detailed information about the precise neuronal and synaptic heterogeneity, coupling, topology, and density would be required. In fact, throughout our analysis, a few degrees of freedom still characterized the definition of the network model. These are represented by the synaptic coupling and by the extent of the neuronal ensembles participating in the population bursting. By making realistic assumptions and compensating for such a lack of experimental details by data reported in the literature (Nakanishi and Kukita 1998Go; Nakanishi et al. 1999Go), we could a posteriori confirm the consistency of some results in terms of the realistic range for the excitatory synaptic coupling Je in vitro. However, we further considered the systematic study of some aspects of the network activity, emerging in the computer simulations, as a function of the single-neuron parameters (see Fig. 9). This kind of analysis might improve our understanding of which features are indeed required for the network-patterned activation and it would be considerably difficult by using a more detailed neuronal model, characterized by several parameters and kinetic mechanisms.



View larger version (26K):
[in this window]
[in a new window]
 
FIG. 9. The analysis of Fig. 8B was repeated by keeping all the model parameters fixed (see Table 1), while changing {alpha} (AB), {tau}arp (CD), and H (EF), respectively. With the aim of outlining the distinct dependency of some features of the network activity on the single-neuron parameters, the most relevant quantities, {Delta}'{theta} and {Delta}''{theta} (left column), as well as the interburst firing rate during the PBs and for stronger synaptic coupling (right column), have been plotted as a function of the synaptic efficacy Je. Finally, by selecting the experimental data from a single cell and comparing the corresponding best-fit parameter sets for the aLIF and the aCLIFF model (GH), it is possible to predict that a network of aCLIFF neurons would perform very similarly to a network of aLIF units, apart from an offset in the synaptic efficacy Je and a slightly different sensitivity in {Delta}'{theta} and {Delta}''{theta}, and not in fburst and frest. Parameters used were: Ne = 100, Cee = 0.38, {tau}e = 5 ms, m0 = 25.1 pA, s0 = 92 pA; (aLIF): {alpha} = 3.18 pA·s, {tau} = 20.13 ms, C = 26.12 pF, H = 16.4 mV, {tau}arp = 47.7 ms; d and (aCLIFF) {alpha} = 3.0 pA·s, {lambda} = 4.2617 mV ms–1, C = 88.6 pF, H = 18.49 mV, {tau}arp = 48.21 ms.

 
From such an analysis, it appears that the most important parameter determining the emergence of PBs is {alpha} (Fig. 9, A and B). As already discussed, only when {alpha} > 0 a PB may occur (see Fig. 8). For increasing values of {alpha} (0; 15] pA·s, the range of Je for the emergence of irregular PBs increases (i.e., {Delta}'{theta} > 0), and the regimes of regular and frequent PBs extend to larger Je (Fig. 9A). The emergence of population breaks is affected by a change in {alpha} only in terms of an offset on the minimal synaptic coupling Je because the slope of {Delta}''{theta} remains unchanged (Fig. 9A). For a given synaptic coupling the amplitude of the population bursting (i.e., the intraburst spiking rate) fburst is insensitive to {alpha} (Fig. 9B), although the burst's duration decreases with {alpha} (not shown).

The effective parameters {tau}arp and H in the aLIF are related to the strength and kinetics of the outward membrane currents that are activated after each spike. By definition, the absolute refractory period {tau}arp is mainly limiting the maximal output firing rate of the single neurons of the network. From Fig. 9, C and D, it can be concluded that the regular/irregular character of the PBs (i.e., {Delta}'{theta}) is only weakly affected by increasing {tau}arp [0.5; 30] ms, which shrinks the range of Je where PBs may occur and increases the probability of population breaks (i.e., {Delta}''{theta} decreases). As expected, the amplitude of the PB strongly depends on {tau}arp (Fig. 9D).

The reset voltage H [0; 19] mV, following the emission of each action potential, does not substantially affect the amplitude and the probability of PB and population breaks, unless it approaches the excitability threshold {theta} (Fig. 9, E and F). Under these conditions, with the distance to cross the spiking threshold reduced in each neuron, the PBs occur for weaker synaptic coupling Je, whereas the probability of population breaks is considerably reduced.

Finally, it is important to underline that, provided that the essential features of the current–frequency response properties are retained, alternative mathematical models selected to quantitatively fit the same experimental response under noisy current clamp are expected to perform in a qualitatively similar way when their network activity is considered. We investigated whether this was indeed the case for the aCLIFF model, which routinely did not pass the fit tests. Instead of considering and comparing the aLIF and the aCLIFF characterized by average values for the parameters, we decided to look at the 2 parameter sets obtained for the same cell. The resulting simulated network activity is very similar (Fig. 9, G and H), with the exception of a different threshold for the emergence of the population bursts and a slightly reduced dynamic range for the synaptic coupling Je corresponding to irregular PBs. Therefore a network of aCLIFF, or even conductance-based model neurons matching the experimental current-to-firing rate response function and the dependency on the amplitude of the fluctuations s, would perform very similarly.

Summary of the results and predictions

With the aim of summarizing the results discussed so far and comparing our predictions to additional in vitro experiments, we report a list of the most important issues. Most of the suggested experimental protocols described below were indeed tested and the predicted resulting network transitions were confirmed by several investigators.

  1. Temporally patterned aperiodic synchronization of the electrical activity may emerge in a network, without ad hoc single-cell pacemaker mechanisms (Maeda et al. 1995Go) or ad hoc active membrane conductances.
  2. Because the predictions about the simulated collective dynamics are in good agreement with in vitro data, the protocol used for identifying single-neuron response properties was adequate to perform the reduction leading to IF models.
  3. As long as the network activity is of concern, on the basis of the extended mean-field hypotheses, only single-neuron firing rates are determining the in vitro collective activity.
  4. The glutamatergic synapses, the sporadic irregular firing of individual neurons, and the finite-size activity fluctuations are sufficient to account for the mechanisms of PBs initiation, which can be assimilated to a noise-dominated first-passage time with dead time.
  5. The mechanism of PBs termination is consistent with the temporal development of intrinsic neuronal adaptation (Robinson et al. 1993Go) and, together with the location of stable and unstable activity configurations (Fig. 8), determines the duration of the PBs and the statistics of the IBI distribution.
  6. Any pharmacological/ionic enhancement of synaptic transmission mediated by AMPA and NMDA receptors (e.g., by removing AMPA-R inactivation by cyclothiazide; see Chiappalone et al. 2003Go), decreasing (increasing) the extracellular concentration of Mg++ (Ca++), etc.), is expected to shift the network activity from asynchronous activity toward a more regular bursting, with longer-lasting PBs (Canepari et al. 1997Go; Nakanishi and Kukita 1998Go).
  7. Any pharmacological/ionic suppression of the AMPAr- and NMDAr-mediated synaptic transmission (e.g., by CNQX, D-APV, increasing the extracellular concentration of Mg++ or decreasing Ca++) is expected to shift the network activity back to asynchronous activity and irregular rare bursting, characterized by shorter PBs and a lower intraburst spike rate (Canepari et al. 1997Go; Kamioka et al. 1996Go; Nakanishi and Kukita 1998Go).
  8. Developmental synaptogenesis (Kamioka et al. 1996Go) as well as the use of a repeated electrical stimulation to induce activity-dependent long-term potentiation (LTP) of glutamatergic synapses, correspond to an increase in Je (Maeda et al. 1998Go). This is expected to change the spontaneous network activity from an asynchronous regime, or from rare PBs composed by a few spikes, to more frequent PBs composed by a higher number of intraburst spikes (see Fig. 8B). The same is expected to hold for electrically evoked PBs, and finally a graded modification of synaptic strength would correspond to a graded increase in PBs frequency and in its amplitude (i.e., the intraburst frequency), according to Fig. 8.
  9. A brief distributed electrical stimulation is expected to elicit a PB, even in a silent culture (Jimbo et al. 2000Go). The existence and location of stable network states could be unmasked in such a way.
  10. By evoking PBs of increasing amplitudes, an indirect evidence for the role of {Delta}' (see RESULTS) might be provided by an analysis of the evoked PB duration, expected to be longer for those stimuli that evoke more intraburst spikes (Maeda et al. 1998Go).
  11. An increase in the single-cell excitability, by raising the extracellular potassium concentration [K+], would correspond to a permanent upregulation of the network response function (i.e., like m0 -> m0 + {Delta}m0; see Figs. 6 and 7). This is expected to increase frest and to decrease {Delta}'{theta}, making the bursting network activity more regular and frequent (Canepari et al. 1997Go).

In this work we investigated the single-neuron response properties that are relevant for the quasi-stationary collective network activity in vitro. We provided a minimal quantitative description of the experimental data set, summarizing the current-to-firing rate response properties of cultured neocortical neurons by just 5 parameters. Although accurate in accounting for the mean firing rate output of a neuron, under realistic re-created network inputs, such a description is simple enough to provide analytical predictions on the spontaneous network activity, emerging in the numerical computer simulations of a population of interacting neurons.

Because such collective activity regimes match the activity experimentally observed in the cultures of dissociated neurons in vitro, our results attempt at answering a more general question: among the richness of biological details, which features should a model of neuronal excitability retain, when the electrical activity, emerging from the synaptic interactions of individual cells, is of concern?

Under such a perspective, we are convinced that many other recent modeling studies, using simplified model neurons (Segev et al. 2001Go; Tal et al. 1998Go; Tateno 2002Go; Tsodyks et al. 1998Go, 2000Go; Wiedemann and Lüthi 2003Go), will greatly benefit from the availability of the effective model parameters, identified in our work, as well as from the observations of intrinsic neuronal properties such as the spike-frequency adaptation and the spike-emission slow cumulative inactivation, to perform large-scale computer simulation with increased realism.

Finally, it is interesting to note that the reduction driven by the experimental data at the single-neuron level did not compromise the richness of phenomena occurring at the level of a population. In particular, the network activity investigated here depends only quantitatively on the model details, provided that the experimental current-to-firing rate response function is reproduced and that the sensitivity to the input fluctuations retained.


 GRANTS
 
 TOP
 ABSTRACT
 INTRODUCTION
 METHODS
 RESULTS
 DISCUSSION
 GRANTS
 ACKNOWLEDGMENTS
 REFERENCES
 
This work was supported by the Swiss National Science Foundation Grants 31-61335.00 to H.-R. Lüscher and 3100-067048.01 to J. Streit and the Human Frontier Science Program Grant LT00561/2001-B to M. Giugliano. M. Arsiero was supported by the "EC Thematic Network in Neuroinformatics and Computational Neuroscience" and by the University of Genoa, Italy.


 ACKNOWLEDGMENTS
 
 TOP
 ABSTRACT
 INTRODUCTION
 METHODS
 RESULTS
 DISCUSSION
 GRANTS
 ACKNOWLEDGMENTS
 REFERENCES
 
We are indebted to S. Fusi, G. La Camera, and A. Rauch, whose previous work motivated the present research and we thank them for constant help. We are grateful to R. Rubli for assistance with neocortical dissociated cultures, to A. Tscherter, H. Köndgen, T. Berger, W. Senn, and V. Crunelli for helpful remarks on the manuscript, and to D. Benzi for discussions. M.G. is particularly grateful to M. Grattarola for having introduced him to the MEAs technique and to the field of computational neuroscience.

Present address of M. Arsiero: Institute of Physiology, University of Bern, Bühlplatz 5, 3012 Bern, Switzerland.


 FOOTNOTES
 
The costs of publication of this article were defrayed in part by the payment of page charges. The article must therefore be hereby marked "advertisement" in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.

Address for reprint requests and other correspondence: M. Giugliano, Institute of Physiology, University of Bern, Bühlplatz 5, 3012 Bern, Switzerland (E-mail: giugliano{at}pyl.unibe.ch).


 REFERENCES
 
 TOP
 ABSTRACT
 INTRODUCTION
 METHODS
 RESULTS
 DISCUSSION
 GRANTS
 ACKNOWLEDGMENTS
 REFERENCES
 
Abbott LF and Dayan P. Theoretical Neuroscience. Cambridge, MA: MIT Press, 2001.

Amit D and Brunel N. Dynamics of a recurrent network of spiking neurons before and following learning. Network Comput Neural Syst 8: 373–404, 1997a.[CrossRef]

Amit DJ and Brunel N. Model of global spontaneous activity and local structured (learned) delay activity during delay. Cereb Cortex 7: 237–252, 1997b.[Abstract/Free Full Text]

Amit DJ and Tsodyks MV. Effective neurons and attractor neural networks in cortical environment. Network 3: 121–137, 1992.[Web of Science]

Bernander O, Douglas RJ, Martin KA, and Koch C. Synaptic background activity influences spatiotemporal integration in single pyramidal cells. Proc Natl Acad Sci USA 88: 11569–11573, 1991.[Abstract/Free Full Text]

Brunel N. Persistent activity and the single cell f–I curve in a cortical network model. Network 11: 261–280, 2000.[Web of Science][Medline]

Brunel N and Sergi S. Firing frequency of leaky integrate-and-fire neurons with synaptic currents dynamic. J Theor Biol 195: 87–95, 1998.[CrossRef][Web of Science][Medline]

Brunel N and Wang XJ. What determines the frequency of fast network oscillations with irregular neural discharges? I. Synaptic dynamics and excitation–inhibition balance. J Neurophysiol 90: 415–430, 2003.[Abstract/Free Full Text]

Buonomano DV and Merzenich MM. Cortical plasticity: from synapses to maps. Annu Rev Neurosci 21: 149–186, 1998.[CrossRef][Web of Science][Medline]

Canepari M, Bove M, Maeda E, Cappello M, and Kawana A. Experimental analysis of neuronal dynamics in cultured cortical networks and transitions between different patterns of activity. Biol Cybern 77: 153–162, 1997.[CrossRef][Web of Science][Medline]

Chance FS, Abbott LF, and Reyes AD. Gain modulation from background synaptic input. Neuron 35: 773–782, 2002.[CrossRef][Web of Science][Medline]

Chiappalone M, Vato A, Tedesco MT, Marcoli M, Davide FA, and Martinoia S. Networks of neurons coupled to microelectrode arrays: a neuronal sensory system for pharmacological applications. Biosens Bioelectron 18: 627–634, 2003.[CrossRef][Web of Science][Medline]

Chicca E, Badoni D, Dante V, D'Andreagiovanni M, Salina G, Carota L, Fusi S, and Del Giudice P. A VLSI recurrent network of integrate-and-fire neurons connected by plastic synapses with long term memory. IEEE Trans Neural Networks (Special Issue on Hardware Implementation) 14: 1297–1307, 2003.[CrossRef]

Compte A, Sanchez-Vives MV, McCormick DA, and Wang X-J. Cellular and network mechanisms of slow oscillatory activity (1 Hz) and wave propagations in a cortical network model. J Neurophysiol 89: 2707–2725, 2003.[Abstract/Free Full Text]

Cox DR and Miller HD. The Theory of Stochastic Processes. London: Chapman & Hall, 1965.

Del Giudice P and Mattia M. Long- and short-term synaptic plasticity and the formation of working memory: a case study. Neurocomputing 38: 1175–1180, 2001.[CrossRef][Web of Science]

Desai NS, Rutherford LC, and Turrigiano GG. Plasticity in the instrinsic excitability of cortical pyramidal neurons. Nat Neurosci 2: 515–520, 1999.[CrossRef][Web of Science][Medline]

Destexhe A, Mainen ZF, and Sejnowski TJ. Synaptic transmission and neuromodulation using a common kinetic formalism. J Comput Neurosci 1: 195–230, 1994.[CrossRef][Medline]

Destexhe A and Paré D. Impact of network activity on the integrative properties of neocortical pyramidal neurons in vivo. J Neurophysiol 81: 1531–1547, 1999.[Abstract/Free Full Text]

Destexhe A and Paré D. A combined computational and intracellular study of correlated synaptic bombardment in neocortical pyramidal neurons in vivo. Neurocomputing 32: 113–119, 2000.[CrossRef]

Destexhe A, Rudolph M, and Paré D. The high-conductance state of neocortical neurons in vivo. Nat Rev 4: 739–751, 2003.[CrossRef]

Fleidervish I, Friedman A, and Gutnick MJ. Slow inactivation of Na+ current and slow cumulative spike adaptation in mouse and guinea-pig neocortical neurons in slices. J Physiol 493: 83–97, 1996.

Fourcaud N and Brunel N. Dynamics of firing probability of noisy integrate-and-fire neurons. Neural Comput 14: 2057–2110, 2002.[CrossRef][Web of Science][Medline]

Fourcaud-Trocmé N, Hansel D, van Vreeswijk C, and Brunel N. How spike generation mechanisms determine the neuronal response to fluctuating inputs. J Neurosci 23: 11628–11640, 2003.[Abstract/Free Full Text]

Froemke RC and Dan Y. Spike-timing-dependent synaptic modification induced by natural spike trains. Nature 416: 433–438, 2002.[CrossRef][Medline]

Fuhrmann G, Markham H, and Tsodyks M. Spike frequency adaptation and neocortical rhythms. J Neurophysiol 88: 761–770, 2002.[Abstract/Free Full Text]

Fusi S. Spike-driven synaptic plasticity for learning correlated patterns of mean firing rates. Rev Neurosci 14: 73–84, 2003.[Web of Science][Medline]

Fusi S and Mattia M. Collective behavior of networks with linear (VLSI) integrate and fire neurons. Neural Comput 11: 633–652, 1999.[CrossRef][Web of Science][Medline]

Gerstner W. Population dynamics of spiking neurons: fast transients, asynchronous states, and locking. Neural Comput 12: 43–89, 2000.[CrossRef][Web of Science][Medline]

Gerstner W and Kistler W. Spiking Neuron Models. Cambridge, UK: Cambridge Univ. Press, 2002.

Giugliano M and La Camera G, Rauch A, Lüscher H-R, and Fusi S. Non-monotonic current-to-rate response function in a novel integrate-and-fire model neuron. Proc ICANN2002, edited by Dorronsoro JR. New York: Springer-Verlag, 2002, p. 141–146.

Giugliano M and Lüscher H-R. Irregular in vivo–like background synaptic activity recreated in in vitro neocortial slices. Proc Sixth IBRO World Conf of Neurosci, Prague, 2003.

Giugliano M, Lüscher H-R, and Berger T. Reverberating network activity induced by distributed electrical stimulation in neocortical slices. Abstract Viewer/Itinerary Planner. Washington, DC: Soc. for Neurosci. Program 173.7, 2003.

Hamill OP, Marty A, Neher E, Sakmann B, and Sigworth FJ. Improved patch-clamp techniques for high-resolution current recording from cells and cell-free membrane patches. Pfluegers Arch 391: 85–100, 1981.[CrossRef][Web of Science][Medline]

Iansek R and Redman SJ. An analysis of the cable properties of spinal motoneurones using a brief intracellular current pulse. J Physiol 234: 613–636, 1973.[Abstract/Free Full Text]

Jimbo Y, Kawana A, Parodi P, and Torre V. The dynamics of a neuronal culture of dissociated cortical neurons of neonatal rats. Biol Cybern 83: 1–20, 2000.[CrossRef][Web of Science][Medline]

Jimbo Y, Tateno T, and Robinson HPC. Simultaneous induction of pathway-specific potentiation and depression in networks of cortical neurons. Biophys J 76: 670–678, 1999.[Web of Science][Medline]

Jolivet R, Lewis TJ, and Gerstner W. Generalized integrate-and-fire models of neuronal activity approximate spike trains of a detailed model to a high degree of accuracy. J Neurophysiol 92: 959–976, 2004.[Abstract/Free Full Text]

Kamioka H, Maeda E, Jimbo Y, Robinson HPC, and Kawana A. Spontaneous periodic synchronized bursting during formation of mature patterns of connections in cortical cultures. Neurosci Lett 206: 109–112, 1996.[CrossRef][Web of Science][Medline]

La Camera G, Rauch A, Senn W, Lüscher H-R, and Fusi S. Firing rate adaptation without losing sensitivity to fluctuations. Proc ICANN2002, edited by Dorronsoro JR. New York: Springer-Verlag, 2002, p. 180–185.

La Camera G, Senn W, and Fusi S. Equivalent networks of conductance- and current-driven neurons. Proc ICANN/ICONIP 2003.

Liu Y-H and Wang X-J. Spike-frequency adaptation of a generalized leaky integrate-and-fire model neuron. J Comput Neurosci 10: 25–45, 2001.[CrossRef][Web of Science][Medline]

Maeda E, Kuroda Y, Robinson HPC, and Kawana A. Modification of parallel activity elicited by propagating bursts in developing networks of rat neocortical neurones. Eur J Neurosci 10: 488–496, 1998.[CrossRef][Web of Science][Medline]

Maeda E, Robinson HP, and Kawana A. The mechanisms of generation and propagation of synchronized bursting in developing networks of cortical neurons. J Neurosci 15: 6834–6845, 1995.[Abstract/Free Full Text]

Mainen ZF and Sejnowski T. Reliability of spike timing in neocortical neurons. Science 268: 1503, 1995.[Abstract/Free Full Text]

Marom S and Shahaf G. Development, learning and memory in large random networks of cortical neurons: lessons beyond anatomy. Q Rev Biophys 35: 63–87, 2002.[CrossRef][Web of Science][Medline]

Mattia M and Del Giudice P. Efficient event-driven simulation of large networks of spiking neurons and dynamical synapses. Neural Comput 12: 2305–2329, 2000.[CrossRef][Web of Science][Medline]

Mattia M and Del Giudice P. Population dynamics of interacting spiking neurons. Phys Rev E 66: 051917, 2002.[CrossRef]

Mongillo G and Amit DJ. Oscillations and irregular emission in networks of linear spiking neurons. J Comput Neurosci 11: 249–261, 2001.[CrossRef][Web of Science][Medline]

Nakanishi K and Kukita F. Functional synapses in synchronized bursting of neocortical neurons in culture. Brain Res 795: 137–146, 1998.[CrossRef][Web of Science][Medline]

Nakanishi K, Kukita F, Asai K, and Kato T. Recurrent subthreshold electrical activities of rat neocortical neurons progress during long-term culture. Neurosci Lett 304: 85–88, 2001.[CrossRef][Web of Science][Medline]

Nakanishi K, Nakanishi M, and Kukita F. Dual intracellular recording of neocortical neurons in a neuron-glia co-culture system. Brain Res Protocols 4: 105–114, 1999.[CrossRef][Medline]

Paninski L, Lau B, and Reyes A. Noise-driven adaptation: in vitro and mathematical analysis. Neurocomputing 52: 877–883, 2003.[CrossRef]

Poliakov AV, Powers RK, and Binder MD. Functional identification of the input–output transforms of motoneurones in the rat and cat. J Physiol 504: 401–424, 1997.[Abstract/Free Full Text]

Potter SM. Distributed processing in cultured neuronal networks. In: Progress in Brain Research. Advances in Neural Population Cording, edited by Nicolelis MAL. Amsterdam: Elsevier, 2001, vol. 130, p. 49–62.

Powers RK, Sawczuk A, Musick JR, and Binder MD. Multiple mechanisms of spike-frequency adaptation in motoneurones. J Physiol (Paris) 93: 101–114, 1999.

Press W, Teukolsky SA, Vetterling WT, and Flannery BP. Numerical Recipes in C: The Art of Scientific Computing. Cambridge, UK: Cambridge Univ. Press, 1992.

Protopapas AD and Bower JM. Spike coding in pyramidal cells of the piriform cortex of rat. J Neurophysiol 86: 1504–1510, 2001.[Abstract/Free Full Text]

Rauch A, La Camera G, Lüscher H-R, Senn W, and Fusi S. Neocortical pyramidal cells respond as integrate-and-fire neurons to in vivo–like input currents. J Neurophysiol 90: 1598–1612, 2003.[Abstract/Free Full Text]

Renart A, Brunel N, and Wang X-J. Mean-field theory of recurrent cortical networks: from irregularly spiking neurons to working memory. In: Computational Neuroscience: A Comprehensive Approach, edited by Feng J. Boca Raton, FL: CRC Press, 2003, chap. 15, p. 431–490.

Reutimann J, Giugliano M, and Fusi S. Event-driven simulation of spiking neurons with stochastic dynamics. Neural Comput 15: 811–830, 2003.[CrossRef][Web of Science][Medline]

Robinson HPC, Kawahara M, Jimbo, Y, Torimitsu K, Kuroda Y, and Kawana A. Periodic synchronized bursting and intracellular calcium transients elicited by low magnesium in cultured cortical neurons. J Neurophysiol 70: 1606–1616, 1993.[Abstract/Free Full Text]

Rudolph M and Destexhe A. A fast-conducting, stochastic integrative mode for neocortical neurons in vivo. J Neurosci 23: 2466–2476, 2003.

Salinas E. Background synaptic activity as a switch between dynamical states in a network. Neural Comput 15: 1439–1475, 2003.[CrossRef][Web of Science][Medline]

Sanchez-Vives MV and McCormick DA. Cellular and network mechanisms of rhythmic recurrent activity in neocortex. Nat Neurosci 3: 1027–1034, 2000.[CrossRef][Web of Science][Medline]

Sanchez-Vives MV, Nowak LG, and McCormick DA. Cellular mechanisms of long-lasting adaptation in visual cortical neurons in vitro. J Neurosci 20: 4286–4299, 2000.

Sawczuk A, Powers RK, and Binder MD. Contribution of outward currents to spike frequency adaptation in hypoglossal motoneurons of the rat. J Physiol 78: 2246–2253, 1997.

Schwindt PC, Spain WJ, and Crill WE. Long-lasting reduction of excitability by a sodium-dependent potassium current in cat neocortical neurons. J Neurophysiol 61: 233–244, 1989.[Abstract/Free Full Text]

Segev R, Shapira Y, Benveniste M, and Ben-Jacob E. Observations and modeling of synchronized bursting in two-dimensional neural networks. Phys Rev E 64: 011920, 2001.[CrossRef]

Shadlen MN and Newsome WT. The variable discharge of cortical neurons: implications for connectivity, computation, and information coding. J Neurosci 18: 3870–3896, 1998.[Abstract/Free Full Text]

Shu Y, Hasenstaub A, and McCormick DA. Turning on and off recurrent balanced cortical activity. Nature 423: 288–293, 2003.[CrossRef][Medline]

Sjöström PJ, Turrigiano GG, an Nelson SB. Rate, timing, and cooperativity jointly determine cortical synaptic plasticity. Neuron 32: 1149–1164, 2001.[CrossRef][Web of Science][Medline]

Steriade M. Impact of network activities on neuronal properties in corticothalamic systems. J Neurophysiol 86: 1–39, 2001.[Abstract/Free Full Text]

Streit J, Tscherter A, Heuschkel MO, and Renaud P. The generation of rhythmic activity in dissociated cultures of rat spinal cord. Eur J Neurosci 14: 191–202, 2001.[CrossRef][Web of Science][Medline]

Svirskis G and Rinzel J. Influence of temporal correlation of synaptic input on the rate and variability of firing in neurons. Biophys J 5: 629–637, 2000.

Tabak J, Senn W, O'Donovan J, and Rinzel J. Modeling of spontaneous activity in developing spinal cord using activity-dependent depression in an excitatory network. J Neurosci 20: 3041–3056, 2000.[Abstract/Free Full Text]

Tal D, Jacobson E, Lyakhov V, and Marom S. Frequency tuning of input–output relation in a rat cortical neuron in-vitro. Neurosci Lett 300: 1–4, 1998.[CrossRef]

Tateno T. Noise-induced effects on period-doubling bifurcation for integrate-and-fire oscillators. Phys Rev E 65: 021901, 2002.[CrossRef]

Tscherter A, Heuschkel MO, Renaud P, and Streit J. Spatiotemporal characterization of rhythmic activity in spinal cord slice cultures. Eur J Neurosci 14: 179–190, 2001.[CrossRef][Web of Science][Medline]

Tsodyks M, Pawelzik K, and Markram H. Neural networks with dynamic synapses. Neural Comput 10: 821–835, 1998.[CrossRef][Web of Science][Medline]

Tsodyks M, Uziel A, and Markram H. Synchrony generation in recurrent networks with frequency-dependent synapses. J Neurosci 20: RC50, 1–5, 2000.[Abstract/Free Full Text]

Tuckwell HC. Introduction to Theoretical Neurobiology. Cambridge, UK: Cambridge Univ. Press, 1988.

Turrigiano GG, Leslie KR, Desai NS, Rutherford LC, and Nelson SB. Activity-dependent scaling of quantal amplitude in neocortical neurons. Nature 391: 892–896, 1998.[CrossRef][Medline]

van Vreeswijk C and Hansel D. Paterns of synchrony in neural networks with spike adaptation. Neural Comput 13: 959–992, 2001.[CrossRef][Web of Science][Medline]

Wang X-J. Synaptic basis of cortical persistent activity: the importance of NMDA receptors to working memory. J Neurosci 19: 9587–9603, 1999.[Abstract/Free Full Text]

Wang X-J. Synaptic reverberation underlying mnemonic persistent activity. Trends Neurosci 24: 455–463, 2001.[CrossRef][Web of Science][Medline]

Wiedemann UA and Lüthi A. Timing of network synchronization by refractory mechanisms. J Neurophysiol 90: 3902–3911, 2003.[Abstract/Free Full Text]

Yakovlev V, Fusi S, Berman E, and Zohary E. Inter-trial neuronal activity in infero-temporal cortex: a putative vehicle to generate long term associations. Nat Neurosci 1: 310–317, 1998.[CrossRef][Web of Science][Medline]




This article has been cited by other articles:


Home page
Cereb CortexHome page
H. Kondgen, C. Geisler, S. Fusi, X.-J. Wang, H.-R. Luscher, and M. Giugliano
The Dynamical Response Properties of Neocortical Neurons to Temporally Modulated Noisy Inputs In Vitro
Cereb Cortex, September 1, 2008; 18(9): 2086 - 2097.
[Abstract] [Full Text] [PDF]


Home page
J. Neurophysiol.Home page
O. Feinerman, M. Segal, and E. Moses
Identification and Dynamics of Spontaneous Burst Initiation Zones in Unidimensional Neuronal Cultures
J Neurophysiol, April 1, 2007; 97(4): 2937 - 2948.
[Abstract] [Full Text] [PDF]


Home page
J. Neurosci.Home page
M. Arsiero, H.-R. Luscher, B. N. Lundstrom, and M. Giugliano
The Impact of Input Fluctuations on the Frequency-Current Relationships of Layer 5 Pyramidal Neurons in the Rat Medial Prefrontal Cortex
J. Neurosci., March 21, 2007; 27(12): 3274 - 3284.
[Abstract] [Full Text] [PDF]


Home page
J. Neurophysiol.Home page
G. La Camera, A. Rauch, D. Thurbon, H.-R. Luscher, W. Senn, and S. Fusi
Multiple Time Scales of Temporal Response in Pyramidal and Fast Spiking Cortical Neurons
J Neurophysiol, December 1, 2006; 96(6): 3448 - 3464.
[Abstract] [Full Text] [PDF]


Home page
J. Neurophysiol.Home page
D. Golomb, A. Shedmi, R. Curtu, and G. B. Ermentrout
Persistent Synchronized Bursting Activity in Cortical Tissues With Low Magnesium Concentration: A Modeling Study
J Neurophysiol, February 1, 2006; 95(2): 1049 - 1067.
[Abstract] [Full Text] [PDF]


Home page
J. Physiol.Home page
F. J. L Arnold, F. Hofmann, C. P. Bengtson, M. Wittmann, P. Vanhoutte, and H. Bading
Microelectrode array recordings of cultured hippocampal networks reveal a simple model for transcription and protein synthesis-dependent plasticity
J. Physiol., April 1, 2005; 564(1): 3 - 19.
[Abstract] [Full Text] [PDF]


This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow All Versions of this Article:
92/2/977    most recent
00067.2004v1
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Right arrow Citation Map
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in Web of Science
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Web of Science (24)
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Giugliano, M.
Right arrow Articles by Streit, J.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Giugliano, M.
Right arrow Articles by Streit, J.


HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
Visit Other APS Journals Online
Copyright © 2004 by the The American Physiological Society.