Does abnormal spinal reciprocal inhibition lead to co-contraction of antagonist motor units? A modeling study

It is suggested that co-contraction of antagonist motor units perhaps due to abnormal disynaptic I(a) reciprocal inhibition is responsible for Parkinsonian rigidity. A neural model of Parkinson's disease bradykinesia is extended to incorporate the effects of spindle feedback on key cortical cells and examine the effects of dopamine depletion on spinal activities. Simulation results show that although reciprocal inhibition is reduced in DA depleted case, it doesn't lead to co-contraction of antagonist motor neurons. Implications to Parkinsonian rigidity are discussed.     

Cortical Resonance Frequencies Emerge from Network Size and Connectivity

Neural oscillations occur within a wide frequency range with different brain regions exhibiting resonance-like characteristics at specific points in the spectrum. At the microscopic scale, single neurons possess intrinsic oscillatory properties, such that is not yet known whether cortical resonance is consequential to neural oscillations or an emergent property of the networks that interconnect them. Using a network model of loosely-coupled Wilson-Cowan oscillators to simulate a patch of cortical sheet, we demonstrate that the size of the activated network is inversely related to its resonance frequency. Further analysis of the parameter space indicated that the number of excitatory and inhibitory connections, as well as the average transmission delay between units, determined the resonance frequency. The model predicted that if an activated network within the visual cortex increased in size, the resonance frequency of the network would decrease. We tested this prediction experimentally using the steady-state visual evoked potential where we stimulated the visual cortex with different size stimuli at a range of driving frequencies. We demonstrate that the frequency corresponding to peak steady-state response inversely correlated with the size of the network. We conclude that although individual neurons possess resonance properties, oscillatory activity at the macroscopic level is strongly influenced by network interactions, and that the steady-state response can be used to investigate functional networks.     

Activity of precentral neurons during torque-triggered hand movements in awake primates.

In monkeys performing a handle-repositioning task involving primarily wrist flexion-extension (F-E) movements after a torque perturbation was delivered to the handle, single units were recorded extracellularly in the contralateral precentral cortex. Precentral neurons were identified by passive somatosensory stimulation, and were classified into five somatotopically organized populations. Based on electromyographic recordings, it was observed that flexors and extensors about the wrist joint were specifically involved in the repositioning of the handle, while many other muscles which act at the wrist and other forelimb joints were involved in the task in a supportive role. In precentral cortex, all neuronal responses observed were temporally correlated to both the sensory stimuli and the motor responses. Visual stimuli, presented simultaneously with torque perturbations, did not affect the early portion of cortical responses to such torque perturbations. In each of the five somatotopically organized neuronal populations, task-related neurons as well as task-unrelated ones were observed. A significantly larger proportion of wrist (F-E) neurons was related to the task, as compared with the other, nonwrist (F-E) populations. The above findings were discussed in the context of a hypothesis for the function of precentral cortex during voluntary limb movement in awake primates. This hypothesis incorporates a relationship between activities of populations of precentral neurons, defined with respect to their responses to peripheral events at or about single joints, and movements about the same joint.     

The Importance of Lateral Connections in the Parietal Cortex for Generating Motor Plans

Substantial evidence has highlighted the significant role of associative brain areas, such as the posterior parietal cortex (PPC) in transforming multimodal sensory information into motor plans. However, little is known about how different sensory information, which can have different delays or be absent, combines to produce a motor plan, such as executing a reaching movement. To address these issues, we constructed four biologically plausible network architectures to simulate PPC: 1) feedforward from sensory input to the PPC to a motor output area, 2) feedforward with the addition of an efference copy from the motor area, 3) feedforward with the addition of lateral or recurrent connectivity across PPC neurons, and 4) feedforward plus efference copy, and lateral connections. Using an evolutionary strategy, the connectivity of these network architectures was evolved to execute visually guided movements, where the target stimulus provided visual input for the entirety of each trial. The models were then tested on a memory guided motor task, where the visual target disappeared after a short duration. Sensory input to the neural networks had sensory delays consistent with results from monkey studies. We found that lateral connections within the PPC resulted in smoother movements and were necessary for accurate movements in the absence of visual input. The addition of lateral connections resulted in velocity profiles consistent with those observed in human and non-human primate visually guided studies of reaching, and allowed for smooth, rapid, and accurate movements under all conditions. In contrast, Feedforward or Feedback architectures were insufficient to overcome these challenges. Our results suggest that intrinsic lateral connections are critical for executing accurate, smooth motor plans.     

The functional role of different neural activation profiles during precision grip: An artificial neural network approach

A dynamic and recurrent artificial neural network was used to investigate the functional properties of firing patterns observed in the primary motor (M1) and the primary somatosensory (S1) cortex of the behaving monkey during control of precision grip force. In the behaving monkey it was found that neurons in M1 and in S1 increase their firing activity with increasing grip force, as do the intrinsic and extrinsic hand muscles implicated in the task. However, some neurons also decreased their activity as a function of increasing force. The functional implication of these latter neurons is not clear and has not been elucidated so far. In order to explore their functional implication, we therefore simulated patterns of neural activity in artificial neural networks that represent cortical, spinal and afferent neural populations and tested whether particular activity profiles would emerge as a function of the input and of the connectivity of these networks. The functional implication of units with emergent or imposed decreasing activity was then explored.Decreasing patterns of activity in M1 units did not emerge from the networks. However, the same networks generated decreasing activity if imposed as target patterns. As indicated by the emerging weight space, M1 projection units with decreasing patterns are functionally less involved in driving alpha motoneurons than units with increasing profiles. Furthermore, these units did not provide significant fusimotor drive, whereas those with increasing profiles did. Fusimotor drive was a function of the (imposed) form of muscle spindle afferent activity: with gamma (fusimotor) drive, muscle spindle afferents provided signals other than muscle length (as observed experimentally). The network solutions thus predict a functional dichotomy between increasing and decreasing M1 neurons: the former primarily drive alpha and gamma motoneurons, the latter only weakly alpha motoneurons.     

Parkinsonian tremor and simplification in network dynamics

We explore the behavior of richly connected inhibitory neural networks under parameter changes that correspond to weakening of synaptic efficacies between network units, and show that transitions from irregular to periodic dynamics are common in such systems. The weakening of these connections leads to a reduction in the number of units that effectively drive the dynamics and thus to simpler behavior. We hypothesize that the multiple interconnecting loops of the brain’s motor circuitry, which involve many inhibitory connections, exhibit such transitions. Normal physiological tremor is irregular while other forms of tremor show more regular oscillations. Tremor in Parkinson’s disease, for example, stems from weakened synaptic efficacies of dopaminergic neurons in the nigro-striatal pathway, as in our general model. The multiplicity of structures involved in the production of symptoms in Parkinson’s disease and the reversibility of symptoms by pharmacological and surgical manipulation of connection parameters suggest that such a neural network model is appropriate. Furthermore, fixed points that can occur in the network models are suggestive of akinesia in Parkinson’s disease. This model is consistent with the view that normal physiological systems can be regulated by robust and richly connected feedback networks with complex dynamics, and that loss of complexity in the feedback structure due to disease leads to more orderly behavior.     

Recurrent neural networks of integrate-and-fire cells simulating short-term memory and wrist movement tasks derived from continuous dynamic networks

Dynamic recurrent neural networks composed of units with continuous activation functions provide a powerful tool for simulating a wide range of behaviors, since the requisite interconnections can be readily derived by gradient descent methods. However, it is not clear whether more realistic integrate-and-fire cells with comparable connection weights would perform the same functions. We therefore investigated methods to convert dynamic recurrent neural networks of continuous units into networks with integrate-and-fire cells. The transforms were tested on two recurrent networks derived by backpropagation. The first simulates a short-term memory task with units that mimic neural activity observed in cortex of monkeys performing instructed delay tasks. The network utilizes recurrent connections to generate sustained activity that codes the remembered value of a transient cue. The second network simulates patterns of neural activity observed in monkeys performing a step-tracking task with flexion/extension wrist movements. This more complicated network provides a working model of the interactions between multiple spinal and supraspinal centers controlling motoneurons. Our conversion algorithm replaced each continuous unit with multiple integrate-and-fire cells that interact through delayed "synaptic potentials". Successful transformation depends on obtaining an appropriate fit between the activation function of the continuous units and the input-output relation of the spiking cells. This fit can be achieved by adapting the parameters of the synaptic potentials to replicate the input-output behavior of a standard sigmoidal activation function (shown for the short-term memory network). Alternatively, a customized activation function can be derived from the input-output relation of the spiking cells for a chosen set of parameters (demonstrated for the wrist flexion/extension network). In both cases the resulting networks of spiking cells exhibited activity that replicated the activity of corresponding continuous units. This confirms that the network solutions obtained through backpropagation apply to spiking networks and provides a useful method for deriving recurrent spiking networks performing a wide range of functions.     

Red nucleus: role in motor control

Experimental reports in the past year have provided a better understanding of the motor functions of excitatory and inhibitory neurotransmitters in the red nucleus, and of the sensorimotor properties of single rubral neurons. These data fit well within the framework of a neural network model of the rubrocerebellar system.     

Physiological and pathological oscillatory networks in the human motor system.

Human brain functions are heavily contingent on neural interactions both at the single neuron and the neural population or system level. Accumulating evidence from neurophysiological studies strongly suggests that coupling of oscillatory neural activity provides an important mechanism to establish neural interactions. With the availability of whole-head magnetoencephalography (MEG) macroscopic oscillatory activity can be measured non-invasively from the human brain with high temporal and spatial resolution. To localise, quantify and map oscillatory activity and interactions onto individual brain anatomy we have developed the 'dynamic imaging of coherent sources' (DICS) method which allows to identify and analyse cerebral oscillatory networks from MEG recordings. Using this approach we have characterized physiological and pathological oscillatory networks in the human sensorimotor system. Coherent 8 Hz oscillations emerge from a cerebello-thalamo-premotor-motor cortical network and exert an 8 Hz oscillatory drive on the spinal motor neurons which can be observed as a physiological tremulousness of the movement termed movement discontinuities. This network represents the neurophysiological substrate of a discrete mode of motor control. In parkinsonian resting tremor we have identified an extensive cerebral network consisting of primary motor and lateral premotor cortex, supplementary motor cortex, thalamus/basal ganglia, posterior parietal cortex and secondary somatosensory cortex, which are entrained in the tremor or twice the tremor rhythm. This low frequency entrapment of motor areas likely plays an important role in the pathophysiology of parkinsonian motor symptoms. Finally, studies on patients with postural tremor in hepatic encephalopathy revealed that this type of tremor results from a pathologically slow thalamocortical and cortico-muscular coupling during isometric hold tasks. In conclusion, the analysis of oscillatory cerebral networks provides new insights into physiological mechanisms of motor control and pathophysiological mechanisms of tremor disorders.     

Neural Coding of Finger and Wrist Movements

Previous work (Schieber and Hibbard, 1993) has shown that single motor cortical neurons do not discharge specifically for a particular flexion-extension finger movement but instead are active with movements of different fingers. In addition, neuronal populations active with movements of different fingers overlap extensively in their spatial locations in the motor cortex. These data suggested that control of any finger movement utilizes a distributed population of neurons. In this study we applied the neuronal population vector analysis (Georgopoulos et al., 1983) to these same data to determine (1) whether single cells are tuned in an abstract, three-dimensional (3D) instructed finger and wrist movement space with hand-like geometry and (2) whether the neuronal population encodes specific finger movements. We found that the activity of 132/176 (75%) motor cortical neurons related to finger movements was indeed tuned in this space. Moreover, the population vector computed in this space predicted well the instructed finger movement. Thus, although single neurons may be related to several disparate finger movements, and neurons related to different finger movements are intermingled throughout the hand area of the motor cortex, the neuronal population activity does specify particular finger movements.     

Driving reservoir models with oscillations: a solution to the extreme structural sensitivity of chaotic networks

A large body of experimental and theoretical work on neural coding suggests that the information stored in brain circuits is represented by time-varying patterns of neural activity. Reservoir computing, where the activity of a recurrently connected pool of neurons is read by one or more units that provide an output response, successfully exploits this type of neural activity. However, the question of system robustness to small structural perturbations, such as failing neurons and synapses, has been largely overlooked. This contrasts with well-studied dynamical perturbations that lead to divergent network activity in the presence of chaos, as is the case for many reservoir networks. Here, we distinguish between two types of structural network perturbations, namely local (e.g., individual synaptic or neuronal failure) and global (e.g., network-wide fluctuations). Surprisingly, we show that while global perturbations have a limited impact on the ability of reservoir models to perform various tasks, local perturbations can produce drastic effects. To address this limitation, we introduce a new architecture where the reservoir is driven by a layer of oscillators that generate stable and repeatable trajectories. This model outperforms previous implementations while being resistant to relatively large local and global perturbations. This finding has implications for the design of reservoir models that capture the capacity of brain circuits to perform cognitively and behaviorally relevant tasks while remaining robust to various forms of perturbations. Further, our work proposes a novel role for neuronal oscillations found in cortical circuits, where they may serve as a collection of inputs from which a network can robustly generate complex dynamics and implement rich computations.     

A dynamical neural network model for motor cortical activity during movement: population coding of movement trajectories

As a dynamical model for motor cortical activity during hand movement we consider an artificial neural network that consists of extensively interconnected neuron-like units and performs the neuronal population vector operations. Local geometrical parameters of a desired curve are introduced into the network as an external input. The output of the model is a time-dependent direction and length of the neuronal population vector which is calculated as a sum of the activity of directionally tuned neurons in the ensemble. The main feature of the model is that dynamical behavior of the neuronal population vector is the result of connections between directionally tuned neurons rather than being imposed externally. The dynamics is governed by a system of coupled nonlinear differential equations. Connections between neurons are assigned in the simplest and most common way so as to fulfill basic requirements stemming from experimental findings concerning the directional tuning of individual neurons and the stabilization of the neuronal population vector, as well as from previous theoretical studies. The dynamical behavior of the model reveals a close similarity with the experimentally observed dynamics of the neuronal population vector. Specifically, in the framework of the model it is possible to describe a geometrical curve in terms of the time series of the population vector. A correlation between the dynamical behavior of the direction and the length of the population vector entails a dependence of the neural velocity on the curvature of the tracing trajectory that corresponds well to the experimentally measured covariation between tangential velocity and curvature in drawing tasks.On leave of absencefrom the Institute of Molecular Genetics, Russian Academy of Sciences, Moscow, Russia.     

A zonal model of cortical functions

A model of cortical functions is developed with the object of simulating the observed behavior of individual neurons organized in unit circuits and functional systems of the cerebellum, the cerebrum and the hippocampal formation. The neuronal model is capable of representing refractory and potentiated states, as well as the firing and lowest resting states. The unit circuits of each system consist of all common types of cells with known synaptic connections. In the cerebral system these unit circuits are interconnected to form columns as well as zones. A new discrete neural network equation, which takes account of interactions with the extracellular field, is proposed to simulate electrical activity in these circuits. A coherent theory of cortical activity and functions is derived that accounts for many of the observed phenomena, including those associated with the development of long-term potentiation and sequential memory. Three appendices are devoted to the theory of extracellular interactions, the derivation of non-linear network equations, and a computer program to simulate learning in the cortex.     

On how network architecture determines the dominant patterns of spontaneous neural activity

In the absence of sensory stimulation, neocortical circuits display complex patterns of neural activity. These patterns are thought to reflect relevant properties of the network, including anatomical features like its modularity. It is also assumed that the synaptic connections of the network constrain the repertoire of emergent, spontaneous patterns. Although the link between network architecture and network activity has been extensively investigated in the last few years from different perspectives, our understanding of the relationship between the network connectivity and the structure of its spontaneous activity is still incomplete. Using a general mathematical model of neural dynamics we have studied the link between spontaneous activity and the underlying network architecture. In particular, here we show mathematically how the synaptic connections between neurons determine the repertoire of spatial patterns displayed in the spontaneous activity. To test our theoretical result, we have also used the model to simulate spontaneous activity of a neural network, whose architecture is inspired by the patchy organization of horizontal connections between cortical columns in the neocortex of primates and other mammals. The dominant spatial patterns of the spontaneous activity, calculated as its principal components, coincide remarkably well with those patterns predicted from the network connectivity using our theory. The equivalence between the concept of dominant pattern and the concept of attractor of the network dynamics is also demonstrated. This in turn suggests new ways of investigating encoding and storage capabilities of neural networks.     

Motor patterns during kicking movements in the locust

Locusts (Schistocerca gregaria) use a distinctive motor pattern to extend the tibia of a hind leg rapidly in a kick. The necessary force is generated by an almost isometric contraction of the extensor tibiae muscle restrained by the co-contraction of the flexor tibiae (co-contraction phase) and aided by the mechanics of the femoro-tibial joint. The stored energy is delivered suddenly when the flexor muscle is inhibited. This paper analyses the activity of motor neurons to the major hind leg muscles during kicking, and relates it to tibial movements and the resultant forces.During the co-contraction phase flexor tibiae motor neurons are driven by apparently common sources of synaptic inputs to depolarized plateaus at which they spike. The two excitatory extensor motor neurons are also depolarized by similar patterns of synaptic inputs, but with the slow producing more spikes at higher frequencies than the fast. Trochanteral depressors spike at high frequency, the single levator tarsi at low frequency, and common inhibitors 2 and 3 spike sporadically. Trochanteral levators, depressor tarsi, and a retractor unguis motor neuron are hyperpolarized.Before the tibia extends all flexor motor neurons are hyperpolarized simultaneously, two common inhibitors, and the levator trochanter and depressor tarsi motor neurons are depolarized. Later, but still before the tibial movement starts, the extensor tibiae and levator tarsi motor neurons are hyperpolarized. After the movement has started, the extensor motor neurons are hyperpolarized further and the depressor trochanteris motor neurons are also hyperpolarized, indicating a contribution of both central and sensory feedback pathways.Variations in the duration of the co-contraction of almost twenty-fold, and in the number of spikes in the fast extensor tibiae motor neuron from 2–50 produce a spectrum of tibial extensions ranging from slow and weak, to rapid and powerful. Flexibility in the networks producing the motor pattern therefore results in a range of movements suited to the fluctuating requirements of the animal.     

A thalamo-cortical neural mass model for the simulation of brain rhythms during sleep

Cortico-thalamic interactions are known to play a pivotal role in many brain phenomena, including sleep, attention, memory consolidation and rhythm generation. Hence, simple mathematical models that can simulate the dialogue between the cortex and the thalamus, at a mesoscopic level, have a great cognitive value. In the present work we describe a neural mass model of a cortico-thalamic module, based on neurophysiological mechanisms. The model includes two thalamic populations (a thalamo-cortical relay cell population, TCR, and its related thalamic reticular nucleus, TRN), and a cortical column consisting of four connected populations (pyramidal neurons, excitatory interneurons, inhibitory interneurons with slow and fast kinetics). Moreover, thalamic neurons exhibit two firing modes: bursting and tonic. Finally, cortical synapses among pyramidal neurons incorporate a disfacilitation mechanism following prolonged activity. Simulations show that the model is able to mimic the different patterns of rhythmic activity in cortical and thalamic neurons (beta and alpha waves, spindles, delta waves, K-complexes, slow sleep waves) and their progressive changes from wakefulness to deep sleep, by just acting on modulatory inputs. Moreover, simulations performed by providing short sensory inputs to the TCR show that brain rhythms during sleep preserve the cortex from external perturbations, still allowing a high cortical activity necessary to drive synaptic plasticity and memory consolidation. In perspective, the present model may be used within larger cortico-thalamic networks, to gain a deeper understanding of mechanisms beneath synaptic changes during sleep, to investigate the specific role of brain rhythms, and to explore cortical synchronization achieved via thalamic influences.     

Modeling the Impact of Lesions in the Human Brain

Lesions of anatomical brain networks result in functional disturbances of brain systems and behavior which depend sensitively, often unpredictably, on the lesion site. The availability of whole-brain maps of structural connections within the human cerebrum and our increased understanding of the physiology and large-scale dynamics of cortical networks allow us to investigate the functional consequences of focal brain lesions in a computational model. We simulate the dynamic effects of lesions placed in different regions of the cerebral cortex by recording changes in the pattern of endogenous (“resting-state”) neural activity. We find that lesions produce specific patterns of altered functional connectivity among distant regions of cortex, often affecting both cortical hemispheres. The magnitude of these dynamic effects depends on the lesion location and is partly predicted by structural network properties of the lesion site. In the model, lesions along the cortical midline and in the vicinity of the temporo-parietal junction result in large and widely distributed changes in functional connectivity, while lesions of primary sensory or motor regions remain more localized. The model suggests that dynamic lesion effects can be predicted on the basis of specific network measures of structural brain networks and that these effects may be related to known behavioral and cognitive consequences of brain lesions.     

Cellularly-driven differences in network synchronization propensity are differentially modulated by firing frequency

Spatiotemporal pattern formation in neuronal networks depends on the interplay between cellular and network synchronization properties. The neuronal phase response curve (PRC) is an experimentally obtainable measure that characterizes the cellular response to small perturbations, and can serve as an indicator of cellular propensity for synchronization. Two broad classes of PRCs have been identified for neurons: Type I, in which small excitatory perturbations induce only advances in firing, and Type II, in which small excitatory perturbations can induce both advances and delays in firing. Interestingly, neuronal PRCs are usually attenuated with increased spiking frequency, and Type II PRCs typically exhibit a greater attenuation of the phase delay region than of the phase advance region. We found that this phenomenon arises from an interplay between the time constants of active ionic currents and the interspike interval. As a result, excitatory networks consisting of neurons with Type I PRCs responded very differently to frequency modulation compared to excitatory networks composed of neurons with Type II PRCs. Specifically, increased frequency induced a sharp decrease in synchrony of networks of Type II neurons, while frequency increases only minimally affected synchrony in networks of Type I neurons. These results are demonstrated in networks in which both types of neurons were modeled generically with the Morris-Lecar model, as well as in networks consisting of Hodgkin-Huxley-based model cortical pyramidal cells in which simulated effects of acetylcholine changed PRC type. These results are robust to different network structures, synaptic strengths and modes of driving neuronal activity, and they indicate that Type I and Type II excitatory networks may display two distinct modes of processing information.     

Associative neural network model for the generation of temporal patterns. Theory and application to central pattern generators.

Cyclic patterns of motor neuron activity are involved in the production of many rhythmic movements, such as walking, swimming, and scratching. These movements are controlled by neural circuits referred to as central pattern generators (CPGs). Some of these circuits function in the absence of both internal pacemakers and external feedback. We describe an associative neural network model whose dynamic behavior is similar to that of CPGs. The theory predicts the strength of all possible connections between pairs of neurons on the basis of the outputs of the CPG. It also allows the mean operating levels of the neurons to be deduced from the measured synaptic strengths between the pairs of neurons. We apply our theory to the CPG controlling escape swimming in the mollusk Tritonia diomedea. The basic rhythmic behavior is shown to be consistent with a simplified model that approximates neurons as threshold units and slow synaptic responses as elementary time delays. The model we describe may have relevance to other fixed action behaviors, as well as to the learning, recall, and recognition of temporally ordered information.