Email updates

Keep up to date with the latest news and content from BMC Systems Biology and BioMed Central.

Open Access Research article

An optimally evolved connective ratio of neural networks that maximizes the occurrence of synchronized bursting behavior

Chao-Yi Dong1 and Kwang-Hyun Cho2*

Author Affiliations

1 Department of Automatic Control, Inner Mongolia University of Technology, Huhhot 010080, People's Republic of China

2 Department of Bio and Brain Engineering, Korea Advanced Institute of Science and Technology (KAIST), Daejeon 305-701, Republic of Korea

For all author emails, please log on.

BMC Systems Biology 2012, 6:23  doi:10.1186/1752-0509-6-23

The electronic version of this article is the complete one and can be found online at: http://www.biomedcentral.com/1752-0509/6/23


Received:14 December 2011
Accepted:31 March 2012
Published:31 March 2012

© 2012 Dong and Cho; licensee BioMed Central Ltd.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

Background

Synchronized bursting activity (SBA) is a remarkable dynamical behavior in both ex vivo and in vivo neural networks. Investigations of the underlying structural characteristics associated with SBA are crucial to understanding the system-level regulatory mechanism of neural network behaviors.

Results

In this study, artificial pulsed neural networks were established using spike response models to capture fundamental dynamics of large scale ex vivo cortical networks. Network simulations with synaptic parameter perturbations showed the following two findings. (i) In a network with an excitatory ratio (ER) of 80-90%, its connective ratio (CR) was within a range of 10-30% when the occurrence of SBA reached the highest expectation. This result was consistent with the experimental observation in ex vivo neuronal networks, which were reported to possess a matured inhibitory synaptic ratio of 10-20% and a CR of 10-30%. (ii) No SBA occurred when a network does not contain any all-positive-interaction feedback loop (APFL) motif. In a neural network containing APFLs, the number of APFLs presented an optimal range corresponding to the maximal occurrence of SBA, which was very similar to the optimal CR.

Conclusions

In a neural network, the evolutionarily selected CR (10-30%) optimizes the occurrence of SBA, and APFL serves a pivotal network motif required to maximize the occurrence of SBA.

Background

In the brain development, neurons are assembled together via numerous synapses to build up complicated neuronal networks performing specific behaviors, such as transient or sporadic activity, synchronized bursting activity (SBA), and hyper-excitable activity. One of the most prominent behaviors in cortical networks is the synchronized bursting spikes occurring in the brain development and maturation [1-3]. The behavior is not only found in ex vivo cultured cortical networks [4] but also in the brain regions of several in vivo animal models like visual cortex [5], hippocampus [6], and auditory neocortex [7]. In particular, under in vivo conditions, SBA is considered highly related to a variety of crucial biophysical functions, such as attentional selection [8-10], cognitive motor processes [11], visual pattern recognition [12], auditory object perception [13].

Although SBA is an unique phenomenon in neuronal networks, characteristics of the neural networks causing SBA remain unknown, in contrast to the study on the function significance of the SBA [14]. Presently, large random ex vivo cortical networks are more appropriate experimental model systems in the studies on the universal mechanisms governing the formation and conservation of neural network activities. Experiments using ex vivo cultured neural networks have demonstrated that the adjustment of synaptic connections is highly correlated with the development of neuronal network behavior such as the evolution of spontaneous electrical activity [15,16]. In the matured phase of an ex vivo cultured neural network, each neuron builds up synaptic connections with 10-30% of other neurons within the neural network [17,18]. Another line of evidence has indicated that the electrical activity of neurons can directly affect the outgrowth of neurites, and such reconfiguration of neuronal networks in turn causes adaptive adjustment of the neuronal electrical activity [19]. This behavior-dependent regulatory mechanism precisely drives and controls networks to grow, prune, and finally converge to a proper connective ratio (CR) (10-30%).

According to the above described connectivity characteristics of ex vivo cultured neural networks, two interesting questions arise: why do such a matured neural network keep its CR within a fixed range (10-30%); and what biological significance and associated implications does this fixed CR have? To answer these intriguing questions, we hypothesize that the CR is associated with the facilitation of synchronized bursting network behaviors, since synaptic connections are always found correlated with network behaviors in ex vivo experiments. Spike response models [20-23] were used to construct randomly connected artificial pulsed neural networks. The connective weights between two neurons were randomly selected, and the CR of the networks was increased progressively to mimic the process of development of cultured neural networks. The correlation between network behavior and structure was investigated using simulations. Subjecting the simulations to parameter perturbations revealed that, for a network with an excitatory ratio (ER) at 80-90% (a realistic ratio for ex vivo networks), the CR of the network always lies in a range of 10-30% when the occurrence of SBA reaches its highest expectation. This value is consistent with the matured CR of ex vivo neuronal networks with the inhibitory synaptic ratio at 10-20% [24,25]. This result reveals that the networks are evolved to form such a CR for optimizing the occurrence of SBA rather than randomly connected.

This study also explored the relationship between the occurrence of SBA and the composition of network motifs in the neural networks [26,27]. We found that SBA can be found only in the networks containing an all-positive-interaction feedback loop (APFL) [2]. For networks containing APFLs, the number of APFLs also demonstrates an optimal range corresponding to the maximized occurrence of SBA, close to the CR. Thus, we infer that the APFL may serve a crucial network motif underlying to maximize the occurrence of SBA.

For a pilot study in real neural networks, we have employed the neural network of nematode worm C. elegans [28,29]. The nervous system of C. elegans consists of 302 neurons and the number of neurons is almost same for different individuals. Each neuron in C. elegans' nervous system has distinct properties in view of morphology, connectivity, and position, and therefore it can be labelled specifically. The neural network of C. elegans is highly clustered like regular lattices and also has small characteristic path lengths like random graphs. So, it is well represented by small-world networks [30,31]. We investigated the egg-laying circuit of C. elegans including 11 neurons or neuron classes to examine our major claims [32,33]. As a result, we found that the egg-laying circuit has 17.3% CR and 10.5% ER which lie within the aforementioned evolved ranges. We also found that three two-node APFLs included in this circuit contribute to inducing a much higher level of SBAs in contrast to the randomly connected networks with the same number of network nodes.

Results

The optimal CR at the maximal occurrence of SBA

To unravel the biological significance of the CR of matured neural networks (10-30%), we first investigated the relationships between CR, ER, and the occurrence of SBA. Spike response models (SRM) can be used to simulate random ex vivo cortical networks so that their fundamental dynamical properties can be modeled [34]. The detailed simulation protocols for the artificial pulsed neural networks constructed by SRM are introduced in the Methods section. The connective weights among neurons were randomly assigned in a certain range to obtain a result that was irrelevant to the specific value of connective weights. For each CR and ER, 1,000 randomly connected pulsed neural networks were constructed to generate various network behaviors which were then further classified into four major categories with respect to the proposed criteria in the section entitled "The typical behaviors of spike neural networks".

The SBA properties of networks were investigated at two different scales: small networks with 12 nodes and large networks with 60 nodes. We recorded and calculated the expectation and standard deviation of SBA occurrence over the 1,000 networks (with a stereotyped CR and ER). Figure 1 shows the relationships among CR, ER, and the occurrence of SBA in the 12-node networks (Figure 1, a-1, a-2, a-3) and the 60-node networks (Figure 1, b-1, b-2, b-3). The mean and standard deviation of the occurrence of SBA were calculated over the 1,000 spike neural networks in each possible combination of CR and ER. The interesting profiles concerning the relationship between CR and the occurrence of SBA are found when the ER was equal to 0.9, which coincides with the experimental observation in cultured neural networks [25]. Figures 1 a-3 and b-3 show an optimal CR of 15% in the 12-node networks and 10% in the 60-node networks, respectively. This relation between CR and SBA suggests that the direction of evolutionary selection of CRs (not only in artificial pulse neural networks but also in in vitro neural networks) is to maximize the possibility for synchronized bursting behavior by networks.

thumbnailFigure 1. The relationship between the occurrence of SBA and CR (or ER) for 12-node networks and 60-node networks. (a-1), (a-2) The relationships among CR, ER, and the occurrence of SBA for 12-node networks from two different perspectives. The occurrence of SBA in 1,000 randomly connected networks was evaluated on each point of a lattice composed of CRs (0:0.05:1) and ERs (0:0.1:1). The surface points represent the mean value of the occurrence of SBA, while the upper and lower red bars show their standard deviations over 1,000 simulations. (a-3) The relationship between CR and the occurrence of SBA for 12-node networks when ER equals 0.9, which is consistent with the actual biological level. (b-1), (b-2) The relationship among CR, ER, and the occurrence of SBA for 60-node networks from two different perspectives. The occurrence of SBA in 1,000 randomly connected networks was evaluated on the same lattice as (a). The surface points represent the mean value of the occurrence of SBA, while the upper and lower red bars show their standard deviations over 1,000 simulations. (b-3) The relationship between CR and the occurrence of SBA for 60-node networks when ER equals 0.9 (the biological level).

The optimal number of APFLs causing maximal occurrence of SBA

The existence of APFLs was shown to be a prerequisite for inducing SBA for 2, 3, and 4-node pulsed neural networks [2]. SRM simulations (2, 3, and 4-node small networks) indicated that only those networks containing APFLs could produce SBA. Other networks containing only negative feedback loops or double-negative positive feedback loops (referred to the section entitled "Definition of network motifs and feedback loops" for its definition) and those without any feedback loop, could not generate SBA irrespective of synaptic efficacy. We were intrigued by these observations and further investigated the relationship between APFLs and the occurrence of SBA in larger-scale networks. Figure 2 provides detailed descriptions of this relationship in both 12-node and 60-node networks. Statistical tests on these relationships were carried out for all 1,000 simulations on the same lattices, ERs (0:0.05:1) × CRs (0:0.1:1), shown in Figure 1 a-1, a-2 and Figure 1 b-1, b-2.

thumbnailFigure 2. The relationship between the total number of 2, 3, and 4-node APFL motifs and the occurrence of SBA or HEA for the 12-node networks and 60-node networks. (a-1) The relationship between the total number of 2, 3, and 4-node APFL motifs and the occurrence of SBA for the 12-node pulsed neural networks. (a-2) The relationship between the total number of 2, 3, and 4-node APFL motifs and the occurrence of 2-channel HEA for the 12-node networks. (b-1) The relationship between the total number of 2, 3, and 4-node APFL motifs and the occurrence of SBA for the 60-node pulsed neural networks. (b-2) The relationship between the total number of 2, 3, 4-node APFL motifs and the occurrence of 2-channel HEA for the 60-node networks.

Figures 2 a-1 and a-2 show the correlation between two typical network behaviors, SBA and hyper-excitable activity (HEA), and the total number of 2, 3, and 4-node APFL motifs in 12-node pulsed neural networks (see Figure 2 a-1, a-2). When the number of APFL motifs increases, the mean of SBA occurrence initially increases, then reaches a peak (maximum SBA occurrence of 4.3, corresponding to four APFLs in the 12-node networks), and finally returns to zero (after 37 APFLs in the 12-node networks, see Figure 2 a-1). However, the occurrence of HEA always increased with the increase in the number of APFL motifs (see Figure 2 a-2). If the number of APFLs exceeded 37 in a 12-node network, the occurrence of 2-channel HEAs (see the section of 'Simulation protocols" for its definition) approached the maximum, <a onClick="popup('http://www.biomedcentral.com/1752-0509/6/23/mathml/M1','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/6/23/mathml/M1">View MathML</a>, which corresponds to the case in which all pairs of channels are hyper-excitable. Thus, the decreased the occurrence of SBA can be explained by the increase of HEA occurrence when APFLs are sufficiently enriched in a network. A similar result was obtained in 60-node networks (see Figures 2 b-1 and b-2). The maximum 2-channel SBA occurrence was 727, corresponding to 400 APFL motifs (2, 3, 4-node) in the 60-node networks. When the number of 2, 3, 4-node APFL motifs exceeded 900, HEA fully dominated, and all other behaviors including SBA vanished from the 60-node networks. Notably, the first points (mean = 0 and standard deviation = 0) in Figures 2 a-1 and b-1 imply that no SBA occurs when no APFL motif is present in the network. Therefore, the APFLs is necessary to trigger SBA in a pulsed neural network.

Figures 1 and 2 clearly indicate that SBA occurs significantly within an optimal range of CR and APFL number. In fact, the simulations demonstrate that the number of APFL motifs increases along with the increase of ER and CR in randomly generated synthetic networks (data not shown). Therefore, we infer that the primary factor inducing the maximal occurrence of SBA may be the formation of a suitable number of APFLs in neural networks.

The relationship between the number of 2, 3, or 4-node APFLs and the occurrence of SBA

We also investigated the relationship between the distribution of each type of APFL and the level of SBA. In the 12-node neural networks, we found that 2-node APFLs are significantly enriched compared to 3-node or 4-node APFLs for all levels of SBA (Figure 3 a-1). Figure 3 b-1 shows the relationship between the number of 2, 3, or 4-node APFLs and the occurrence of SBA for 60-node pulsed neural networks. The number of 2-node APFLs significantly exceeds the number of 3-node or 4-node APFLs only when the occurrence of SBA exceeds a high level (1,200 events). Therefore, the number of 2-node APFLs dominates the total number of APFLs subject to a high level of SBA occurrence in both small- and large-scale networks.

thumbnailFigure 3. The influence of 2, 3, or 4-node APFL motifs on inducing SBA. (a-1) The number of 2, 3, or 4-node APFL motifs with respect to the level of SBA (12-node networks). (a-2) The occurrence of SBA along with the absence of particular APFL motifs or combinations of motifs in 12-node networks. (b-1) The number of 2, 3, or 4-node APFL motifs with respect to the level of SBA (60-node networks). (b-2) The occurrence of SBA along with the absence of particular APFL motifs or their combinations in 60-node networks.

How is SBA inhibited when each type of APFL motif is absent from the pulsed neural networks? Figure 3 a-2 shows the occurrence of SBA when one type or a combination of APFL motifs is not present in the 12-node networks. With the exclusion of more types of APFL motifs, both the mean and standard deviation approached zero. For example, if we take the absence of 2, 3, 4, 5, and 6-node APFLs into account, the occurrence of SBA is only 3.418 × 10-5 ± 0.0101 (mean frequency and standard deviation). Thus, the loss of more types of APFL motif gradually inhibits the occurrence of SBA. In 60-node networks, the observed trend slightly differs in that the exclusion of 2-node APFLs completely prohibits the occurrence of SBA (see the first point in Figure 3 b-2). This fact implies that 2-node APFLs may function dominantly in the inhibition of SBA in large-scale networks, compared with 3-node and 4-node APFLs.

A case study of the egg-laying circuit of C. elegans

We investigated the egg-laying circuit of C. elegans composed of 11 neurons or neuron classes for a pilot study of real neural networks [32,33]. We identified that this circuit contains three two-node APFLs as indicated by the bi-directional red arrows (Figure 4). In addition, we found that the CR of this circuit is 17.3% and the ER of it is 10.5%.

thumbnailFigure 4. The egg-laying circuit of C. elegans. Each node represents a neuron or a neuron class. The arrows represent excitatory synaptic connections and the line bars denotes inhibitory synaptic connections. In this network diagram, all outer connections are excluded for simplicity.

We carried out simulations over two network groups for 20,000 times with different synaptic weight perturbations. One group of networks are randomly connected with any CR between 0 and 1, and the other group of networks have the same topological structure as shown in Figure 4. Two-sample t-test was carried out for these two groups to examine whether the real biological neural network induce higher SBAs compared to the randomly connected neural networks. Details are as follows: Let vector x denote the SBAs of 20,000 randomly connected networks and vector y represent the SBAs of the egg-laying neural network of C. elegans with the 20,000 synaptic weight perturbations. Assuming that the variances of x and y are unknown, it becomes the Behrens-Fisher problem [35]. The t-test statistic is <a onClick="popup('http://www.biomedcentral.com/1752-0509/6/23/mathml/M2','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/6/23/mathml/M2">View MathML</a> ~t(nx + ny - 2) where sx = 4.3695 and sy = 6.9421 are the standard deviations of x and y; <a onClick="popup('http://www.biomedcentral.com/1752-0509/6/23/mathml/M3','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/6/23/mathml/M3">View MathML</a> and ȳ = 4.9804 are the means of x and y; nx and ny are the numbers of data of x and y. The p-value for the null hypothesis <a onClick="popup('http://www.biomedcentral.com/1752-0509/6/23/mathml/M4','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/6/23/mathml/M4">View MathML</a> is less than 0.0001 which is much less than the significance level α = 0.05. Therefore, the null hypothesis H0 should be rejected and the alternative is accepted.

Discussion

The present study unraveled the direction of neural network development to facilitate a relatively high level of SBA. Thus, the CR range of a mature cultured neural network may represent a delicate design and not the result of random selection. In addition, such biological interpretation of the optimal CR may be further applicable to in vivo situations, since the distribution of cell types in ex vivo networks is often similar to those found in vivo [25,36]. Some evidence indicates that neural networks first develop toward certain connective structures and then form specific functions by adjusting their synaptic efficacies according to the external stimulus [37]. Our simulations suggest that neurons may connect with each other at a 10-30% CR to achieve the highest possibility of SBA occurrence in the early stage, and then, based on such an optimal CR range, the constructed networks further recruit and control SBA by chemically adjusting their synaptic efficacies.

We showed that our main results are quite robust to variations of network scales, network topological properties, and simulation parameters. We carried out simulations (see Simulation protocols) for a variety of neural networks with 10-300 nodes and found that the mean value of all the optimal CRs is 13% (with a standard deviation 0.0181) which lies within the evolutionarily selected range of CR (10-30%). In addition, note that the networks used for simulations in the early part (1,000 networks were constructed for each CR and ER) were based on random connections and therefore various possible topological structures were already taken into account. So, we confirmed that our results hold regardless of particular connective forms. We have also investigated the possible influences by perturbation of parameters {τs, τm, Δax, τ} (over a range from -95% to 250% perturbations of the original simulation parameters). For instance, the mean of all the optimal CRs for 10-node networks with the parameter perturbations was 18.67% (with a standard deviation 0.0183). In this way, we could also confirm that our results, the evolutionarily selected range of CR of 10-30%, still hold against the parameter perturbations.

Conclusions

In this study, we investigated the underlying cause of the evolutionarily selected CRs of neural networks. Artificial pulsed neural network simulation has shown that an optimal CR range (10-30%) maximizes the occurrence of synchronized bursting behaviors (when ER = 0.9), which is consistent with previous ex vivo experimental observation, in which the CRs of cultured cortical networks consistently lie in a range of 10-30% with an ER of 80-90%.

Employing time-series data from multi-electrode array experiments, we identified some APFL motifs in cultured cortical networks of E18 Sprague-Dawley rats [2]. To further unravel the crucial role of the APFL motifs, we investigated the relationship between specific network structures (network motifs) and network behaviors in artificial pulsed neural networks. This study can readily be used to capture the fundamental dynamical characteristics of cultured neural networks. We found that the existence of APFL motifs is a necessary condition of SBA, not only for small-scale networks but also for large-scale networks. To recruit a high level of SBA, networks must have an optimal number of APFL motifs. Therefore, we infer that the formation of the appropriate number of APFLs is related to the maximal occurrence of SBA, whereas the optimal CR is only a necessary condition to achieve the required APFLs.

Furthermore, we investigated the distribution of each type of APFL motif (2, 3, or 4-node) at different SBA levels. In both 12-node and 60-node networks, the 2-node APFL motif dominated among APFL motifs at high SBA levels. More importantly, the contribution of each type of APFL motif to SBA was demonstrated by comparing the inhibitory effect of each APFL motif against SBA. For large-scale networks, the exclusion of 2-node APFLs almost fully prohibits the occurrence of SBA, implying that compared to other APFL motifs, 2-node APFL motifs may be crucial for neural networks to produce SBA.

Methods

Definitions of network motifs and feedback loops

A network motif is defined as an enriched sub-network pattern in complex networks that occurs more frequently than in randomized networks [38-41]. Here, this concept was extended to a more general definition. A motif refers to any sub-network with a particular structure. To relate the structure of various distinct network motifs and their dynamic behaviors, a range of different network structures can be considered, and their correlations to specific dynamic behaviors investigated. This paper focused on synchronized bursting activities.

A feedback loop (FBL) is defined as a network motif composed of network nodes (neurons) and closed directed paths (synapses). The example network shown in Figure 5a contains five FBLs (Figure 5b). FBL1, FBL2 and FBL3 are positive feedback loops (PFLs); FBL4 and FBL5 are NFLs. The sign of a feedback loop is determined by (-1), where q is the total number of inhibitory interactions contained in the loop. Moreover, if a positive feedback loop contains only positive interactions like FBL1 and FBL2, it is referred to as an all-positive-interaction feedback loop (APFL).

thumbnailFigure 5. Illustration of feedback loops and APFLs. (a) An example network composed of five different feedback loops. An interaction between nodes is represented by arrows to denote excitatory regulation and blunt lines for inhibitory regulation. The sub-networks colored in red denote the two APFL motifs. (b) Five feedback loops (FBL) contained in the network shown in (a).

Pulsed neural networks

To infer the relationship between a type of feedback motif and its network behaviors, typical network motifs were constructed based on pulsed neural networks, and their network responses to randomly assigned initial states and simulation parameters were observed. The pulsed neural networks, also called the third generation of artificial neural networks, are based on spiking neurons, or "integrate-and-fire" neurons [22,23]. These neurons utilize recent insights from neurophysiology, specifically the use of temporal coding to pass information between neurons [11,42,43], which closely mimic realistic communication between neurons. Therefore, pulsed neural networks are commonly applied to study the properties of neural networks.

For a spiking neuron i, the membrane voltage can be denoted by a state variable xi. Once xi reaches a threshold δ, the neuron is fired; the moment of crossing the threshold is represented by a firing time <a onClick="popup('http://www.biomedcentral.com/1752-0509/6/23/mathml/M5','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/6/23/mathml/M5">View MathML</a>. The set of all firing times of neuron i, commonly called a spike train, is described as

<a onClick="popup('http://www.biomedcentral.com/1752-0509/6/23/mathml/M6','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/6/23/mathml/M6">View MathML</a>

(1)

where <a onClick="popup('http://www.biomedcentral.com/1752-0509/6/23/mathml/M7','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/6/23/mathml/M7">View MathML</a> is the most recent spike before the current time t. Two different processes contribute to the value of xi. The first contribution is a negative-value function <a onClick="popup('http://www.biomedcentral.com/1752-0509/6/23/mathml/M8','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/6/23/mathml/M8">View MathML</a> indicating an immediate "reset" after each firing time in Φi. In the biological context, Ψi is used to account for neuronal refractoriness. The second contribution is inputs from pre-synaptic neurons j Λi where

<a onClick="popup('http://www.biomedcentral.com/1752-0509/6/23/mathml/M9','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/6/23/mathml/M9">View MathML</a>

(2)

A pre-synaptic spike at time <a onClick="popup('http://www.biomedcentral.com/1752-0509/6/23/mathml/M10','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/6/23/mathml/M10">View MathML</a> increases (or decreases) the state xi of post-synaptic neuron i for <a onClick="popup('http://www.biomedcentral.com/1752-0509/6/23/mathml/M11','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/6/23/mathml/M11">View MathML</a> by summing up a weighted kernel function as <a onClick="popup('http://www.biomedcentral.com/1752-0509/6/23/mathml/M12','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/6/23/mathml/M12">View MathML</a> The signs can be reflected in synaptic efficacy, wij, using wij > 0 for excitatory synapses and wij < 0 for inhibitory synapses. The kernel εij describes the response of xi due to pre-synaptic potentials at <a onClick="popup('http://www.biomedcentral.com/1752-0509/6/23/mathml/M10','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/6/23/mathml/M10">View MathML</a>, which can be viewed as a combined effect of the axonal transmission and membrane transmission properties of neurons. Therefore, the state of neuron i at current time t is given by a linear superposition of the two main previously mentioned contributions,

<a onClick="popup('http://www.biomedcentral.com/1752-0509/6/23/mathml/M13','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/6/23/mathml/M13">View MathML</a>

(3)

The models described by (1)-(3) are referred to as SRMs [22]. They, together with the connectivity topology of neural networks, form a simple mechanism of simulating biological neural networks. Frequently, noise was introduced into the SRM by adding an effect of a stochastic current <a onClick="popup('http://www.biomedcentral.com/1752-0509/6/23/mathml/M14','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/6/23/mathml/M14">View MathML</a> to the right-hand side of (3). Then (3) can be altered to

<a onClick="popup('http://www.biomedcentral.com/1752-0509/6/23/mathml/M15','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/6/23/mathml/M15">View MathML</a>

(4)

where the kernel function e(s) mimics the dynamic from the local noise current stimulation to the membrane voltage of neuron i. As usual, several typical mathematical formulations were adopted (illustrated in Figure 6) to describe refractoriness Ψi, post-synaptic potential εij, and membrane dynamics ei. For instance, let

<a onClick="popup('http://www.biomedcentral.com/1752-0509/6/23/mathml/M16','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/6/23/mathml/M16">View MathML</a>

(5)

thumbnailFigure 6. The functions describing the dynamics of a spike neuron. (a) The kernel εij(t) describing the response of xi(t) caused by a pre-synaptic spike at t = 0, with Δax = 50 ms, τs = 3.5 ms, and τm = 8 ms. (b) The function ψi (t) reflecting refractoriness after a spike emitted at t = 0, with δ = 0.1 and τ = 40 ms. (c) The kernel ei(t) representing the dynamics from local current stimulation to the membrane voltage of a neuron. The time constant is identical to τm in (a). (d) The membrane voltage xi (t) firing at time <a onClick="popup('http://www.biomedcentral.com/1752-0509/6/23/mathml/M5','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/6/23/mathml/M5">View MathML</a> when it reaches a threshold voltage δ = 0.1. After firing, it is reset by the function ψi(t) and then re-accumulates via the pre-synaptic spike inputs <a onClick="popup('http://www.biomedcentral.com/1752-0509/6/23/mathml/M12','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/6/23/mathml/M12">View MathML</a>.

where τs and τm are time constants describing axonal transmission dynamics and membrane dynamics, respectively, and Δax is axonal transmission delay. H(t - Δax) is the Heaviside step function which vanishes for t ≤ Δax, and set t > Δax equal to 1.

One typical membrane voltage reset function is

<a onClick="popup('http://www.biomedcentral.com/1752-0509/6/23/mathml/M17','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/6/23/mathml/M17">View MathML</a>

(6)

where Trefractoryis the absolute refractory period of neuron i. During such a period, the neuron cannot be fired regardless of the membrane voltage. For the membrane dynamical function, the following equation can be used:

<a onClick="popup('http://www.biomedcentral.com/1752-0509/6/23/mathml/M18','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/6/23/mathml/M18">View MathML</a>

(7)

Networks of different sizes can exhibit similar network behaviors (or dynamics) if their neurons are supplied with the same average inputs [22]. To make networks of two sizes comparable with respect to the same average input of each neuron, the networks must have weight scopes scaled by the number of network nodes. For example, take an n1-node network as a nominal case with an allowed weight scope of [-Wmax, Wmax]. Then the weight scope of an n2-node network should be assigned as <a onClick="popup('http://www.biomedcentral.com/1752-0509/6/23/mathml/M19','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/6/23/mathml/M19">View MathML</a>. In practice, the smaller network (n1 = 12) was set as the nominal network. Thus, the weight scope of the other network (n2 = 60) was scaled by the factor (60 - 1)/(12 - 1) ≈ 5.

The typical behaviors of spike neural networks

Four typical spontaneous network behaviors appeared during the simulations: transient response activity (TRA), SBA, asynchronized bursting activities (ASBA), and HEA [44-47]. TRA is a non-sustained firing phenomenon, while both SBA and ASBA demonstrate regularly separated clusters of spikes during their full durations. The two cases can be differentiated using a synchrony index (SI) [48,49]. SI is defined based on a cross-correlation coefficient, which is often used to quantify the temporal relationship of a pair of neurons. Suppose the spike trains of the two cells in a duration of T seconds are denoted by x(t) and y(t) (0 < t < T). Their discrete-time forms x(k) and y(k) can be obtained by dividing T into n bins (T/bin width; k = 1 . . . n) and then counting the number of spikes in each bin width. In our simulations, T = 2 s, and bin width = 20 ms. The simulations showed that our results are quite robust with respect to arbitrary choices of T and bin width. The cross-correlation coefficient r between x(k) and y(k) is calculated as follows:

<a onClick="popup('http://www.biomedcentral.com/1752-0509/6/23/mathml/M20','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/6/23/mathml/M20">View MathML</a>

(8)

where <a onClick="popup('http://www.biomedcentral.com/1752-0509/6/23/mathml/M21','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/6/23/mathml/M21">View MathML</a>, <a onClick="popup('http://www.biomedcentral.com/1752-0509/6/23/mathml/M22','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/6/23/mathml/M22">View MathML</a>, and <a onClick="popup('http://www.biomedcentral.com/1752-0509/6/23/mathml/M23','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/6/23/mathml/M23">View MathML</a>.

By its definition, r is a value in the range of [0, 1]. If r exceeds a threshold (0.7 was used in our simulations), the two spike trains are considered synchronized with each other. Otherwise, they are considered asynchronized. In the case of HEA, all nodes permanently fire with an interval of absolute refractory period Trefractory. Notably, hyper-excited spike trains also have a high SI, which is quite similar to SBA. However, they fire as closely as possible without any perceivable intermittence during their final steady state. Figure 7 illustrates these four typical activities using a simple 2-node PFL motif. For a larger-scale network with n nodes, k-channel SBAs (2 ≤ k n and k ) can be calculated by the same criteria. For simplicity, only 2-channel SBA was selected in this study to represent the level of k-channel SBAs, considering that measurements greater than 2-channel do not change the final results and major conclusions of the study. Similarly, 2-channel HEA was selected in to represent the level of k-channel HEAs either.

thumbnailFigure 7. Four possible network behaviors for a 2-node PFL motif with various synaptic efficacies. The initial state of the network is taken as (1,1). (a) The scheme of a 2-node PFL motif. (b) Transient response activity. (c) Synchronized bursting activity. (d) Asynchronized bursting activity. (e) Hyper-excitable activity.

Simulation protocols

Simulations were carried out using SRMs for the randomly connected networks where all neurons were assumed to have an identical parameter set {τs, τm, Δax, τ, δ}, and the synaptic efficacies wij were randomly chosen from a uniform distribution [-1, 1] (see the section of 'Pulsed neural networks' for further details on SRM). For a randomly generated network with n neurons and m synaptic connections (suppose that me denotes excitatory connections and m-me represents inhibitory connections), the CR is defined as the percentage of the number of existing synaptic connections of the network divided by that of the fully connected network with node number n, i.e., <a onClick="popup('http://www.biomedcentral.com/1752-0509/6/23/mathml/M24','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/6/23/mathml/M24">View MathML</a>. Here, only a single connection between two different neurons is allowed for simplification. Hence, the possibility of self-connection of one neuron and multiple connections between any pair of neurons are excluded. The ER is referred to as a quotient of excitatory synaptic number over the total number of synaptic connections, i.e., <a onClick="popup('http://www.biomedcentral.com/1752-0509/6/23/mathml/M25','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/6/23/mathml/M25">View MathML</a>.

To investigate the relationships among CR, ER, and the occurrence of SBA, simulations with both synaptic efficacies and network structures randomly perturbed were carried out using different combinations of CR and ER. Classification of four typical network behaviors (TRA, SBA, ASBA, and HEA) can be found in the section entitled "The typical behaviors of spike neural networks". For each ratio pair (CR, ER), 1,000 randomly connected artificial neural networks were constructed, and simulations based on these networks were carried out. For each constructed network (corresponding to one simulation), the total number of APFL motifs (2, 3, and 4-node) and the occurrences of k-channel SBA and HEA (2 ≤ k n and k ) were investigated to show the correlations between APFL number and such typical behaviors. Our simulations demonstrated that the change of k or inclusion of more (greater than 4) node motifs had no effect on evaluating the dynamical characteristics of networks. In practice, only the 2-channel SBA and 2-channel HEA were evaluated and used to measure the levels of synchronization and hyper-excitation in networks. The correlation between the occurrence of 2-channel SBA or 2-channel HEA and network topological characteristics was obtained from the statistics of the occurrence of such behaviors with various CRs, ERs, and randomly assigned synaptic efficacies.

Abbreviations

APFL: All-positive-interaction feedback loop; TRA: Transient response activities; SBA: Synchronized bursting activities; ASBA: Asynchronous bursting activities; HEA: Hyper-excitable activities; CR: Connective ratio; ER: Excitatory ratio; SI: Synchrony index; SRM: Spike response model.

Competing interests

The authors declare that they have no competing interests.

Authors' contributions

KHC designed the study, CYD performed the simulations, CYD and KHC analyzed the data, and CYD and KHC wrote the paper. All authors read and approved the final manuscript.

Acknowledgements

We thank Dongkwan Shin and Chaoxuan Dong for their critical reading of this paper. This work was supported by the National Re-search Foundation of Korea (NRF) grants funded by the Korea Government, the Ministry of Education, Science & Technology (MEST) (2009-0086964, 2010-0017662, and 2011-0006314). It was also supported by the WCU (World Class University) program (R32-2008-000-10218-0) through the NRF funded by MEST.

References

  1. Aradi I, Maccaferri G: Cell type-specific synaptic dynamics of synchronized bursting in the juvenile CA3 rat hippocampus.

    J Neurosci 2004, 24:9681. PubMed Abstract | Publisher Full Text OpenURL

  2. Dong C-Y, Lim J, Nam Y, Cho K-H: Systematic analysis of synchronized oscillatory neuronal networks reveals an enrichment for coupled direct and indirect feedback motifs.

    Bioinformatics 2009, 25:1680-1685. PubMed Abstract | Publisher Full Text OpenURL

  3. Klemm K, Bornholdt S: Topology of biological networks and reliability of information processing.

    Proc Natl Acad Sci 2005, 102:18414-18419. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  4. Kamioka H, Maeda E, Jimbo Y, Robinson HPC, Kawana A: Spontaneous periodic synchronized bursting during formation of mature patterns of connections in cortical cultures.

    Neurosci Lett 1996, 206:109-112. PubMed Abstract | Publisher Full Text OpenURL

  5. Gray CM, Singer W: Stimulus-specific neuronal oscillations in orientation columns of cat visual cortex.

    Proc Natl Acad Sci 1989, 86:1698-1702. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  6. Bragin A, Jando G, Nadasdy Z, Hetke J, Wise K, Buzsaki G: Gamma (40-100 Hz) oscillation in the hippocampus of the behaving rat.

    J Neurosci 1995, 15:47-60. PubMed Abstract | Publisher Full Text OpenURL

  7. Traub R, Contreras D, Cunningham M, Murray H, LeBeau F, Roopun A, Bibbig A, Wilent W, Higley M, Whittington M: Single-column thalamocortical network model exhibiting gamma oscillations, sleep spindles, and epileptogenic bursts.

    J Neurophysiol 2005, 93:2194-2232. PubMed Abstract | Publisher Full Text OpenURL

  8. Womelsdorf T, Fries P: The role of neuronal synchronization in selective attention.

    Curr Opin Neurobiol 2007, 17:154-160. PubMed Abstract | Publisher Full Text OpenURL

  9. Siegel M, Donner TH, Oostenveld R, Fries P, Engel AK: Neuronal synchronization along the dorsal visual pathway reflects the focus of spatial attention.

    Neuron 2008, 60:709-719. PubMed Abstract | Publisher Full Text OpenURL

  10. Engel AK, Fries P, Singer W: Dynamic predictions: oscillations and synchrony in top-down processing.

    Nat Rev Neurosci 2001, 2:704-716. PubMed Abstract | Publisher Full Text OpenURL

  11. Riehle A, Grun S, Diesmann M, Aertsen A: Spike synchronization and rate modulation differentially involved in motor cortical function.

    Science 1950, 1997:278. OpenURL

  12. Gray CM, Konig P, Engel AK, Singer W: Oscillatory responses in cat visual cortex exhibit inter-columnar synchronization which reflects global stimulus properties.

    Nature 1989, 338:334-337. PubMed Abstract | Publisher Full Text OpenURL

  13. deCharms RC, Merzenich MM: Primary cortical representation of sounds by the coordination of action-potential timing.

    Nature 1996, 381:610-613. PubMed Abstract | Publisher Full Text OpenURL

  14. Kim JR, Shin D, Jung SH, Heslop-Harrison P, Cho KH: A design principle underlying the synchronization of oscillations in cellular systems.

    J Cell Sci 2010, 123:537. PubMed Abstract | Publisher Full Text OpenURL

  15. Muramoto K, Ichikawa M, Kawahara M, Kobayashi K, Kuroda Y: Frequency of synchronous oscillations of neuronal activity increases during development and is correlated to the number of synapses in cultured cortical neuron networks.

    Neurosci Lett 1993, 163:163-165. PubMed Abstract | Publisher Full Text OpenURL

  16. Nakanishi K, Nakanishi M, Kukita F: Dual intracellular recording of neocortical neurons in a neuron-glia co-culture system.

    Brain Res Protoc 1999, 4:105-114. Publisher Full Text OpenURL

  17. Jimbo Y, Tateno T, Robinson HPC: Simultaneous induction of pathway-specific Potentiation and depression in networks of cortical neurons.

    Biophys J 1999, 76:670-678. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  18. Marom S, Shahaf G: Development, learning and memory in large random networks of cortical neurons: lessons beyond anatomy.

    Q Rev Biophys 2002, 35:63-87. PubMed Abstract OpenURL

  19. Fields RD, Nelson PG: Activity-dependent development of the vertebrate nervous system.

    Int Rev Neurobiol 1992, 34:133-214. PubMed Abstract OpenURL

  20. Segev R, Shapira Y, Benveniste M, Ben-Jacob E: Observations and modeling of synchronized bursting in two-dimensional neural networks.

    Phys Rev E 2001, 64:11920. OpenURL

  21. Izhikevich EM, Gally JA, Edelman GM: Spike-timing dynamics of neuronal groups.

    Cereb Cortex 2004, 14:933. PubMed Abstract | Publisher Full Text OpenURL

  22. Maass W, Bishop CM: Pulsed Neural Networks. Cambridge: MIT Press; 1999. OpenURL

  23. Maass W: Networks of spiking neurons: the third generation of neural network models.

    Neural Netw 1997, 10:1659-1671. Publisher Full Text OpenURL

  24. Eckenstein F, Thoenen H: Cholinergic neurons in the rat cerebral cortex demonstrated by immunohistochemical localization of choline acetyltransferase.

    Neurosci Lett 1983, 36:211. PubMed Abstract | Publisher Full Text OpenURL

  25. Huettner JE, Baughman RW: Primary culture of identified neurons from the visual cortex of postnatal rats.

    J Neurosci 1986, 6:3044-3060. PubMed Abstract | Publisher Full Text OpenURL

  26. Lodato I, Boccaletti S, Latora V: Synchronization properties of network motifs.

    EPL (Europhysics Letters) 2007, 78:28001. Publisher Full Text OpenURL

  27. DHuys O, Vicente R, Erneux T, Danckaert J, Fischer I: Synchronization properties of network motifs: Influence of coupling delay and symmetry.

    Chaos 2008, 18:037116. PubMed Abstract | Publisher Full Text OpenURL

  28. Chatterjee N, Sinha S: Understanding the mind of a worm: hierarchical network structure underlying nervous system function in C. elegans.

    Prog Brain Res 2007, 168:145-153. OpenURL

  29. Hobert O: Specification of the nervous system.

    WormBook: the online review of C elegans biology 2005, 12:1. OpenURL

  30. Watts DJ, Strogatz SH: Collective dynamics of'small-world' networks.

    Nature(London) 1998, 393:440-442. Publisher Full Text OpenURL

  31. Varshney LR, Chen BL, Paniagua E, Hall DH, Chklovskii DB: Structural Properties of the Caenorhabditis elegans Neuronal Network.

    PLoS Comput Biol 2011, 7:e1001066. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  32. Ambros V, Horvitz HR: Heterochronic mutants of the nematode Caenorhabditis elegans.

    Science 1984, 226:409. PubMed Abstract | Publisher Full Text OpenURL

  33. Trent C, Tsung N, Horvitz HR: Egg-laying defective mutants of the nematode Caenorhabditis elegans.

    Genetics 1983, 104:619. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  34. Gerstner W, Kistler W: Spiking Neuron Models: An Introduction. NY: Cambridge University Press New York; 2002. OpenURL

  35. Kim SH, Cohen AS: On the Behrens-Fisher problem: a review.

    J Educ and Behav Stat 1998, 23:356-377. OpenURL

  36. Nakanishi K, Kukita F: Intracellular [Cl-] modulates synchronous electrical activity in rat neocortical neurons in culture by way of GABAergic inputs.

    Brain Res 2000, 863:192-204. PubMed Abstract | Publisher Full Text OpenURL

  37. Hagan MT, Demuth HB, Beale MH: Neural network design. MA: PWS Boston; 1996. OpenURL

  38. Kim J-R, Yoon Y, Cho K-H: Coupled feedback loops form dynamic motifs of cellular networks.

    Biophys J 2008, 94:359-365. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  39. Kim TH, Kim J, Heslop-Harrison P, Cho KH: Evolutionary design principles and functional characteristics based on kingdom-specific network motifs.

    Bioinformatics 2010, 27:245. PubMed Abstract | Publisher Full Text OpenURL

  40. Kwon Y-K, Cho K-H: Quantitative analysis of robustness and fragility in biological networks based on feedback dynamics.

    Bioinformatics 2008, 24:987-994. PubMed Abstract | Publisher Full Text OpenURL

  41. Milo R, Shen-Orr S, Itzkovitz S, Kashtan N, Chklovskii D, Alon U: Network motifs: simple building blocks of complex networks.

    Science 2002, 298:824-827. PubMed Abstract | Publisher Full Text OpenURL

  42. Buracas GT, Zador AM, DeWeese MR, Albright TD: Efficient discrimination of temporal patterns by motion-sensitive neurons in primate visual cortex.

    Neuron 1998, 20:959-969. PubMed Abstract | Publisher Full Text OpenURL

  43. Thorpe S, Fize D, Marlot C: Speed of processing in the human visual system.

    Nature 1996, 381:520-522. PubMed Abstract | Publisher Full Text OpenURL

  44. Mehring C, Hehl U, Kubo M, Diesmann M, Aertsen A: Activity dynamics and propagation of synchronous spiking in locally connected random networks.

    Biol Cybern 2003, 88:395. PubMed Abstract | Publisher Full Text OpenURL

  45. Brunel N: Dynamics of networks of randomly connected excitatory and inhibitory spiking neurons.

    J Physiol Paris 2000, 94:445. PubMed Abstract | Publisher Full Text OpenURL

  46. Silberberg G, Bethge M, Markram H, Pawelzik K, Tsodyks M: Dynamics of population rate codes in ensembles of neocortical neurons.

    J Neurophysiol 2004, 91:704. PubMed Abstract | Publisher Full Text OpenURL

  47. Bonneau R: Learning biological networks: from modules to dynamics.

    Nat Chem Biol 2008, 4:658-664. PubMed Abstract | Publisher Full Text OpenURL

  48. Gerstein GL, Kiang NYS: An approach to the quantitative analysis of electrophysiological data from single neurons.

    Biophys J 1960, 1:15-28. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  49. Wylie DR, De Zeeuw CI, Simpson JI: Temporal relations of the complex spike activity of Purkinje cell pairs in the vestibulocerebellum of rabbits.

    J Neurosci 1995, 15:2875. PubMed Abstract | Publisher Full Text OpenURL