Search

In their seminal works on squid giant axons, Hodgkin and Huxley approximated the membrane leak current as Ohmic, i.e. linear, since in their preparation, sub-threshold current rectification due to the influence of ionic concentration is negligible. Most studies on mammalian neurons have made the same, largely untested, assumption. Here we show that the membrane time constant and input resistance of mammalian neurons (when other major voltage-sensitive and ligand-gated ionic currents are discounted) varies non-linearly with membrane voltage, following the prediction of a Goldman-Hodgkin-Katz-based passive membrane model. The model predicts that under such conditions, the time constant/input resistance-voltage relationship will linearize if the concentration differences across the cell membrane are reduced. These properties were observed in patch-clamp recordings of cerebellar Purkinje neurons (in the presence of pharmacological blockers of other background ionic currents) and were more prominent in the sub-threshold region of the membrane potential. Model simulations showed that the non-linear leak affects voltage-clamp recordings and reduces temporal summation of excitatory synaptic input. Together, our results demonstrate the importance of trans-membrane ionic concentration in defining the functional properties of the passive membrane in mammalian neurons as well as other excitable cells.

AbstractClimbing fiber (CF) triggered complex spikes (CS) are massive depolarization bursts in the cerebellar Purkinje cell, showing several high frequency spikelet components (±600 Hz). Since its early observations, the CS is known to vary in shape. In this study we describe CS waveforms, extracellularly recorded in awake primates (Macaca mulatta) performing saccades. Every Purkinje cell analyzed showed a range of CS shapes with profoundly different duration and number of spikelets. The initial part of the CS was rather constant but the later part differed greatly, with a pronounced jitter of the last spikelets causing a large variation in total CS duration. Waveforms did not effect the following pause duration in the simple spike (SS) train, nor were SS firing rates predictive of the waveform shapes or vice versa. The waveforms did not differ between experimental conditions nor was there a preferred sequential order of CS shapes throughout the recordings. Instead, part of their variability, the timing jitter of the CS’s last spikelets, strongly correlated with interval length to the preceding CS: shorter CS intervals resulted in later appearance of the last spikelets in the CS burst, and vice versa. A similar phenomenon was observed in rat Purkinje cells recorded in vitro upon repeated extracellular stimulation of CFs at different frequencies in slice experiments. All together these results strongly suggest that the variability in the timing of the last spikelet is due to CS frequency dependent changes in Purkinje cell excitability.

Neuronal morphologies are pivotal for brain functioning: physical overlap between dendrites and axons constrain the circuit topology, and the precise shape and composition of dendrites determine the integration of inputs to produce an output signal. At the same time, morphologies are highly diverse and variant. The variance, presumably, originates from neurons developing in a densely packed brain substrate where they interact (e.g., repulsion or attraction) with other actors in this substrate. However, when studying neurons their context is never part of the analysis and they are treated as if they existed in isolation.Here we argue that to fully understand neuronal morphology and its variance it is important to consider neurons in relation to each other and to other actors in the surrounding brain substrate, i.e., their context. We propose a context-aware computational framework, NeuroMaC, in which large numbers of neurons can be grown simultaneously according to growth rules expressed in terms of interactions between the developing neuron and the surrounding brain substrate.As a proof of principle, we demonstrate that by using NeuroMaC we can generate accurate virtual morphologies of distinct classes both in isolation and as part of neuronal forests. Accuracy is validated against population statistics of experimentally reconstructed morphologies. We show that context-aware generation of neurons can explain characteristics of variation. Indeed, plausible variation is an inherent property of the morphologies generated by context-aware rules. We speculate about the applicability of this framework to investigate morphologies and circuits, to classify healthy and pathological morphologies, and to generate large quantities of morphologies for large-scale modeling.

There is growing interest in understanding calcium dynamics in dendrites, both experimentally and computationally. Many processes influence these dynamics, but in dendrites there is a strong contribution of morphology because the peak calcium levels are strongly determined by the surface to volume ratio of each branch, which is inversely related to branch diameter. In this study we explore the predicted variance of dendritic calcium concentrations due to local changes in dendrite diameter and how this is affected by the modeling approach used. We investigate this in a model of dendritic calcium spiking in different reconstructions of cerebellar Purkinje cells and in morphological analysis of neocortical and hippocampal pyramidal neurons. We report that many published models neglect diameter-dependent effects on calcium concentration and show how to implement this correctly in the NEURON simulator, both for phenomenological pool based models and for implementations using radial 1D diffusion. More detailed modeling requires simulation of 3D diffusion and we demonstrate that this does not dissipate the local concentration variance due to changes of dendritic diameter. In many cases 1D diffusion of models of calcium buffering give a good approximation provided an increased morphological resolution is implemented.

STEPS is a stochastic reaction-diffusion simulation engine that implements a spatial extension of Gillespie’s Stochastic Simulation Algorithm (SSA) in complex tetrahedral geometries. An extensive Python-based interface is provided to STEPS so that it can interact with the large number of scientific packages in Python. However, a gap existed between the interfaces of these packages and the STEPS user interface, where supporting toolkits could reduce the amount of scripting required for research projects. This paper introduces two new supporting toolkits that support geometry preparation and visualization for STEPS simulations.

The complex spike (CS) in cerebellar Purkinje Cells (PC) is not an all-or-nothing phenomena as originally proposed, but shows variability depending on the spiking behavior of the Inferior Olive and intrinsic variability in the number and shape of spikelets. The spikelets are repolarised by a sole channel, Kv3.3b, which has been proposed to undergo developmental changes during the postnatal PC maturation. We address here the regulation of the intrinsic CS variability by the expression of inactivating Kv3.3 channels in PCs by combining patch-clamp recordings and single-cell PCR methods on the same neurons, using a technique that we recently optimized to correlate single cell transcription levels with membrane ion channel electrophysiology. We show that while the inactivating TEA sensitive Kv3.3 current peak intensity increases with postnatal age, the channel density does not, arguing against postnatal developmental changes of Kv3.3b expression. Real time PCR of Kv3.3b showed a high variability from cell to cell, correlated with the Kv3.3 current density and suggesting that there are no mechanisms regulating these currents beyond the mRNA pool. We show a significant correlation between normalized quantity of Kv3.3b mRNA and both the number of CS spikelets and their rate of voltage fluctuation, linking the intrinsic CS shape directly to the Kv3.3b mRNA pool. Comparing the observed cell-to-cell variance with studies on transcriptional noise suggests that fluctuations of the Kv3.3b mRNA pool are possibly not regulated but represent merely transcriptional noise, resulting in intrinsic variability of the CS.

We describe a novel method for calculating the quasi-static electrical potential on tetrahedral meshes, which we call E-Field. The E-Field method is implemented in STEPS, which performs stochastic spatial reaction-diffusion computations in tetrahedral-based cellular geometry reconstructions. This provides a level of integration between electrical excitability and spatial molecular dynamics in realistic cellular morphology not previously achievable. Deterministic solutions are also possible. By performing the Rallpack tests we demonstrate the accuracy of the E-Field method. Efficient node ordering is an important practical consideration, and we find that a breadth-first search provides the best solutions, although principal axis ordering suffices for some geometries. We discuss potential applications and possible future directions, and predict that the E-Field implementation in STEPS will play an important role in the future of multiscale neural simulations.