Home | About | Journals | Submit | Contact Us | Français |

**|**PLoS One**|**v.5(11); 2010**|**PMC2975677

Formats

Article sections

Authors

Related links

PLoS One. 2010; 5(11): e13651.

Published online 2010 November 8. doi: 10.1371/journal.pone.0013651

PMCID: PMC2975677

Aldo Rustichini, Editor^{}

University of Minnesota, United States of America

* E-mail: ac.awattou@saijemj

Conceived and designed the experiments: HJK JT. Performed the experiments: JFM. Analyzed the data: JFM. Contributed reagents/materials/analysis tools: JFM HJK JT. Wrote the paper: JFM HJK JT.

Received 2010 August 31; Accepted 2010 October 2.

Copyright Mejias et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

This article has been cited by other articles in PMC.

Complex coherent dynamics is present in a wide variety of neural systems. A typical example is the voltage transitions between up and down states observed in cortical areas in the brain. In this work, we study this phenomenon via a biologically motivated stochastic model of up and down transitions. The model is constituted by a simple bistable rate dynamics, where the synaptic current is modulated by short-term synaptic processes which introduce stochasticity and temporal correlations. A complete analysis of our model, both with mean-field approaches and numerical simulations, shows the appearance of complex transitions between high (up) and low (down) neural activity states, driven by the synaptic noise, with permanence times in the up state distributed according to a power-law. We show that the experimentally observed large fluctuation in up and down permanence times can be explained as the result of sufficiently noisy dynamical synapses with sufficiently large recovery times. Static synapses cannot account for this behavior, nor can dynamical synapses in the absence of noise.

Neural systems, even in the absence of external stimuli, can exhibit a wide variety of coherent collective behaviors, as *in vivo* and *in vitro* experiments show [1]–[3]. One of the most prominent examples is the spontaneous transition between two different voltage states, namely up and down states, observed in simultaneous individual single neuron recordings as well as in local field measures. Such behavior, which is generated within the cortex, may provide a framework for neural computations [4], and could also coordinate some sleep rhythms into a coherent rhythmic sequence of recurring cortical and thalamocortical activities [3], [5], [6]. The phenomenon of *up and down transitions* has been measured in a number of situations, such as in the primary visual cortex of anesthetized animals [7], [8], during slow-wave sleep [1], [5], [6], in the somatosensory cortex of awake animals [9], or in slice preparation under different experimental protocols [3], [10], [11], to name a few. The origin of such structured neuronal activity is still unclear, although several studies have shown that both intrinsic cell properties [12]–[14] and the high level of recurrency present in actual neural circuits [3], [15], [16] may contribute to the generation of up and down transitions. In particular, the contribution that reverberations in recurrent neural networks may have in the appearance of these transitions could depend strongly on synaptic properties. It is known, for instance, that excitatory synapses with a slow dynamics (such as synapses mediated by NMDA receptors) may play a relevant role in the generation of persistent activity or up cortical states [17]. On the other hand, several modeling studies indicate that activity-dependent synaptic mechanisms, such as short-term synaptic depression and facilitation, can induce voltage transitions between up and down neural states as well [16], [18]–[20].

Many crucial points about the understanding of up and down transitions are, however, still lacking. For instance, *in vivo* experiments in the cat visual cortex show that the permanence times in the depolarized or up state present a high variability, and can range from a scale of milliseconds to seconds [7]. A similar level of irregularity has also been recently found in *in vivo* recordings of up-down transitions in the rat auditory cortex [21], as well as in sleep-wake transitions [15], [22], [23], where power-law distributions in the duration of wake states have been measured. Such complexity in the time series of the neuron membrane potentials remains far to be explained, and could reflect scale invariance in permanence times, which could in turn be a (preliminary) indicative of criticality. In fact, there are many recent studies that have shown criticality in different contexts in the brain [24], [25], as well as in neural network models which present self-organization and criticality properties [26]–[28], and even it has been reported to occur in sleep-wake transitions in *in vivo* conditions [22], [23]. Although it is worth noting that the irregularity of the dynamics of up and down states is not a sufficient condition for criticality, a concrete characterization of such irregularity may be a convenient starting point for future works on this topic.

To study in detail the relevant issue of irregular up and down cortical dynamics, we propose in this work a minimal model for up and down transitions in neural media. We consider a simple bistable rate model whose stable solutions represent two possible voltage states of the mean membrane potential of the network. More precisely, such states correspond, respectively, to high and low levels of activity in the network (that is, the up and down cortical states). In addition, we consider that the synaptic connections between neurons of the network present short-term synaptic depression (STD) mechanisms, which introduce temporal correlations, as well as synaptic stochasticity, in the dynamics of the system [29]–[32]. A complete analysis of this simple mathematical model depicts, both numerically and within a theoretical probabilistic approach, the appearance of power-law dependences in the distribution of permanence times in the up state. Our results show that the appearance of such scale free distributions is due to the complex interplay between several factors including synaptic stochasticity and the temporal correlations introduced by STD. The emergence of power-law dependences could, indeed, explain the high variability in permanence times in the up state suggested by experiments [7], [21].

Our starting point is a bistable rate model, which mimics the dynamics of the electrical activity of a population of interconnected excitatory neurons (although it can be easily extended to other situations) with two stable levels of activity. The model has the form [33]

(1)

where is the mean firing rate of the (homogeneous) neural population, is the maximum level of activity which can be reached by the population (in absence of noise), is the synaptic coupling strength in absence of STD, and is the firing threshold of the neurons in the population. The variable is a Gaussian white noise of zero mean and standard deviation , which takes into account the inner stochasticity of the neural population (caused by other sources of uncontrolled noise in the system). The parameter is the population time constant, which may be assumed to be around the duration of the synaptic current pulse [34], [35]. For generality purposes, we set , and therefore time and frequency are given in units of and , respectively. The term represents the transduction function, which gives the nonlinear effect that the mean postsynaptic current (coming from recurrent connections of the neural population) induces in the network mean firing rate. Employing this form for , the up and down stable levels of activity correspond to and , respectively.

On the other hand, the variable in equation (1) takes into account the dynamical modification of the strength of the synaptic connections during short time scales due to high network activity, and it is usually named short-term synaptic plasticity. Based on the model proposed in [29], [36] for short-term depression, and following previous studies concerning the dynamics of neural populations [16], we assume that evolves according to

(2)

where is the characteristic time scale of the STD mechanism, and is a parameter related with the reliability of the synaptic transmission. According to experimental measurements for these parameters in the somatosensory cortex of the rat [36], we set and unless specified otherwise. Assuming, for instance, a population time constant of , which would approximately correspond to the duration of a fast synaptic current pulse mediated by AMPA receptors, we obtain , which is within the physiological range measured in [36]. The last term on the right hand side of equation (2) is added to the original model in [36] to include some level of stochasticity in this, otherwise, deterministic description of synaptic transmission. The inclusion of such term constitutes a simple manner of considering the stochasticity due, for instance, to the unreliability of synaptic transmission [31], [32], the stochastic properties of receptor-transmitter interactions [37], the sparse connectivity of cortical circuits [38], [39], or other sources of noise not yet considered (see the Discussion Section for more details). The parameter controls the strength of this fluctuating term, and is a Gaussian white noise with zero mean and variance one.

Equations (1) and (2) constitute our minimal model of an excitatory neural network with stochastic depressing synapses. The simplifications assumed by such model allows to obtain some analytical derivations for the quantities of interest, and concretely for the probability distributions of permanence times in the up state, denoted by . Bistable systems in the presence of different sources of noise have been theoretically studied in detail in many works [40]–[44]. Here, however, we have employed a probabilistic approach which is very appropriate for the computation of the distribution of permanence times. In the following, we will derive an approximate expression for within this approach. First, we obtain the potential function and the conditions in which the dynamics of the system is driven by the variable . After that, we compute the probability distribution of ruin times of which, as we will see, leads to the probability distribution of permanence times in the up state, namely .

In order to compute the potential function of the dynamics (1,2) (namely ) one can see that, for realistic values of , the dynamics of is very slow compared to that of . We therefore can write equation (1) as

(3)

where we have adiabatically eliminated from the dynamics of . The extrema of are given by the solutions of the equation

(4)

In the following, we choose , with and . With this choice, one can easily check from equation (3) that the potential becomes symmetric in around when .

Equation (4) may have one or three solutions, depending on the slope of the hyperbolic tangent and on the value of . In order to obtain three solutions of (4) (that is, the bistable regime) the maximal slope of the hyperbolic tangent must be large enough, concretely the condition must be fulfilled. In addition, the threshold term must be not too small or too large so that has three crossing points with the straight line rather than one. This last condition can be written, as a first approach, as and , where are the values where the curvature of the hyperbolic tangent is maximal and minimal, respectively. The points can be easily computed from the third derivative of :

(5)

By setting we obtain

(6)

Using now these values for , the conditions and can be written as

(7)

which implies that, in order to have one maxima and two minima in , the variable must be in the range , where

(8)

From equation (8), one can see that the range of that allows to have three extrema in the potential is

(9)

The condition implies which is, therefore, a sufficient condition to obtain a double well potential for some value of . One can find, however, a small discrepancy between this approximate prediction and the actual properties of . The discrepancy appears because we have assumed that a sufficient condition for the existence of the three fixed point solutions of equation (4) is that and , and such assumption is only approximately correct. Plotting directly the potential as a function of reveals that the condition to obtain a double well potential for is , rather than .

Assuming that the above condition () is satisfied, three different shapes for the potential function can be found, as the figure 1A illustrates. When the potential function presents only one minimum, located around . Similarly, for the potential presents also a single minimum, but now located around . Finally, for the potential will take a double well shape, with the maximum being located around and the minima located around and , respectively.

It is worth noting that , with being the mean value of . Due to this, if the range is small compared with the fluctuations of , namely , the potential function will spend most of the time in the regimes and , with the double well regime appearing only when the system tries to jump from one of these regimes to the other (that is, when ). A direct consequence of this is that the mean firing rate will basically be switching between the up and down states (that is, and ), and that this switching will be driven by the dynamics of , as the figure 2 illustrates. Therefore, one expects that the distribution of permanence times of in the up (down) state, becomes approximately equal to the distribution of permanence times of in the () regime, as long as is satisfied. Due to this equivalence, in order to compute we only need to compute the distribution of permanence times of the variable in the regime, denoted as .

On the other hand, it should be noted that, since is a fraction of available neurotransmitters, its value should be kept within the range . In practice, this means that the value of must not be too large, so in order to make one has to restrict to small. In the results presented here, remain in its realistic range of values, and imposing *ad hoc* restrictions in such a way that is always within the range does not affect the results obtained here.

In order to compute the distribution of permanence times of in the (or ) regime, one can assume that the firing rate takes its mean value in equation (2). This is a reasonable approach since is much slower than for realistic values of the parameters. Considering this approach, and after the rescaling , equation (2) can be written as

(10)

which is the equation of the Ornstein-Uhlenbeck (OU) process (see [45] for details), with being the correlation time and . Therefore, computing the distribution of permanence times in the up state for our system is equivalent to obtain the distribution of the so called *ruin times* for the OU process [46], [47], which may be defined as follows: if we consider a stochastic process starting at from , the ruin time is the interval , where is the time at which returns to for the first time. Since is a stochastic process, the ruin times are stochastic quantities which follow a certain probability distribution.

The strategy employed here to calculate the distribution of ruin times is based on the relation between the ruin time and the *first passage time*, which is the typical time that a stochastic process needs to arrive at a certain threshold value when starting from a certain initial condition [47]. Because of the symmetry of the OU process, the distribution of ruin times are equivalent when considering excursions of the variable in the region or in the region. If we consider excursions in the region, we can set a small positive threshold near zero (that is, ), in such a way that the typical ruin time will be approximately equal to the corresponding first passage time, as the figure 1B illustrates. The excursions in the region typically lead to very short first passage times (since is too small) which we will not take into account in our calculations by considering only large enough ruin times.

The first passage time for the OU process with a small threshold can be performed by using the relation

(11)

where is the conditional probability distribution of the OU process, and is the first passage time distribution. This equation can be solved by taking into account the following property of the Laplace transformation

(12)

where is the Laplace transform of . By solving the Fokker-Planck equation associated with equation (10), one can obtain the conditional probability for the OU process

(13)

where , and being the standard deviation of . From expression (13), and assuming that is large enough (more precisely, assuming that , which is a valid hypothesis since most of the permanence times in the up state are much lower than ), one arrives at

(14)

We denote and . Employing the Laplace transformation in and the following expressions are obtained

(15)

Now, taking into account the property (12) in equation (11), the expression for is

(16)

Finally, for small one can approximate . With this approximation, one can easily perform the inverse Laplace transformation to equation (16) and obtain the distribution of first passage times for the OU process

(17)

A similar expression may be obtained if one considers more classical derivations of the first passage time of the OU process (see, for instance, [48]). In order to obtain the distribution of ruin times of the OU process, one has to consider a small (but positive) value of , which leads to . The distribution of ruin times of the variable , namely , and therefore, the distribution of permanence times in the up state, namely , for our system are also given by

(18)

which corresponds to a power-law probability distribution for .

Summarizing, the three following conditions must be fulfilled to obtain a power-law dependence in with exponent :

- Large enough values of . With this condition, we ensure that the dynamics of is much slower than that of .
- Large enough values of . In particular, we must have , according to the condition and the definitions and . This condition can be achieved even with very small values of , since can be arbitrarily small (by increasing , for instance).
- The condition must hold to ensure the existence of two well defined (up-down) states.

All these conditions may be easily achieved (up to some point) with realistic values of the model parameters, indicating that power-law distributions of the permanence times in the up state are plausible to be found in actual cortical media.

As we have stated in the previous sections, equations (1–2) govern the dynamics of our simplified neural system. A typical time series of the dynamics of this model, for the case of deterministic synapses (that is, ), is depicted in figure 2A. In this case, the mean firing rate of the population is characterized by a periodic switching between up and down states. This type of periodic behavior was already found and analyzed in previous theoretical studies [14], [16], [18] and yields bimodal histograms for the mean firing rate of the neural population (see figure 2B), as the experiments indicate [3]. However, these approaches ignore the stochastic nature of synaptic transmission, and other forms of stochasticity at the synaptic level, which seem to be crucial for information processing in neural systems [31], [32], [49]. Considering a certain level of synaptic stochasticity in addition to STD in our model, one obtains a qualitatively different emergent behavior, as is shown in figure 2C for . The mean firing rate presents then a complex switching between up and down states, and in particular involves a high variability in the permanence times in the up state.

When deterministic synapses are considered (that is, ) the dynamics of the mean firing rate becomes quasi periodic, as it was reported in [16], [18], [19], for instance. This type of dynamics naturally leads to exponential distributions for the permanence times. More precisely, for our model is similar (except for the term ) to the one analyzed in [16], which shows periodic oscillations of the network mean firing rate. In our case, however, the term introduces certain level of stochasticity which turns these periodic oscillations into quasi-periodic oscillations. This leads to the exponential distributions for the permanence times in the up state. When is increased, on the other hand, the stochasticity of the synapses leads to the appearance of power-law distributions for the duration of the up states. This behavior is shown in figure 3A, where low values of corresponds to exponential distributions for , while larger values of give as predicted by our theoretical calculations. Such power-law distributions may explain the high variability of permanence times in the up state, which has been observed in a number of *in vivo* experiments, such as in the cat visual cortex [7] and rat auditory cortex [21], to name a few. Interestingly, similar power-law dependences have been observed during sleep-wake transitions *in vivo* when one measures the distribution of permanence times in the wake state [22], [23]. On the other hand, exponential-like distributions, obtained for the case of having , are not able to explain this variability of the duration of up states.

By looking at the data for in figure 3A, one can observe the existence of a small deviation of the numerical results (blue points) with respect to the theoretically predicted slope (solid line) for very large values of . Such deviation is due to the fact that the separation of timescales between the dynamics of and (a necessary condition to obtain power-law dependences in ) is only approximate when considering realistic values of the parameters (and in particular, realistic values for ). More precisely, the approximation fails when the activity of the system falls in the occasional periods of very long permanence times in the up state (that is, for large enough values of , comparable with ). In order to study the effect of the separation of timescales between and , we have computed for different (increasing) values of while keeping fixed values for and (this can be done by properly modifying and with , respectively). As a consequence of this, the only effect of increasing will be a clearer separation of timescales between and . The results are shown in figure 3B, where one can see that larger values of (that is, clearer separation of timescales) lead to a displacement of the effective cut-off towards higher values of , as expected, and a clearer power-law distribution emerges.

It is worth noting that the appearance of an effective cut-off in for realistic conditions does not represent an unrealistic feature of the model, but rather it constitutes a prediction about the effective range of permanence times which are expected to occur in actual neural systems. Indeed, for realistic values of the parameters, our results predict permanence times in the up state up to , which is the maximum permanence time observed in experimental realizations [7]. Larger permanence times in the up state (of about seconds, for instance) should be expected to appear only as a consequence of input driven mechanisms (such as persistent activity associated with working memory tasks [50], [51]), and not as a consequence of spontaneous transitions between different voltage levels, which are the matter of interest in this work.

For a better characterization of the dynamics of the system, one can use, for instance, other statistical magnitudes such as the autocorrelation function of , which can be defined as

(19)

Here, indicates a temporal average. The autocorrelation function is depicted in figure 4A for the case of deterministic depressing synapses () and stochastic depressing synapses (). presents, for , two well located peaks at , which indicates a strong periodicity of the time series (as can be seen in figure 2A). On the contrary, the inclusion of a certain level of intrinsic stochasticity in the dynamics of introduces more pronounced temporal correlations in the dynamics of the system. This fact reflects the existence of long permanence stays in the up state, which occurs with more probability for high enough values of , as we have already discussed.

The spectral properties of the dynamics can be analyzed as well, via the power spectrum defined as

(20)

As one could expect, the power spectrum of the case presents a pronounced peak around a certain frequency, which in the particular case presented in the figure 4B is . The power spectrum for higher values of shows however different properties than the case . For instance, the figure 4B (which considers ) indicates an approximated power-law behavior for the power spectrum, with . This scale-free dependence can be understood by considering that, if is algebraic with exponent , the corresponding power spectrum becomes also algebraic with exponent , where the equation relates both exponents [44]. In our particular case, since , one obtains a theoretical prediction of for the exponent of the power spectrum. The theoretical relation between and exposed above, however, is only valid under the so called *single interval approximation*, which implies that the integration variable in equation (20) is smaller than the permanence time (see [44] for details). This condition does not strictly hold for our system (where ranges over several scales), and therefore it may introduce deviations in the theoretically predicted value of (which is around ) with respect to the value found in simulations (of around ).

Besides the level of synaptic stochasticity, i.e. , other parameters of the model could have an important effect on the dynamics as well. The parameter , for instance, controls the level of stochasticity of the dynamics of , and therefore one should expect that increasing its value could strongly influence the probability distribution . This is shown in figure 5A, where an increase of disrupts the appearance of power-law dependences, and exponential distributions appear instead. This change in is due to the fact that high levels of the additive noise make the system to jump more frequently from one state to the other, and therefore long stays in the up state (and thus distributions with long power-law tails) rarely occur.

The parameters involving the dynamics of also affect the probability distributions . The parameter , for instance, is responsible for the modulation of via the mean firing rate (see equation (2)), and therefore it can influence both the dynamics of and . As one may see in figure 5B, when takes low values a bump in emerges for high . Such deviation from the power-law dependence indicates that long stays in the up state occur more frequently than in the power-law case. Attending at equation (2), one can see that an increase of the mean firing rate decreases the variable via the parameter . Therefore, if takes lower values the decrement of will be smaller. As a consequence, the stays of in the regime (see Methods Section) will last longer, and the stays of the system in the up state will also last longer, causing the observed deviation from the power-law tendency. It should be noted, however, that the values of which allow the appearance of power-law dependences in for our model agree with the values of measured in actual cortical media where up and down transitions are observed [36].

We have also analyzed in detail the effect that varying has on the probability distribution of permanence times. Note that, contrarily to the previous study presented above, we have now varied the parameter while all the other parameters are kept fixed. This implies that the modification of will now have an effect on the separation of timescales between and , but also on the concrete value of and on the amplitude of the noisy term of equation (2) (namely ). The results are shown in figure 5C, where one can distinguish three different regimes as a function of the particular value of . For low (red region in the figure), the probability distributions show an exponential decay for large permanence times. The reason for this decay is that, for low , the variable does not perform long excursions in the region (see Methods Section), and therefore the probability to have large values of decreases and the power law behavior for is not obtained. As is increased, long excursions for begin to occur, and we obtain a power law behavior (green region in the figure). Finally, one can appreciate that, for even larger values of (blue region in the figure), the probability distribution of permanence times in the up state presents a power law dependence with being an increasing function of . Such dependence can not be explained by our previous theoretical predictions, based in the assumption that the system is in the bistable regime, and deserves a detailed analysis which will be exposed in the next section.

In the Methods Section, we established several conditions which had to be fulfilled in order to obtain power law dependences for . In particular, our previous analysis indicates that the condition must hold in order to have a potential function with three extrema (bistable regime). However, as we will see in the following, power law expressions for may appear even if the potential function has only one extremum in (concretely, one minimum), although the origin of such power law distributions is different from the one considered in previous sections, as we will see.

When (which occurs for or , for instance), the potential function has only one minimum in , whose location strongly depends on . An approximated expression for the location of this minimum as a function of can be obtained by expanding the hyperbolic tangent of the fixed point expression of (see equation (4)) around its argument (which is small in this limit), yielding

(21)

where is the value of which corresponds to the minimum of the potential function. Therefore as varies around , the location of the minimum of the potential also varies in the same way around . As an example, time series of both and are shown in figure 6A for a given set of parameters which satisfies . In this time series, the variable fluctuates around the value , which is fully determined by (that is, the variable becomes a slave variable of ). The predictions of equation (21) agree approximately well with simulations and with the numerical evaluation of the fixed points of equation (1), as the figure 6B shows.

Since behaves now as a stochastic variable which does not present a clear bistable dynamics, the numerical computation of the distribution of the permanence times will depend on the exact value of above which the system is considered to be in the up state. As we have seen before, this *threshold* value takes the form (see caption of figure 3), where usually may take a value between and . While the results presented for (that is, the bistable regime) are quite robust for different values of , in the regime this parameter has indeed some effect on , which indicates the difficulty to accurately analyze the up and down dynamics in this case.

In figure 7A, one observes that the distribution shows also a power law behavior for and different values of , for a set of parameter values which satisfies (that is the monostable regime). The concrete value of depends strongly on and it has also a weaker dependence with , as the figure 7B illustrates. This type of power-law behavior appearing in the monostable regime corresponds to the blue region in figure 5C, as well.

It is worth noting that actual recordings of up and down transitions does not present a clear distinction between up and down states, and several nontrivial methods are commonly employed to discriminate between both states [52]. Therefore, the results found for the regime could indeed reflect the behavior of actual cortical up-down transitions, showing power law dependences in with and indicating that the concrete nature of the transitions is a synaptic-driven monostable dynamics.

For a complete characterization of the model, one can summarize all the observed behaviors in a phase plot such as the one presented in figure 8A. A total of four different behaviors can be found in the space. The first one concerns the dynamics of whose permanence times in the up state follows an exponential distribution (labeled as “E” in the figure). If the noise amplitude is sufficiently high, one can increase the value of to reach the regime “C”, in which the dependence is obtained. By increasing even more, the probability distribution takes the form , with (regime denoted by “S”), as we have already seen in figure 6. Finally, we also observe that when the depression time scale is not large enough (and ), a regime of quasi-periodic time series of is obtained, with a well-defined duration of up states (regime denoted by “P”). The lines between the different regimes have been obtained by visual inspection of for different values of and . In particular, the regime “P” is characterized by the appearance of a *bump* in the probability distribution for some value of (which reflects a preferred duration of the up state), and the existence of such bump has been used as a criterion to distinguish between regimes “P” and “E”. Similarly, we assumed that the regimes “C” and “S” correspond to the situation in which a power-law behavior that extends for two decades or more is found for . Such criterion, together with an estimation of the slope of the power-law via standard Levenberg-Marquardt fitting algorithms, allows to distinguish between regimes “E”, “C” and “S”.

It must be clarified, however, that actual up and down cortical transitions present most likely a richer repertoire of dynamical regimes than the one obtained with our simplified model. It is known, for instance, that attractor neural networks with dynamic synapses may exhibit different dynamics corresponding to memory, non-memory and switching regimes [18], [19]. In this work, we have extensively explored different regimes of switching behavior, and its implications for the up and down dynamics observed in the cortex. The memory and non-memory regimes, however, can be also found in our simplified model by assuming that . After taking these limits, the system will be in the memory regime if the potential function is bistable, or in the non-memory regime if is monostable.

We have shown that the experimentally observed large fluctuations in up and down permanence times can be explained as the result of sufficiently noisy dynamical synapses with sufficiently large recovery times. Our study suggests that a power-law distribution for these permanence times may emerge as a consequence of these two ingredients. Static synapses cannot account for this behavior, nor can dynamical synapses in the absence of noise.

The origin of up and down cortical transitions is still unclear, although different factors that may influence their occurrence have been recently reported. It is known, for instance, that inhibitory GABAergic currents strongly contribute to the temporal coding and spike timing precision of cortical networks during up states of activity [3], [53], [54]. Several modeling studies also show the relevance of inhibitory interneurons in the generation of many types of oscillations in the brain (see for instance [55]). However, other studies indicate that most of the main features of up and down transitions depends strongly on synaptic plasticity mechanisms, both of long-term and short-term ones [16], [56], and that the transitions appear even in the absence of inhibition [16]. In this work we have made the common assumption that the effects of inhibition can be treated as additive and can be incorporated in the threshold of the neuron. This is known to be a valid approximation in mean field neural network analysis, but may fail when precise timing and details of the dynamical aspects of the neuron affect the inhibition [57], [58].

Regarding to synaptic characteristics, recent works show that synaptic fluctuations could have an important role in the generation of transitions between up and down states [14], [59], [60]. Since our model introduces stochasticity in the synaptic dynamics in a highly simplified manner, however, the last term in equation (2) should not be associated only with ureliability in synaptic transmission. Indeed, we have assumed that other sources of stochasticity may be contributing to this fluctuating term in the mean-field quantities and . For instance, it is widely known that connectivity in actual cortical media is highly sparse. Such feature implies that, in order to obtain the mean-field quantity , the average over synapses must be performed over a number of synapses, with ranging over connections per neuron [38]. In this situation, the fluctuations of would be of order , which leads to a range of for the values given above. As we have seen, our results state that a value of is enough to obtain power-law distributions (see figure 3), which lies within this range. Therefore, topology-induced fluctuations constitute an important source of stochasticity which could be responsible of the appearance of power-law distributions in . Other sources of stochasticity at synaptic level, such as the stochastic properties of receptor-transmitter interactions, may also contribute to the last term of equation (2). Moreover, the low activity rates typical from cortical media lead to a poor time-averaging of the incoming input, and therefore the fluctuations at the postsynaptic level will be large at these short-time scales (of the order of the typical synaptic integration time constant).

On the other hand, the amplitude of the noisy term, , does not need to be very high to induce the appearance of power-law distributions in . As we have stated above, a sparse connectivity already induces a level of stochasticity which is within the desired range, for instance. Furthermore, the noisy term could even be arbitrarily small: attending to our theoretical predictions, a necessary condition to have power-law distributions is that fluctuations of must be much larger than (see Methods Section). Since may be lowered to arbitrary levels (by increasing , for instance), even a small noisy term in the dynamics of may induce power-law distributions.

It is also known that short-term synaptic mechanisms, such as short-term depression and facilitation, usually play a role in the efficient processing of information. In particular, they may be relevant in many tasks, such as in signal detection and coding [29], [61]–[63] or switching between different activity patterns previously stored [19], [64]. However, their role on the transitions between cortical states has been pointed out only by a few studies [15], [16], [65], and their possible effects on the statistics of the transitions, which is the focus of our work, have been ignored. To the best of our knowledge, the present study is the first one which analyzes, even in a simplified manner, the strong effect of synaptic stochasticity– in a general sense– and dynamic synapses in the statistics of the up and down transitions. The possible role of other short-term synaptic mechanisms, such as STF, has not been addressed yet and constitutes a interesting issue still open.

In our analysis we assumed that the dynamics is symmetric in the up and down states. This is in contradiction with experimental evidences [66] which shows that power-law distributions are obtained for permanence times in the up state, while permanence times in the down state are exponentially distributed. However, this discrepancy disappears when one considers a more realistic transduction function which gives an asymmetric potential for the dynamics, and as a consequence the up-down symmetry is broken. More detailed studies considering, for instance, some of the biologically realistic aspects discussed above, should be performed to test our predictions. In particular, a more elaborated study considering realistic neuron models (such as Hodgkin-Huxley model [67]) and stochastic STD models (see [29], [49], for instance) is necessary, as well as more detailed experimental studies which may confirm our predictions.

From a general point of view, evidences of criticality have been recently found in an increasing number of neural systems, such as in the functional connectivity of the living human brain [24], in critical avalanches of neuronal activity [25], or in sleep-wake transitions [23], to name a few. According to the results presented in this work, transitions between up and down cortical states could also present some relevant properties typical of systems at criticality. Some of these properties have been already measured in experiments, such as a high sensitivity of the system to external stimuli [8], or the presence of power-law dependences in the power spectra of the neural dynamics [53].

It is worth noting that other kind of probability distributions for , such as a log-normal distribution, could also satisfactorily explain the irregularity in the up states found in experiments. Our study shows the importance of some biophysical factors, such as the neurotransmitter recovery time and the inherent synaptic stochasticity, and predicts a power-law dependence on as a consequence of such factors. However, further study is needed to investigate other mechanisms, not taken in account in this work, which could influence the permanence times in the up state. In a more general sense, our results may proportionate a new perspective of the phenomena of up and down transitions (and a theoretical framework) that could serve to conciliate the main experimental findings, and that could help for a deep understanding of this complex dynamics of the brain activity.

We thank useful discussions with Sebastiano de Franciscis and Samuel Johnson.

**Competing Interests: **The authors have declared that no competing interests exist.

**Funding: **This work was supported by the MEyC-FEDER project FIS2009-08451 and the Junta de Andalucia project P07–FQM–02725. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

1. Steriade M, McCormick DA, Sejnowski TJ. Thalamocortical oscillations in the sleeping and aroused brain. Science. 1993;262:679–685. [PubMed]

2. Arieli A, Sterkin A, Grinvald A, Aertsen A. Dynamics of ongoing activity: explanation of the large variability in evoked cortical responses. Science. 1996;273:1868–1871. [PubMed]

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

4. McCormick DA. Neuronal networks: flip-flops in the brain. Curr Biol. 2005;15:R294–R296. [PubMed]

5. Steriade M, Nunez A, Amzica F. A novel slow (<1hz) oscillation of neocortical neurons in vivo: depolarizing and hyperpolarizing components. J Neurosci. 1993;13:3252–3265. [PubMed]

6. Steriade M, Nunez A, Amzica F. Intracellular analysis of relations between the slow (<1hz) neocortical oscillation and other sleep rhythms of the electroencefalogram. J Neurosci. 1993;13:3266–3283. [PubMed]

7. Lampl I, Reichova I, Ferster D. Synchronous membrane potential fluctuations in neurons of the cat visual cortex. Neuron. 1999;22:361–374. [PubMed]

8. Anderson J, Lampl I, Reichova I, Carandini M, Ferster D. Stimulus dependence of two-state fluctuations of membrane potential in cat visual cortex. Nat Neurosci. 2000;3:617–621. [PubMed]

9. Petersen CCH, Hahn TTG, Mehta M, Grin-vald A, Sakmann B. Interaction of sensory responses with spontaneous depolarization in layer 2/3 barrel cortex. Proc Nat Acad Sciences USA. 2003;100:13638–13643. [PubMed]

10. Cossart R, Aronov D, Yuste R. Attractor dynamics of network up states in the neocortex. Nature. 2003;423:283–288. [PubMed]

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

12. Mayor G, Tank D. Persistent neural activity: prevalence and mechanisms. Curr Opin Neurobiol. 2004;14:675–684. [PubMed]

13. Loewenstein Y, Mahon S, Chadderton P, Kitamura K, Sompolinsky H, et al. Bistability of cerebellar Purkinje cells modulated by sensory stimulation. Nat Neurosci. 2005;8:202–211. [PubMed]

14. Parga N, Abbott LF. Network model of spontaneous activity exhibiting synchronous transitions between up and down states. Front Neurosci. 2007;1:57–66. [PMC free article] [PubMed]

15. Steriade M, Timofeev I, Grenier F. Natural waking and sleep states: a view from inside neocortical neurons. J Neurophysiol. 2001;85:1969–1985. [PubMed]

16. Holcman D, Tsodyks M. The emergence of up and down states in cortical networks. PLoS Comput Biol. 2006;2:174–181. [PubMed]

17. Wang XJ. Synaptic basis of cortical persistent activity: the importance of nmda receptors to working memory. J Neurosci. 1999;19:9587–9603. [PubMed]

18. Pantic L, Torres JJ, Kappen HJ, Gielen SCAM. Associative memory with dynamic synapses. Neural Comput. 2002;14:2903–2923. [PubMed]

19. Torres JJ, Cortes J, Marro J, Kappen H. Competition between synaptic depression and facilitation in attractor neural networks. Neural Comp. 2007;19:2739–2755. [PubMed]

20. Melamed O, Barak O, Silberberg G, Markram H, Tsodyks M. Slow oscillations in neural networks with facilitating synapses. J Comput Neurosci. 2008;25:308–316. [PubMed]

21. Deco G, Marti D, Ledberg A, Reig R, Sanchez-Vives MV. Effective reduced diffusion-models: a data driven approach to the analysis of neuronal dynamics. PLoS Comput Biol. 2009;5:e1000587. [PMC free article] [PubMed]

22. Pearlmutter BA, Houghton CJ. A new hypothesis for sleep: tuning for criticality. Neural Comput. 2009;21:1622–1641. [PubMed]

23. Lo CC, Chou T, Penzel T, Scammell TE, Strecker RE, et al. Common scale-invariant patterns of sleep-wake transitions across mammalian species. Proc Nat Acad Sci USA. 2004;101:17545–17548. [PubMed]

24. Eguiluz VM, Chialvo DR, Cecchi GA, Baliki M, Apkarian AV. Scale-free brain functional networks. Phys Rev Lett. 2005;94:018102. [PubMed]

25. Beggs JM, Plenz D. Neuronal avalanches in neocortical circuits. J Neurosci. 2003;23:11167–11177. [PubMed]

26. Lazar A, Pipa G, Triesch J. Fading memory and time series prediction in recurrent networks with different forms of plasticity. Neural Networks. 2007;20:312–322. [PubMed]

27. Gomez V, Kaltenbrunner A, Lopez V, Kappen HJ. Self-organization using synaptic plasticity. 2009. pp. 513–520. In: Advances in neural information processing systems. volume 22.

28. Lazar A, Pipa G, Triesch J. SORN: a self-organizing recurrent neural network. Front Comput Neurosci. 2009;3 doi: 10.3389/neuro.10.023.2009. [PMC free article] [PubMed]

29. Abbott LF, Valera JA, Sen K, Nelson SB. Synaptic depression and cortical gain control. Science. 1997;275:220–224. [PubMed]

30. Zador AM, Dobrunz LE. Dynamic synapses in the cortex. Neuron. 1997;19:1–4. [PubMed]

31. Dobrunz LE, Stevens CF. Heterogeneity of release probability, facilitation, and depletion at central synapses. Neuron. 1997;18:995–1008. [PubMed]

32. Zador A. Impact of synaptic unreliability on the information transmitted by spiking neurons. J Neurophysiol. 1998;79:1219–1229. [PubMed]

33. Wilson HR, Cowan JD. Excitatory and inhibitory interactions in localized populations of model neurons. Biophys Journal. 1972;12:1–24. [PubMed]

34. Gerstner W. Population dynamics of spiking networks: fast transients, asynchronous states, and locking. Neural Comput. 2000;12:43–89. [PubMed]

35. Gerstner W, Kistler W. Spiking neuron models- Single neurons, population, plasticity. Cambridge Univ. Press; 2002.

36. Tsodyks MV, Markram H. The neural code between neocortical pyramidal neurons depends on neurotransmitter release probability. Proc Natl Acad Sci USA. 1997;94:719–723. [PubMed]

37. Faber DS, Young WS, Legendre P, Korn H. Intrinsic quantal variability due to stochastic properties of receptor-transmitter interactions. Science. 1992;258:1494–1498. [PubMed]

38. Kandel ER, Schwartz JH, Jessell TM. Principles of neural science, 4th ed. McGraw-Hill, New York; 2000.

39. Scannell JW, Blakemore C, Young MP. Analysis of connectivity in the cat cerebral cortex. J Neurosci. 1995;15:1463–1483. [PubMed]

40. Pechukas P, Hanggi P. Rates of activated processes with fluctuating barriers. Phys Rev Lett. 1994;73:2772–2775. [PubMed]

41. Madureira AJR, Hanggi P, Buonomano V, Rodrigues WA., Jr Escape from a fluctuating double well. Phys Rev E. 1995;51:3849–3861. [PubMed]

42. Ping Z. Dynamical properties o a bistable system driven by cross-correlated additive and multiplicative colored noises. Chines J Phys. 2006;44:117–126.

43. Dong-Xi L, Wei X, Yong-Feng G, Gao-Jie L. Transient properties of a bistable system with delay time driven by non-gaussian and gaussian noises: mean first-passage time. Commun Theor Phys. 2008;50:669–673.

44. Tu Y, Grinstein G. How white noise generates power-law switching in bacterial flagellar motors. Phys Rev Lett. 2005;94:208101. [PubMed]

45. van Kampen NG. Stochastic processes in physics and chemistry. North-Holland; 1990.

46. Chandrasekhar S. Stochastic problems in physics and astronomy. Rev Mod Phys. 1943;15:1–89.

47. Newman MEJ. Power laws, Pareto distributions and Zipf's law. Contemporary Phys. 2005;46:323–351.

48. Alili L, Patie P, Pedersen JL. Representations of the first hitting time density of an Ornstein-Uhlenbeck process. Stochastic Models. 2005;21:967–980.

49. de la Rocha J, Parga N. Short-term synaptic depression causes a non-monotonic response to correlated stimuli. J Neurosci. 2005;25:8416–8431. [PubMed]

50. Fuster JM, Alexander GE. Neuron activity related to short-term memory. Science. 1971;173:652–654. [PubMed]

51. Goldman-Rakic PS. Cellular basis of working memory. Neuron. 1995;14:477–485. [PubMed]

52. Seamari Y, Narvaez JA, Vico FJ, Lobo D, Sanchez-Vives MV. Robust off- and online separation of intracellularly recorded up and down cortical states. PLoS One. 2007;2:e888. [PMC free article] [PubMed]

53. Hasenstaub A, Shu Y, Haider B, Kraushaar U, Duque A, et al. Inhibitory postsynaptic potentials carry synchronized frequency information in active cortical networks. Neuron. 2005;47:423–435. [PubMed]

54. Compte A. Computational and in vitro studies of persistent activity: edging towards cellular and synaptic mechanisms of working memory. Neurosci. 2006;139:135–151. [PubMed]

55. Brunel N. Dynamics of sparsely connected networks of excitatory and inhibitory spiking neurons. J Comp Neurosci. 2000;8:183–208. [PubMed]

56. Kang S, Kitano K, Fukai T. Structure of spontaneous up and down transitions self-organizing in a cortical network model. PLoS Comput Biol. 2008;4:e100022. [PMC free article] [PubMed]

57. Rotstein HG, Pervouchine DD, Acker CD, Gillies MJ, White JA, et al. Slow and fast inhibition and an h-current interact to create a theta rhythm in a model of CA1 interneuron network. J Neurophysiol. 2005;94:1509–1518. [PubMed]

58. Vierling-Claassen D, Siekmeier P, Stufflebeam S, Kopell N. Modeling GABA alterations in schizophrenia: a link between impaired inhibition and altered gamma and beta range auditory entrainment. J Neurophysiol. 2008;99:2656–2671. [PMC free article] [PubMed]

59. Cortes JM, Torres JJ, Marro J, Garrido PL, Kappen HJ. Effects of fast presynaptic noise in attractor neural networks. Neural Comput. 2006;18:614–633. [PubMed]

60. Johnson S, Marro J, Torres JJ. Functional optimization in complex excitable networks. Europhys Lett. 2008;83:46006.

61. Puccini GD, Sanchez-Vives MV, Compte A. Selective detection of abrupt input changes by integration of spike-frequency adaptation and synaptic depression in a computational network model. J Physiol-Paris. 2006;100:1–15. [PubMed]

62. Mejias JF, Torres JJ. The role of synaptic facilitation in spike coincidence detection. J Comput Neurosci. 2008;24:222–234. [PubMed]

63. Mejias JF, Torres JJ. Emergence of resonances in neural systems: the interplay between adaptive threshold and short-term synaptic plasticity. 2010. Submitted (arXiv:09060756) [PMC free article] [PubMed]

64. Mejias JF, Torres JJ. Maximum memory capacity on neural networks with short-term synaptic depression and facilitation. Neural Comp. 2009;21:851–871. [PubMed]

65. Timofeev I, Grenier F, Steriade M. Impact of intrinsic properties and synaptic factors on the activity of neocortical networks in vivo. J Physiol-Paris. 2000;94:343–355. [PubMed]

66. de Franciscis S, Torres JJ, Sanchez-Vives MV. 2010. In preparation.

67. Hodgkin AL, Huxley AF. A quantitative description of membrane current and its application to conduction and excitation in nerve. J Physiol. 1952;117:500–544. [PubMed]

Articles from PLoS ONE are provided here courtesy of **Public Library of Science**

PubMed Central Canada is a service of the Canadian Institutes of Health Research (CIHR) working in partnership with the National Research Council's national science library in cooperation with the National Center for Biotechnology Information at the U.S. National Library of Medicine(NCBI/NLM). It includes content provided to the PubMed Central International archive by participating publishers. |