We investigate how network structure can influence the tendency for a neuronal network to synchronize, or its synchronizability, independent of the dynamical model for each neuron. The synchrony analysis takes advantage of the framework of second order networks, which defines four second order connectivity statistics based on the relative frequency of two-connection network motifs. The analysis identifies two of these statistics, convergent connections, and chain connections, as highly influencing the synchrony. Simulations verify that synchrony decreases with the frequency of convergent connections and increases with the frequency of chain connections. These trends persist with simulations of multiple models for the neuron dynamics and for different types of networks. Surprisingly, divergent connections, which determine the fraction of shared inputs, do not strongly influence the synchrony. The critical role of chains, rather than divergent connections, in influencing synchrony can be explained by their increasing the effective coupling strength. The decrease of synchrony with convergent connections is primarily due to the resulting heterogeneity in firing rates.

High frequency deep-brain stimulation of the subthalamic nucleus (deep brain stimulation, DBS) relieves many of the symptoms of Parkinson's disease in humans and animal models. Although the treatment has seen widespread use, its therapeutic mechanism remains paradoxical. The subthalamic nucleus is excitatory, so its stimulation at rates higher than its normal firing rate should worsen the disease by increasing subthalamic excitation of the globus pallidus. The therapeutic effectiveness of DBS is also frequency and intensity sensitive, and the stimulation must be periodic; aperiodic stimulation at the same mean rate is ineffective. These requirements are not adequately explained by existing models, whether based on firing rate changes or on reduced bursting. Here we report modeling studies suggesting that high frequency periodic excitation of the subthalamic nucleus may act by desynchronizing the firing of neurons in the globus pallidus, rather than by changing the firing rate or pattern of individual cells. Globus pallidus neurons are normally desynchronized, but their activity becomes correlated in Parkinson's disease. Periodic stimulation may induce chaotic desynchronization by interacting with the intrinsic oscillatory mechanism of globus pallidus neurons. Our modeling results suggest a mechanism of action of DBS and a pathophysiology of Parkinsonism in which synchrony, rather than firing rate, is the critical pathological feature.

Phase resetting curves (PRCs) provide a measure of the sensitivity of oscillators to perturbations. In a noisy environment, these curves are themselves very noisy. Using perturbation theory, we compute the mean and the variance for PRCs for arbitrary limit cycle oscillators when the noise is small. Phase resetting curves and phase dependent variance are fit to experimental data and the variance is computed using an ad-hoc method. The theoretical curves of this phase dependent method match both simulations and experimental data significantly better than an ad-hoc method. A dual cell network simulation is compared to predictions using the analytical phase dependent variance estimation presented in this paper. We also discuss how entrainment of a neuron to a periodic pulse depends on the noise amplitude.

A tonic-clonic seizure transitions from high frequency asynchronous activity to low frequency coherent oscillations, yet the mechanism of transition remains unknown. We propose a shift in network synchrony due to changes in cellular response. Here we use phase-response curves (PRC) from Morris-Lecar (M-L) model neurons with synaptic depression and gradually decrease input current to cells within a network simulation. This method effectively decreases firing rates resulting in a shift to greater network synchrony illustrating a possible mechanism of the transition phenomenon. PRCs are measured from the M-L conductance based model cell with a range of input currents within the limit cycle. A large network of 3000 excitatory neurons is simulated with a network topology generated from second-order statistics which allows a range of population synchrony. The population synchrony of the oscillating cells is measured with the Kuramoto order parameter, which reveals a transition from tonic to clonic phase exhibited by our model network. The cellular response shift mechanism for the tonic-clonic seizure transition reproduces the population behavior closely when compared to EEG data.

No abstract