Daniel Keller

Reconstruction and Simulation of Neocortical Microcircuitry

UNLABELLEDnWe present a first-draft digital reconstruction of the microcircuitry of somatosensory cortex of juvenile rat. The reconstruction uses cellular and synaptic organizing principles to algorithmically reconstruct detailed anatomy and physiology from sparse experimental data. An objective anatomical method defines a neocortical volume of 0.29 ± 0.01 mm(3) containing ~31,000 neurons, and patch-clamp studies identify 55 layer-specific morphological and 207 morpho-electrical neuron subtypes. When digitally reconstructed neurons are positioned in the volume and synapse formation is restricted to biological bouton densities and numbers of synapses per connection, their overlapping arbors form ~8 million connections with ~37 million synapses. Simulations reproduce an array of in vitro and in vivo experiments without parameter tuning. Additionally, we find a spectrum of network states with a sharp transition from synchronous to asynchronous activity, modulated by physiological mechanisms. The spectrum of network states, dynamically reconfigured around this transition, supports diverse information processing strategies.nnnPAPERCLIPnVIDEO ABSTRACT.

The neocortical microcircuit collaboration portal: a resource for rat somatosensory cortex.

We have established a multi-constraint, data-driven process to digitally reconstruct, and simulate prototypical neocortical microcircuitry, using sparse experimental data. We applied this process to reconstruct the microcircuitry of the somatosensory cortex in juvenile rat at the cellular and synaptic levels. The resulting reconstruction is broadly consistent with current knowledge about the neocortical microcircuit and provides an array of predictions on its structure and function. To engage the community in exploring, challenging, and refining the reconstruction, we have developed a collaborative, internet-accessible facility-the Neocortical Microcircuit Collaboration portal (NMC portal; https://bbp.epfl.ch/nmc-portal). The NMC portal allows users to access the experimental data used in the reconstruction process, download cellular and synaptic models, and analyze the predicted properties of the microcircuit: six layers, similar to 31,000 neurons, 55 morphological types, 11 electrical types, 207 morpho-electrical types, 1941 unique synaptic connection types between neurons of specific morphological types, predicted properties for the anatomy and physiology of similar to 40 million intrinsic synapses. It also provides data supporting comparison of the anatomy and physiology of the reconstructed microcircuit against results in the literature. The portal aims to catalyzee consensus on the cellular and synaptic organization of neocortical microcircuitry (ion channel, neuron and synapse types and distributions, connectivity, etc.). Community feedback will contribute to refined versions of the reconstruction to be released periodically. We consider that the reconstructions and the simulations they enable represent a major step in the development of in silica neuroscience.

An Exclusion Zone for Ca2+ Channels around Docked Vesicles Explains Release Control by Multiple Channels at a CNS Synapse

The spatial arrangement of Ca2+ channels and vesicles remains unknown for most CNS synapses, despite of the crucial importance of this geometrical parameter for the Ca2+ control of transmitter release. At a large model synapse, the calyx of Held, transmitter release is controlled by several Ca2+ channels in a domain overlap mode, at least in young animals. To study the geometrical constraints of Ca2+ channel placement in domain overlap control of release, we used stochastic MCell modelling, at active zones for which the position of docked vesicles was derived from electron microscopy (EM). We found that random placement of Ca2+ channels was unable to produce high slope values between release and presynaptic Ca2+ entry, a hallmark of domain overlap, and yielded excessively large release probabilities. The simple assumption that Ca2+ channels can be located anywhere at active zones, except below a critical distance of ~ 30 nm away from docked vesicles (exclusion zone), rescued high slope values and low release probabilities. Alternatively, high slope values can also be obtained by placing all Ca2+ channels into a single supercluster, which however results in significantly higher heterogeneity of release probabilities. We also show experimentally that high slope values, and the sensitivity to the slow Ca2+ chelator EGTA-AM, are maintained with developmental maturation of the calyx synapse. Taken together, domain overlap control of release represents a highly organized active zone architecture in which Ca2+ channels must obey a certain distance to docked vesicles. Furthermore, domain overlap can be employed by near-mature, fast-releasing synapses.

Effect of Ionic Diffusion on Extracellular Potentials in Neural Tissue

Recorded potentials in the extracellular space (ECS) of the brain is a standard measure of population activity in neural tissue. Computational models that simulate the relationship between the ECS potential and its underlying neurophysiological processes are commonly used in the interpretation of such measurements. Standard methods, such as volume-conductor theory and current-source density theory, assume that diffusion has a negligible effect on the ECS potential, at least in the range of frequencies picked up by most recording systems. This assumption remains to be verified. We here present a hybrid simulation framework that accounts for diffusive effects on the ECS potential. The framework uses (1) the NEURON simulator to compute the activity and ionic output currents from multicompartmental neuron models, and (2) the electrodiffusive Kirchhoff-Nernst-Planck framework to simulate the resulting dynamics of the potential and ion concentrations in the ECS, accounting for the effect of electrical migration as well as diffusion. Using this framework, we explore the effect that ECS diffusion has on the electrical potential surrounding a small population of 10 pyramidal neurons. The neural model was tuned so that simulations over ∼100 seconds of biological time led to shifts in ECS concentrations by a few millimolars, similar to what has been seen in experiments. By comparing simulations where ECS diffusion was absent with simulations where ECS diffusion was included, we made the following key findings: (i) ECS diffusion shifted the local potential by up to ∼0.2 mV. (ii) The power spectral density (PSD) of the diffusion-evoked potential shifts followed a 1/f2 power law. (iii) Diffusion effects dominated the PSD of the ECS potential for frequencies up to several hertz. In scenarios with large, but physiologically realistic ECS concentration gradients, diffusion was thus found to affect the ECS potential well within the frequency range picked up in experimental recordings.

An Alien Divalent Ion Reveals a Major Role for Ca2+ Buffering in Controlling Slow Transmitter Release

Ca2+-dependent transmitter release occurs in a fast and in a slow phase, but the differential roles of Ca2+ buffers and Ca2+ sensors in shaping release kinetics are still controversial. Replacing extracellular Ca2+ by Sr2+ causes decreased fast release but enhanced slow release at many synapses. Here, we established presynaptic Sr2+ uncaging and made quantitative Sr2+- and Ca2+-imaging experiments at the mouse calyx of Held synapse, to reveal the interplay between Ca2+ sensors and Ca2+ buffers in the control of fast and slow release. We show that Sr2+ activates the fast, Synaptotagmin-2 (Syt2) sensor for vesicle fusion with sixfold lower affinity but unchanged high cooperativity. Surprisingly, Sr2+ also activates the slow sensor that remains in Syt2 knock-out synapses with a lower efficiency, and Sr2+ was less efficient than Ca2+ in the limit of low concentrations in wild-type synapses. Quantitative imaging experiments show that the buffering capacity of the nerve terminal is markedly lower for Sr2+ than for Ca2+ (∼5-fold). This, together with an enhanced Sr2+ permeation through presynaptic Ca2+ channels (∼2-fold), admits a drastically higher spatially averaged Sr2+ transient compared with Ca2+. Together, despite the lower affinity of Sr2+ at the fast and slow sensors, the massively higher amplitudes of spatially averaged Sr2+ transients explain the enhanced late release. This also allows us to conclude that Ca2+ buffering normally controls late release and prevents the activation of the fast release sensor by residual Ca2+.

Network-timing-dependent plasticity.

Bursts of activity in networks of neurons are thought to convey salient information and drive synaptic plasticity. Here we report that network bursts also exert a profound effect on Spike-Timing-Dependent Plasticity (STDP). In acute slices of juvenile rat somatosensory cortex we paired a network burst, which alone induced long-term depression (LTD), with STDP-induced long-term potentiation (LTP) and LTD. We observed that STDP-induced LTP was either unaffected, blocked or flipped into LTD by the network burst, and that STDP-induced LTD was either saturated or flipped into LTP, depending on the relative timing of the network burst with respect to spike coincidences of the STDP event. We hypothesized that network bursts flip STDP-induced LTP to LTD by depleting resources needed for LTP and therefore developed a resource-dependent STDP learning rule. In a model neural network under the influence of the proposed resource-dependent STDP rule, we found that excitatory synaptic coupling was homeostatically regulated to produce power law distributed burst amplitudes reflecting self-organized criticality, a state that ensures optimal information coding.

Norepinephrine stimulates glycogenolysis in astrocytes to fuel neurons with lactate

The mechanism of rapid energy supply to the brain, especially to accommodate the heightened metabolic activity of excited states, is not well-understood. We explored the role of glycogen as a fuel source for neuromodulation using the noradrenergic stimulation of glia in a computational model of the neural-glial-vasculature ensemble (NGV). The detection of norepinephrine (NE) by the astrocyte and the coupled cAMP signal are rapid and largely insensitive to the distance of the locus coeruleus projection release sites from the glia, implying a diminished impact for volume transmission in high affinity receptor transduction systems. Glucosyl-conjugated units liberated from glial glycogen by NE-elicited cAMP second messenger transduction winds sequentially through the glycolytic cascade, generating robust increases in NADH and ATP before pyruvate is finally transformed into lactate. This astrocytic lactate is rapidly exported by monocarboxylate transporters to the associated neuron, demonstrating that the astrocyte-to-neuron lactate shuttle activated by glycogenolysis is a likely fuel source for neuromodulation and enhanced neural activity. Altogether, the energy supply for both astrocytes and neurons can be supplied rapidly by glycogenolysis upon neuromodulatory stimulus.

Comparison of discrete- and continuous-state stochastic methods to model neuronal signal transduction

Several stochastic methods have been developed for the simulation of biochemical reactions. The best known stochastic reaction method is the Gillespie stochastic simulation algorithm which, in this study, is compared to two types of stochastic differential equation models. As a test case, we use a neuronal signal transduction network of 110 reactions and 63 chemical species. We concentrate on showing when stochastic methods are especially needed and how distributions from different stochastic methods differ. We conclude that stochastic differential equations are not suitable for use in volumes on the order of spines, and that even in dendritic and soma regions a portion of the chemical species will exhibit negative excursions. For this reason, a hybrid deterministic and stochastic method may be most applicable to the larger volumes while the Gillespie stochastic simulation algorithm is needed for spines and subspinal volumes. In addition, a grid computing solution is needed for larger volumes to reduce the computation time to tractable levels.

A Process for Digitizing and Simulating Biologically Realistic Oligocellular Networks Demonstrated for the Neuro-Glio-Vascular Ensemble

One will not understand the brain without an integrated exploration of structure and function, these attributes being two sides of the same coin: together they form the currency of biological computation. Accordingly, biologically realistic models require the re-creation of the architecture of the cellular components in which biochemical reactions are contained. We describe here a process of reconstructing a functional oligocellular assembly that is responsible for energy supply management in the brain and creating a computational model of the associated biochemical and biophysical processes. The reactions that underwrite thought are both constrained by and take advantage of brain morphologies pertaining to neurons, astrocytes and the blood vessels that deliver oxygen, glucose and other nutrients. Each component of this neuro-glio-vasculature ensemble (NGV) carries-out delegated tasks, as the dynamics of this system provide for each cell-type its own energy requirements while including mechanisms that allow cooperative energy transfers. Our process for recreating the ultrastructure of cellular components and modeling the reactions that describe energy flow uses an amalgam of state-of the-art techniques, including digital reconstructions of electron micrographs, advanced data analysis tools, computational simulations and in silico visualization software. While we demonstrate this process with the NGV, it is equally well adapted to any cellular system for integrating multimodal cellular data in a coherent framework.

Cell Densities in the Mouse Brain: A Systematic Review

The mouse brain is the most extensively studied brain of all species. We performed an exhaustive review of the literature to establish our current state of knowledge on cell numbers in mouse brain regions, arguably the most fundamental property to measure when attempting to understand a brain. The synthesized information, collected in one place, can be used by both theorists and experimentalists. Although for commonly-studied regions cell densities could be obtained for principal cell types, overall we know very little about how many cells are present in most brain regions and even less about cell-type specific densities. There is also substantial variation in cell density values obtained from different sources. This suggests that we need a new approach to obtain cell density datasets for the mouse brain.