Jeanette Hellgren Kotaleski

Modelling the molecular mechanisms of synaptic plasticity using systems biology approaches

Synaptic plasticity is thought to underlie learning and memory, but the complexity of the interactions between the ion channels, enzymes and genes that are involved in synaptic plasticity impedes a deep understanding of this phenomenon. Computer modelling has been used to investigate the information processing that is performed by the signalling pathways involved in synaptic plasticity in principal neurons of the hippocampus, striatum and cerebellum. In the past few years, new software developments that combine computational neuroscience techniques with systems biology techniques have allowed large-scale, kinetic models of the molecular mechanisms underlying long-term potentiation and long-term depression. We highlight important advancements produced by these quantitative modelling efforts and introduce promising approaches that use advancements in live-cell imaging.

Transient Calcium and Dopamine Increase PKA Activity and DARPP-32 Phosphorylation

Reinforcement learning theorizes that strengthening of synaptic connections in medium spiny neurons of the striatum occurs when glutamatergic input (from cortex) and dopaminergic input (from substantia nigra) are received simultaneously. Subsequent to learning, medium spiny neurons with strengthened synapses are more likely to fire in response to cortical input alone. This synaptic plasticity is produced by phosphorylation of AMPA receptors, caused by phosphorylation of various signalling molecules. A key signalling molecule is the phosphoprotein DARPP-32, highly expressed in striatal medium spiny neurons. DARPP-32 is regulated by several neurotransmitters through a complex network of intracellular signalling pathways involving cAMP (increased through dopamine stimulation) and calcium (increased through glutamate stimulation). Since DARPP-32 controls several kinases and phosphatases involved in striatal synaptic plasticity, understanding the interactions between cAMP and calcium, in particular the effect of transient stimuli on DARPP-32 phosphorylation, has major implications for understanding reinforcement learning. We developed a computer model of the biochemical reaction pathways involved in the phosphorylation of DARPP-32 on Thr34 and Thr75. Ordinary differential equations describing the biochemical reactions were implemented in a single compartment model using the software XPPAUT. Reaction rate constants were obtained from the biochemical literature. The first set of simulations using sustained elevations of dopamine and calcium produced phosphorylation levels of DARPP-32 similar to that measured experimentally, thereby validating the model. The second set of simulations, using the validated model, showed that transient dopamine elevations increased the phosphorylation of Thr34 as expected, but transient calcium elevations also increased the phosphorylation of Thr34, contrary to what is believed. When transient calcium and dopamine stimuli were paired, PKA activation and Thr34 phosphorylation increased compared with dopamine alone. This result, which is robust to variation in model parameters, supports reinforcement learning theories in which activity-dependent long-term synaptic plasticity requires paired glutamate and dopamine inputs.

GABAergic Circuits Control Spike-Timing-Dependent Plasticity

The spike-timing-dependent plasticity (STDP), a synaptic learning rule for encoding learning and memory, relies on relative timing of neuronal activity on either side of the synapse. GABAergic signaling has been shown to control neuronal excitability and consequently the spike timing, but whether GABAergic circuits rule the STDP remained unknown. Here we show that GABAergic signaling governs the polarity of STDP, because blockade of GABAA receptors was able to completely reverse the temporal order of plasticity at corticostriatal synapses in rats and mice. GABA controls the polarity of STDP in both striatopallidal and striatonigral output neurons. Biophysical simulations and experimental investigations suggest that GABA controls STDP polarity through depolarizing effects at distal dendrites of striatal output neurons by modifying the balance of two calcium sources, NMDARs and voltage-sensitive calcium channels. These findings establish a central role for GABAergic circuits in shaping STDP and suggest that GABA could operate as a Hebbian/anti-Hebbian switch.

Run-Time Interoperability Between Neuronal Network Simulators Based on the MUSIC Framework

MUSIC is a standard API allowing large scale neuron simulators to exchange data within a parallel computer during runtime. A pilot implementation of this API has been released as open source. We provide experiences from the implementation of MUSIC interfaces for two neuronal network simulators of different kinds, NEST and MOOSE. A multi-simulation of a cortico-striatal network model involving both simulators is performed, demonstrating how MUSIC can promote inter-operability between models written for different simulators and how these can be re-used to build a larger model system. Benchmarks show that the MUSIC pilot implementation provides efficient data transfer in a cluster computer with good scaling. We conclude that MUSIC fulfills the design goal that it should be simple to adapt existing simulators to use MUSIC. In addition, since the MUSIC API enforces independence of the applications, the multi-simulation could be built from pluggable component modules without adaptation of the components to each other in terms of simulation time-step or topology of connections between the modules.

Gap Junctions between Striatal Fast-Spiking Interneurons Regulate Spiking Activity and Synchronization as a Function of Cortical Activity

Striatal fast-spiking (FS) interneurons are interconnected by gap junctions into sparsely connected networks. As demonstrated for cortical FS interneurons, these gap junctions in the striatum may cause synchronized spiking, which would increase the influence that FS neurons have on spiking by the striatal medium spiny (MS) neurons. Dysfunction of the basal ganglia is characterized by changes in synchrony or periodicity, thus gap junctions between FS interneurons may modulate synchrony and thereby influence behavior such as reward learning and motor control. To explore the roles of gap junctions on activity and spike synchronization in a striatal FS population, we built a network model of FS interneurons. Each FS connects to 30–40% of its neighbors, as found experimentally, and each FS interneuron in the network is activated by simulated corticostriatal synaptic inputs. Our simulations show that the proportion of synchronous spikes in FS networks with gap junctions increases with increased conductance of the electrical synapse; however, the synchronization effects are moderate for experimentally estimated conductances. Instead, the main tendency is that the presence of gap junctions reduces the total number of spikes generated in response to synaptic inputs in the network. The reduction in spike firing is due to shunting through the gap junctions; which is minimized or absent when the neurons receive coincident inputs. Together these findings suggest that a population of electrically coupled FS interneurons may function collectively as input detectors that are especially sensitive to synchronized synaptic inputs received from the cortex.

Simple cellular and network control principles govern complex patterns of motor behavior

The vertebrate central nervous system is organized in modules that independently execute sophisticated tasks. Such modules are flexibly controlled and operate with a considerable degree of autonomy. One example is locomotion generated by spinal central pattern generator networks (CPGs) that shape the detailed motor output. The level of activity is controlled from brainstem locomotor command centers, which in turn, are under the control of the basal ganglia. By using a biophysically detailed, full-scale computational model of the lamprey CPG (10,000 neurons) and its brainstem/forebrain control, we demonstrate general control principles that can adapt the network to different demands. Forward or backward locomotion and steering can be flexibly controlled by local synaptic effects limited to only the very rostral part of the network. Variability in response properties within each neuronal population is an essential feature and assures a constant phase delay along the cord for different locomotor speeds.

Neural mechanisms potentially contributing to the intersegmental phase lag in lamprey.II. Hemisegmental oscillations produced by mutually coupled excitatory neurons.

Abstract. Most previous models of the spinal central pattern generator (CPG) underlying locomotion in the lamprey have relied on reciprocal inhibition between the left and right side for oscillations to be produced. Here, we have explored the consequences of using self-oscillatory hemisegments. Within a single hemisegment, the oscillations are produced by a network of recurrently coupled excitatory neurons (E neurons) that by themselves are not oscillatory but when coupled together through N-methyl-d-aspartate (NMDA) and α-amino-3-hydroxy-5-methyl-4-isoxazolepropionicacid (AMPA)/kainate transmission can produce oscillations. The bursting mechanism relies on intracellular accumulation of calcium that activates Ca2+-dependent K+. The intracellular calcium is modeled by two different intracellular calcium pools, one of which represents the calcium entry following the action potential, CaAP pool, and the other represents the calcium inflow through the NMDA channels, CaNMDA pool. The Ca2+-dependent K+ activated by these two calcium pools are referred to as KCaAP and KCaNMDA, respectively, and their relative conductances are modulated and increase with the background activation of the network. When changing the background stimulation, the bursting activity in this network can be made to cover a frequency range of 0.5–5.5 Hz with reasonable burst proportions if the adaptation is modulated with the activity. When a chain of such hemisegments are coupled together, a phase lag along the chain can be produced. The local oscillations as well as the phase lag is dependent on the axonal conduction delay as well as the types of excitatory coupling that are assumed, i.e. AMPA/kainate and/or NMDA. When the caudal excitatory projections are extended further than the rostral ones, and assumed to be of approximately equal strength, this kind of network is capable of reproducing several experimental observations such as those occurring during strychnine blockade of the left-right reciprocal inhibition. Addition of reciprocally coupled inhibitory neurons in such a network gives rise to antiphasic activity between the left and right side, but not necessarily to any change of the frequency if the burst proportion of the hemisegmental bursts is well below 50%. Prolongation of the C neuron projection in the rostrocaudal direction restricts the phase lag produced by only the excitatory hemisegmental network by locking together the interburst intervals at different levels of the spinal cord.

Modeling a vertebrate motor system : pattern generation, steering and control of body orientation

The lamprey is one of the few vertebrates in which the neural control system for goal-directed locomotion including steering and control of body orientation is well described at a cellular level. In this report we review the modeling of the central pattern-generating network, which has been carried out based on detailed experimentation. In the same way the modeling of the control system for steering and control of body orientation is reviewed, including neuromechanical simulations and robotic devices.

Neural mechanisms potentially contributing to the intersegmental phase lag in lamprey I. : Segmental oscillations dependent on reciprocal inhibition

Abstract. Factors contributing to the production of a phase lag along chains of oscillatory networks consisting of Hodgkin-Huxley type neurons are analyzed by means of simulations. Simplified network configurations are explored consisting of the basic building blocks of the spinal central pattern generator (CPG) generating swimming in the lamprey. It consists of reciprocally coupled crossed inhibitory C interneurons and ipsilateral excitatory E interneurons that activate C neurons and other E neurons. Oscillatory activity in the model network can, in the simplest case, be produced by a pair of reciprocally coupled C interneurons oscillating through an escape mechanism. Different levels of tonic excitation drive the network over a wide burst frequency range. In this type of network, powerful frequency-regulating factors are the effective inhibition produced by the active side, in combination with the tendency of the inactive side to escape from the inhibition. These two mechanisms can be affected by several factors, e.g. spike frequency adaptation (calcium-dependent K+ channels), N-methyl-D-aspartate membrane properties as well as presence of low-voltage activated calcium channels. A rostrocaudal phase lag can be produced either by extending the contralateral inhibitory projections or the ipsilateral excitatory projections relatively more in the caudal than the rostral direction, since both an increased inhibition and a phasic excitation slow down the receiving network. The phase lag becomes decreased if the length of the intersegmental projections is increased or if the projections are extended symmetrically in both the rostral and the caudal directions. The simulations indicate that the conditions in the ends of an oscillator chain may significantly affect sign, magnitude and constancy of the phase lag. Also, with short and relatively weak intersegmental connections, the network remains robust against perturbations as well as intrinsic frequency differences along the chain. The phase lag (percentage of cycle duration) increases, however, with burst frequency also when the coupling strength is comparatively weak. The results are discussed and compared with previous “phase pulling” models as well as relaxation oscillators.

Integrative neuroscience: linking levels of analyses

If we are to understand how the brain performs different integrated functions in cellular terms, we need both to understand all relevant levels of analysis from the molecular to the behavioural and cognitive levels and to realize an integration of such levels. This is currently a major challenge for neuroscience. Most research, whether dealing with perception, action or learning, focuses on a few levels of organization, for instance the molecular level and brain imaging, and leaves other crucial areas practically untouched. To reach the level of understanding that we desire, a multi-level approach is required in which the different levels link into each other. It is possible to bridge across the different levels for one system, and this has been demonstrated, for example, in the lamprey in generation of goal-directed locomotion. It can be argued that an integrated analysis of any neural system cannot be performed without the aid of a close interaction between experiments and modelling. The dynamic processing within any neural system is such that an intuitive interpretation is rarely sufficient.