Pacemaker and network mechanisms of rhythm generation: cooperation and competition

The origin of rhythmic activity in brain circuits and CPG-like motor networks is still not fully understood. The main unsolved questions are (i) What are the respective roles of intrinsic bursting and network based dynamics in systems of coupled heterogeneous, intrinsically complex, even chaotic, neurons? (ii) What are the mechanisms underlying the coexistence of robustness and flexibility in the observed rhythmic spatio-temporal patterns? One common view is that particular bursting neurons provide the rhythmogenic component while the connections between different neurons are responsible for the regularisation and synchronisation of groups of neurons and for specific phase relationships in multi-phasic patterns. We have examined the spatio-temporal rhythmic patterns in computer-simulated motif networks of H-H neurons connected by slow inhibitory synapses with a non-symmetric pattern of coupling strengths. We demonstrate that the interplay between intrinsic and network dynamics features either cooperation or competition, depending on three basic control parameters identified in our model: the shape of intrinsic bursts, the strength of the coupling and its degree of asymmetry. The cooperation of intrinsic dynamics and network mechanisms is shown to correlate with bistability, i.e., the coexistence of two different attractors in the phase space of the system corresponding to different rhythmic spatio-temporal patterns. Conversely, if the network mechanism of rhythmogenesis dominates, monostability is observed with a typical pattern of winnerless competition between neurons. We analyse bifurcations between the two regimes and demonstrate how they provide robustness and flexibility to the network performance.     

Coherency and connectivity in oscillating neural networks: linear partialization analysis

 This paper studies the relation between the functional synaptic connections between two artificial neural networks and the correlation of their spiking activities. The model neurons had realistic non-oscillatory dynamic properties and the networks showed oscillatory behavior as a result of their internal synaptic connectivity. We found that both excitation and inhibition cause phase locking of the oscillating activities. When the two networks excite each other the oscillations synchronize with zero phase lag, whereas mutual inhibition between the networks resulted in an anti-phase (half period phase difference) synchronization. Correlations between the activities of the two networks can also be caused by correlated external inputs driving the systems (common input). Our analysis shows that when the networks exhibit oscillatory behavior and the rate of the common input is smaller than a characteristic network oscillator frequency, the cross-correlation functions between the activities of two systems still carry information about the mutual synaptic connectivity. This information can be retrieved with linear partialization, removing the influence of the common input. We further explored the network responses to periodic external input. We found that when the input is of a frequency smaller than a certain threshold, the network responds with bursts at the same frequency as the input. Above the threshold, the network responds with a fraction of the input frequency. This frequency threshold, characterizing the oscillatory properties of the network, is also found to determine the limit to which linear partialization works. Received: 20 October 1995 / Accepted in revised form: 20 May 1996     

Linking dynamics of the inhibitory network to the input structure

Networks of inhibitory interneurons are found in many distinct classes of biological systems. Inhibitory interneurons govern the dynamics of principal cells and are likely to be critically involved in the coding of information. In this theoretical study, we describe the dynamics of a generic inhibitory network in terms of low-dimensional, simplified rate models. We study the relationship between the structure of external input applied to the network and the patterns of activity arising in response to that stimulation. We found that even a minimal inhibitory network can generate a great diversity of spatio-temporal patterning including complex bursting regimes with non-trivial ratios of burst firing. Despite the complexity of these dynamics, the network’s response patterns can be predicted from the rankings of the magnitudes of external inputs to the inhibitory neurons. This type of invariant dynamics is robust to noise and stable in densely connected networks with strong inhibitory coupling. Our study predicts that the response dynamics generated by an inhibitory network may provide critical insights about the temporal structure of the sensory input it receives.     

Parallel processing by a homogeneous group of coupled model neurons can enhance, reduce and generate signal correlations

 Correlated activities have been proposed as correlates of flexible association and assembly coding. We addressed the basic question of how signal correlations on parallel pathways are enhanced, reduced and generated by homogeneous groups of coupled neurons, and how this depends on the input activities and their interactions with internal coupling processes. For this we simulated a fully connected group of identical impulse-coded neurons with dynamic input and threshold processes and additive or multiplicative lateral coupling. Input signals were Gaussian white noise (GWN), completely independent or partially correlated on a subgroup of the parallel inputs. We show that in states of high average spike rates input-output correlations were weak while the network could generate correlated activities of stochastic, oscillatory and rhythmic bursting types depending exclusively on lateral coupling strength. In states of low average spike rates input-output correlations were high and the network could effectively enhance or reduce differences in spatial correlation applied to its parallel inputs. The correlation differences were more pronounced with multiplicative lateral coupling than with the additive interactions commonly used. As the different modes of correlation processing emerged already by global changes in the average spike rate and lateral coupling strength, we assume that in real cortical circuits changes in correlational processing may also be induced by unspecific modulations of activation and lateral coupling. Received: 11 December 1995 / Accepted in revised form: 29 November 1996     

Synchronized oscillations in the visual cortex — a synergetic model

 We present an oscillator network model for the synchronization of oscillatory neuronal activity underlying visual processing. The single neuron is modeled by means of a limit cycle oscillator with an eigenfrequency corresponding to visual stimulation. The eigenfrequency may be time dependent. The mutual coupling strengths are unsymmetrical and activity dependent, and they scatter within the network. Synchronized clusters (groups) of neurons emerge in the network due to the visual stimulation. The different clusters correspond to different visual stimuli. There is no limitation of the number of stimuli. Distinct clusters do not perturb each other, although the coupling strength between all model neurons is of the same order of magnitude. Our analysis is not restricted to weak coupling strength. The scatter of the couplings causes shifts of the cluster frequencies. The model’s behavior is compared with the experimental findings. The coupling mechanism is extended in order to model the influence of bicucullin upon the neural network. We additionally investigate repulsive couplings, which lead to constant phase differences between clusters of the same frequency. Finally, we consider the problem of selective attention from the viewpoint of our model. Received: 15 February 1995/Accepted in revised form: 18 July 1995     

Synchronizing a multicellular system by external input: an artificial control strategy

MOTIVATION: Although there are significant advances on elucidating the collective behaviors on biological organisms in recent years, the essential mechanisms by which the collective rhythms arise remain to be fully understood, and further how to synchronize multicellular networks by artificial control strategy has not yet been well explored. RESULTS: A control strategy is developed to synchronize gene regulatory networks in a multicellular system when spontaneous synchronization cannot be achieved. We first construct an impulsive control system to model the process of periodically injecting coupling substances with constant or random impulsive control amounts into the common extracellular medium, and further study its effects on the dynamics of individual cells. We derive the threshold of synchronization induced by the periodic substance input. Therefore, we can synchronize the multicellular network to a specific collective behavior by changing the frequency and amplitude of the periodic stimuli. Moreover, a two-stage scheme is proposed to facilitate the synchronization in this paper. We show that the presence of the external input may also initiate different dynamics. The multicellular network of coupled repressilators is used to show the effectiveness of the proposed method. The results not only provide a perspective to understand the interactions between external stimuli and intrinsic physiological rhythms, but also may lead to development of realistic artificial control strategy and medical therapy. AVAILABILITY: CONTACT: aihara@sat.t.u-tokyo.ac.jp.     

Dynamics of spiking neurons connected by both inhibitory and electrical coupling

We study the dynamics of a pair of intrinsically oscillating leaky integrate-and-fire neurons (identical and noise-free) connected by combinations of electrical and inhibitory coupling. We use the theory of weakly coupled oscillators to examine how synchronization patterns are influenced by cellular properties (intrinsic frequency and the strength of spikes) and coupling parameters (speed of synapses and coupling strengths). We find that, when inhibitory synapses are fast and the electrotonic effect of the suprathreshold portion of the spike is large, increasing the strength of weak electrical coupling promotes synchrony. Conversely, when inhibitory synapses are slow and the electrotonic effect of the suprathreshold portion of the spike is small, increasing the strength of weak electrical coupling promotes antisynchrony (see Fig. 10). Furthermore, our results indicate that, given a fixed total coupling strength, either electrical coupling alone or inhibition alone is better at enhancing neural synchrony than a combination of electrical and inhibitory coupling. We also show that these results extend to moderate coupling strengths.     

Synchronized oscillations in the visual cortex — a synergetic model

We present an oscillator network model for the synchronization of oscillatory neuronal activity underlying visual processing. The single neuron is modeled by means of a limit cycle oscillator with an eigenfrequency corresponding to visual stimulation. The eigenfrequency may be time dependent. The mutual coupling strengths are unsymmetrical and activity dependent, and they scatter within the network. Synchronized clusters (groups) of neurons emerge in the network due to the visual stimulation. The different clusters correspond to different visual stimuli. There is no limitation of the number of stimuli. Distinct clusters do not perturb each other, although the coupling strength between all model neurons is of the same order of magnitude. Our analysis is not restricted to weak coupling strength. The scatter of the couplings causes shifts of the cluster frequencies. The model's behavior is compared with the experimental findings. The coupling mechanism is extended in order to model the influence of bicucullin upon the neural network. We additionally investigate repulsive couplings, which lead to constant phase differences between clusters of the same frequency. Finally, we consider the problem of selective attention from the viewpoint of our model.     

The generation of phase differences and frequency changes in a network model of inferior olive subthreshold oscillations

It is commonly accepted that the Inferior Olive (IO) provides a timing signal to the cerebellum. Stable subthreshold oscillations in the IO can facilitate accurate timing by phase-locking spikes to the peaks of the oscillation. Several theoretical models accounting for the synchronized subthreshold oscillations have been proposed, however, two experimental observations remain an enigma. The first is the observation of frequent alterations in the frequency of the oscillations. The second is the observation of constant phase differences between simultaneously recorded neurons. In order to account for these two observations we constructed a canonical network model based on anatomical and physiological data from the IO. The constructed network is characterized by clustering of neurons with similar conductance densities, and by electrical coupling between neurons. Neurons inside a cluster are densely connected with weak strengths, while neurons belonging to different clusters are sparsely connected with stronger connections. We found that this type of network can robustly display stable subthreshold oscillations. The overall frequency of the network changes with the strength of the inter-cluster connections, and phase differences occur between neurons of different clusters. Moreover, the phase differences provide a mechanistic explanation for the experimentally observed propagating waves of activity in the IO. We conclude that the architecture of the network of electrically coupled neurons in combination with modulation of the inter-cluster coupling strengths can account for the experimentally observed frequency changes and the phase differences.     

Models of central pattern generators for quadruped locomotion I. Primary gaits

In this paper we continue the analysis of a network of symmetrically coupled cells modeling central pattern generators for quadruped locomotion proposed by Golubitsky, Stewart, Buono, and Collins. By a cell we mean a system of ordinary differential equations and by a coupled cell system we mean a network of identical cells with coupling terms. We have three main results in this paper. First, we show that the proposed network is the simplest one modeling the common quadruped gaits of walk, trot, and pace. In doing so we prove a general theorem classifying spatio-temporal symmetries of periodic solutions to equivariant systems of differential equations. We also specialize this theorem to coupled cell systems. Second, this paper focuses on primary gaits; that is, gaits that are modeled by output signals from the central pattern generator where each cell emits the same waveform along with exact phase shifts between cells. Our previous work showed that the network is capable of producing six primary gaits. Here, we show that under mild assumptions on the cells and the coupling of the network, primary gaits can be produced from Hopf bifurcation by varying only coupling strengths of the network. Third, we discuss the stability of primary gaits and exhibit these solutions by performing numerical simulations using the dimensionless Morris-Lecar equations for the cell dynamics.     

A linear model for characterization of synchronization frequencies of neural networks

The synchronization frequency of neural networks and its dynamics have important roles in deciphering the working mechanisms of the brain. It has been widely recognized that the properties of functional network synchronization and its dynamics are jointly determined by network topology, network connection strength, i.e., the connection strength of different edges in the network, and external input signals, among other factors. However, mathematical and computational characterization of the relationships between network synchronization frequency and these three important factors are still lacking. This paper presents a novel computational simulation framework to quantitatively characterize the relationships between neural network synchronization frequency and network attributes and input signals. Specifically, we constructed a series of neural networks including simulated small-world networks, real functional working memory network derived from functional magnetic resonance imaging, and real large-scale structural brain networks derived from diffusion tensor imaging, and performed synchronization simulations on these networks via the Izhikevich neuron spiking model. Our experiments demonstrate that both of the network synchronization strength and synchronization frequency change according to the combination of input signal frequency and network self-synchronization frequency. In particular, our extensive experiments show that the network synchronization frequency can be represented via a linear combination of the network self-synchronization frequency and the input signal frequency. This finding could be attributed to an intrinsically-preserved principle in different types of neural systems, offering novel insights into the working mechanism of neural systems.     

Synaptic information transfer in computer models of neocortical columns

Understanding the direction and quantity of information flowing in neuronal networks is a fundamental problem in neuroscience. Brains and neuronal networks must at the same time store information about the world and react to information in the world. We sought to measure how the activity of the network alters information flow from inputs to output patterns. Using neocortical column neuronal network simulations, we demonstrated that networks with greater internal connectivity reduced input/output correlations from excitatory synapses and decreased negative correlations from inhibitory synapses, measured by Kendall’s τ correlation. Both of these changes were associated with reduction in information flow, measured by normalized transfer entropy (nTE). Information handling by the network reflected the degree of internal connectivity. With no internal connectivity, the feedforward network transformed inputs through nonlinear summation and thresholding. With greater connectivity strength, the recurrent network translated activity and information due to contribution of activity from intrinsic network dynamics. This dynamic contribution amounts to added information drawn from that stored in the network. At still higher internal synaptic strength, the network corrupted the external information, producing a state where little external information came through. The association of increased information retrieved from the network with increased gamma power supports the notion of gamma oscillations playing a role in information processing.     

Stability switches,oscillatory multistability,and spatio-temporal patterns of nonlinear oscillations in recurrently delay coupled neural networks

A model of time-delay recurrently coupled spatially segregated neural assemblies is here proposed. We show that it operates like some of the hierarchical architectures of the brain. Each assembly is a neural network with no delay in the local couplings between the units. The delay appears in the long range feedforward and feedback inter-assemblies communications. Bifurcation analysis of a simple four-units system in the autonomous case shows the richness of the dynamical behaviors in a biophysically plausible parameter region. We find oscillatory multistability, hysteresis, and stability switches of the rest state provoked by the time delay. Then we investigate the spatio-temporal patterns of bifurcating periodic solutions by using the symmetric local Hopf bifurcation theory of delay differential equations and derive the equation describing the flow on the center manifold that enables us determining the direction of Hopf bifurcations and stability of the bifurcating periodic orbits. We also discuss computational properties of the system due to the delay when an external drive of the network mimicks external sensory input.     

Rigid‐body motions of interacting proteins dominate multispecific binding of ubiquitin in a shape‐dependent manner

To understand the dynamic aspects of multispecificity of ubiquitin, we studied nine ubiquitin–ligand (partner protein) complexes by normal mode analysis based on an elastic network model. The coupling between ubiquitin and ligand motions was analyzed by decomposing it into rigid‐body (external) and vibrational (internal) motions of each subunit. We observed that in total the external motions in one of the subunits largely dominated the coupling. The combination of external motions of ubiquitin and the ligands showed general trends of rotations and translations. Moreover, we observed that the rotational motions of ubiquitin were correlated to the ligand orientations. We also identified ubiquitin atomic vibrations that differentiated the orientation of the ligand molecule. We observed that the extents of coupling were correlated to the shapes of the ligands, and this trend was more pronounced when the coupling involved vibrational motions of the ligand. In conclusion, an intricate interplay between internal and external motions of ubiquitin and the ligands help understand the dynamics of multispecificity, which is mostly guided by the shapes of the ligands and the complex. Proteins 2014; 82:77–89. © 2013 Wiley Periodicals, Inc.     

Synchrony and Asynchrony for Neuronal Dynamics Defined on Complex Networks

We describe and analyze a model for a stochastic pulse-coupled neuronal network with many sources of randomness: random external input, potential synaptic failure, and random connectivity topologies. We show that different classes of network topologies give rise to qualitatively different types of synchrony: uniform (Erdős–Rényi) and “small-world” networks give rise to synchronization phenomena similar to that in “all-to-all” networks (in which there is a sharp onset of synchrony as coupling is increased); in contrast, in “scale-free” networks the dependence of synchrony on coupling strength is smoother. Moreover, we show that in the uniform and small-world cases, the fine details of the network are not important in determining the synchronization properties; this depends only on the mean connectivity. In contrast, for scale-free networks, the dynamics are significantly affected by the fine details of the network; in particular, they are significantly affected by the local neighborhoods of the “hubs” in the network.     

Boolean network-based analysis of the apoptosis network: Irreversible apoptosis and stable surviving

To understand the design principles of the molecular interaction network associated with the irreversibility of cell apoptosis and the stability of cell surviving, we constructed a Boolean network integrating both the intrinsic and extrinsic pro-apoptotic pathways with pro-survival signal transduction pathways. We performed statistical analyses of the dependences of cell fate on initial states and on input signals. The analyses reproduced the well-known pro- and anti-apoptotic effects of key external signals and network components. We found that the external GF signal by itself did not change the apoptotic ratio from randomly chosen initial states when there is no external TNF signal, but can significantly offset apoptosis induced by the TNF signal. While a complete model produces the expected irreversibility of the apoptosis process, alternative models missing one or more of four selected inter-component connections indicate that the feedback loops directly involving the caspase 3 are essential for maintaining irreversibility of apoptosis. The feedback loops involving P53 showed compensating effects when those involving caspase 3 have been removed. The GF signal significantly increases the stability of the surviving states of the network. The apoptosis network seems to use different modules by design to control the irreversibility of the apoptosis process and the stability of the surviving states. Such a design may accommodate the needed plasticity for the network to adapt to different cellular environments: depending on the strength of external pro-surviving signals, apoptosis can be induced either easily or difficultly by pro-apoptotic signal of varying strengths, but proceed with invariable irreversibility.     

Motor pattern selection by combinatorial code of interneuronal pathways

We use a modeling approach to examine ideas derived from physiological network analyses, pertaining to the switch of a motor control network between two opposite control modes. We studied the femur–tibia joint control system of the insect leg, and its switch between resistance reflex in posture control and “active reaction” in walking, both elicited by the same sensory input. The femur–tibia network was modeled by fitting the responses of model neurons to those obtained in animals. The strengths of 16 interneuronal pathways that integrate sensory input were then assigned three different values and varied independently, generating a database of more than 43 million network variants. We demonstrate that the same neural network can produce the two different behaviors, depending on the combinatorial code of interneuronal pathways. That is, a switch between behaviors, such as standing to walking, can be brought about by altering the strengths of selected sensory integration pathways. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Degree of Contracture Related to Residual Muscle Shoulder Strength in Children with Obstetric Brachial Plexus Lesions

Background and Objectives

 Little is known about the relation between residual muscle strength and joint contracture formation in neuromuscular disorders. This study aimed to investigate the relation between residual muscle strength and shoulder joint contractures in children with sequelae of obstetric brachial plexus lesion (OBPL). In OBPL a shoulder joint contracture is a frequent finding. We hypothesize that residual internal and external rotator strength and their balance are related to the extent of shoulder joint contracture.

Methods

 Clinical assessment was performed in 34 children (mean 10.0 years) with unilateral OBPL and Narakas classes I–III. External and internal rotation strengths were measured with the shoulder in neutral position using a handheld dynamometer. Strength on the affected side was given as percentage of the normal side. Contracture was assessed by passive internal and external rotations in degrees (in 0° abduction). Mallet classification was used for active shoulder function.

Results

 External and internal rotation strengths on the affected side were approximately 50% of the normal side and on average both equally affected: 56% (SD 18%) respectively 51% (SD 27%); r = 0.600, p = 0.000. Residual strengths were not related to passive internal or external rotation (p > 0.200). Internal rotation strength (r =  − 0.425, p <0.05) was related to Narakas class. Mallet score was related to external and internal rotation strengths (r = 0.451 and r = 0.515, respectively; p < 0.01).

Conclusion

 The intuitive notion that imbalances in residual muscle strength influence contracture formation cannot be confirmed in this study. Our results are of interest for the understanding of contracture formation in OBPL.     

Stability and structural constraints of random brain networks with excitatory and inhibitory neural populations

The stability of brain networks with randomly connected excitatory and inhibitory neural populations is investigated using a simplified physiological model of brain electrical activity. Neural populations are randomly assigned to be excitatory or inhibitory and the stability of a brain network is determined by the spectrum of the network’s matrix of connection strengths. The probability that a network is stable is determined from its spectral density which is numerically determined and is approximated by a spectral distribution recently derived by Rajan and Abbott. The probability that a brain network is stable is maximum when the total connection strength into a population is approximately zero and is shown to depend on the arrangement of the excitatory and inhibitory connections and the parameters of the network. The maximum excitatory and inhibitory input into a structure allowed by stability occurs when the net input equals zero and, in contrast to networks with randomly distributed excitatory and inhibitory connections, substantially increases as the number of connections increases. Networks with the largest excitatory and inhibitory input allowed by stability have multiple marginally stable modes, are highly responsive and adaptable to external stimuli, have the same total input into each structure with minimal variance in the excitatory and inhibitory connection strengths, and have a wide range of flexible, adaptable, and complex behavior.     

脑电事件相关去同步化和同步化的神经元群模型

利用基于丘脑-皮层网络的神经元群模型,研究被试者在某种认知状态下脑功能区的连接状态。模型包括三个模块,分别对应脑电头皮电极C3、Cz、C4记录的三个皮质区。模型外部输入包括用高斯白噪声表示的上行传入感受器信号、用直流偏移表示的皮质对丘脑的兴奋性输入、用指数衰减表示的来自脑千和前脑基底神经元的调制信号。模型输出的兴奋性神经元群的平均膜电位反映脑电记录的局部电位。改变模型输入,进行多次仿真试验并进行线性和非线性分析。研究结果显示：仿真输出信号的alpha频带功率谱有与实际脑机接口实验一致的事件相关去同步化和同步化现象;模型中功能相近的区域间有更强的耦合,随着耦合强度的增加,输出信号间的相关性和同步性均增加。