The role of ongoing dendritic oscillations in single-neuron dynamics

How neuronal morphology impacts the synchronisation state of neuronal networks

The biophysical properties of neurons not only affect how information is processed within cells, they can also impact the dynamical states of the network. Specifically, the cellular dynamics of action-potential generation have shown relevance for setting the (de)synchronisation state of the network. The dynamics of tonically spiking neurons typically fall into one of three qualitatively distinct types that arise from distinct mathematical bifurcations of voltage dynamics at the onset of spiking. Accordingly, changes in ion channel composition or even external factors, like temperature, have been demonstrated to switch network behaviour via changes in the spike onset bifurcation and hence its associated dynamical type. A thus far less addressed modulator of neuronal dynamics is cellular morphology. Based on simplified and anatomically realistic mathematical neuron models, we show here that the extent of dendritic arborisation has an influence on the neuronal dynamical spiking type and therefore on the (de)synchronisation state of the network. Specifically, larger dendritic trees prime neuronal dynamics for in-phase-synchronised or splayed-out activity in weakly coupled networks, in contrast to cells with otherwise identical properties yet smaller dendrites. Our biophysical insights hold for generic multicompartmental classes of spiking neuron models (from ball-and-stick-type to anatomically reconstructed models) and establish a connection between neuronal morphology and the susceptibility of neural tissue to synchronisation in health and disease.

Higuchi Fractal Dimension as a Method for Assessing Response to Sound Stimuli in Patients with Diffuse Axonal Brain Injury

The aim of the research was to study the fractal dimension of the EEG signal by the Higuchi's method in patients with diffuse axonal injury (DAI) of the brain.

Materials and Methods

The study was performed in 28 patients with DAI of different severity and 13 sex- and age-matched controls. The Higuchi's method of fractal dimension was used to investigate brain response to sound stimuli of different emotional coloring as well as the features of the EEG signal in the resting state.

Results

The EEG data demonstrated the highest values of fractal dimension in patients with DAI in the resting state. The values of fractal dimension in different emotional states considerably differ both in healthy subjects and in those with DAI. An increase in fractal dimension in response to stimuli occurs predominantly at the frequency of the theta rhythm in the control group and the frequency of the alpha rhythm in the patients with severe DAI.

Conclusion

Higuchi fractal dimension can be used as a complementary diagnostic tool that allows differentiating perception of emotionally significant audio information in patients with brain injury.

The glutamatergic synapse: a complex machinery for information processing

Being the most abundant synaptic type, the glutamatergic synapse is responsible for the larger part of the brain's information processing. Despite the conceptual simplicity of the basic mechanism of synaptic transmission, the glutamatergic synapse shows a large variation in the response to the presynaptic release of the neurotransmitter. This variability is observed not only among different synapses but also in the same single synapse. The synaptic response variability is due to several mechanisms of control of the information transferred among the neurons and suggests that the glutamatergic synapse is not a simple bridge for the transfer of information but plays an important role in its elaboration and management. The control of the synaptic information is operated at pre, post, and extrasynaptic sites in a sort of cooperation between the pre and postsynaptic neurons which also involves the activity of other neurons. The interaction between the different mechanisms of control is extremely complicated and its complete functionality is far from being fully understood. The present review, although not exhaustively, is intended to outline the most important of these mechanisms and their complexity, the understanding of which will be among the most intriguing challenges of future neuroscience.

Episodic Memories: How do the Hippocampus and the Entorhinal Ring Attractors Cooperate to Create Them?

The brain is capable of registering a constellation of events, encountered only once, as an episodic memory that can last for a lifetime. As evidenced by the clinical case of the patient HM, memories preserving their episodic nature still depend on the hippocampal formation, several years after being created, while semantic memories are thought to reside in neocortical areas. The neurobiological substrate of one-time learning and life-long storing in the brain, that must exist at the cellular and circuit level, is still undiscovered. The breakthrough is delayed by the fact that studies jointly investigating the rodent hippocampus and entorhinal cortex are mostly targeted at understanding the spatial aspect of learning. Here, we present the concept of an entorhinal cortical module, termed EPISODE module, that could explain how the representations of different elements constituting episodic memories can be linked together at the stage of encoding. The new model that we propose here reconciles the structural and functional observations made in the entorhinal cortex and explains how the downstream hippocampal processing organizes the representations into meaningful sequences.

A Model of Memory Linking Time to Space

The storage of temporally precise spike patterns can be realized by a single neuron. A spiking neural network (SNN) model is utilized to demonstrate the ability to precisely recall a spike pattern after presenting a single input. We show by using a simulation study that the temporal properties of input patterns can be transformed into spatial patterns of local dendritic spikes. The localization of time-points of spikes is facilitated by phase-shift of the subthreshold membrane potential oscillations (SMO) in the dendritic branches, which modifies their excitability. In reference to the points in time of the arriving input, the dendritic spikes are triggered in different branches. To store spatially distributed patterns, two unsupervised learning mechanisms are utilized. Either synaptic weights to the branches, spatial representation of the temporal input pattern, are enhanced by spike-timing-dependent plasticity (STDP) or the oscillation power of SMOs in spiking branches is increased by dendritic spikes. For retrieval, spike bursts activate stored spatiotemporal patterns in dendritic branches, which reactivate the original somatic spike patterns. The simulation of the prototypical model demonstrates the principle, how linking time to space enables the storage of temporal features of an input. Plausibility, advantages, and some variations of the proposed model are also discussed.

Low-rate firing limit for neurons with axon, soma and dendrites driven by spatially distributed stochastic synapses

Analytical forms for neuronal firing rates are important theoretical tools for the analysis of network states. Since the 1960s, the majority of approaches have treated neurons as being electrically compact and therefore isopotential. These approaches have yielded considerable insight into how single-cell properties affect network activity; however, many neuronal classes, such as cortical pyramidal cells, are electrically extended objects. Calculation of the complex flow of electrical activity driven by stochastic spatio-temporal synaptic input streams in these structures has presented a significant analytical challenge. Here we demonstrate that an extension of the level-crossing method of Rice, previously used for compact cells, provides a general framework for approximating the firing rate of neurons with spatial structure. Even for simple models, the analytical approximations derived demonstrate a surprising richness including: independence of the firing rate to the electrotonic length for certain models, but with a form distinct to the point-like leaky integrate-and-fire model; a non-monotonic dependence of the firing rate on the number of dendrites receiving synaptic drive; a significant effect of the axonal and somatic load on the firing rate; and the role that the trigger position on the axon for spike initiation has on firing properties. The approach necessitates only calculating the mean and variances of the non-thresholded voltage and its rate of change in neuronal structures subject to spatio-temporal synaptic fluctuations. The combination of simplicity and generality promises a framework that can be built upon to incorporate increasing levels of biophysical detail and extend beyond the low-rate firing limit treated in this paper.

Multisynaptic cooperation shapes single glutamatergic synapse response

The activity of thousands of excitatory synapse in the dendritic tree produces variations of membrane potential which, while can produce the spike generation at soma (hillock), can also influence the output of a single glutamatergic synapse. We used a model of synaptic diffusion and EPSP generation to simulate the effect of different number of active synapses on the output of a single one. Our results show that, also in subthreshold conditions, the excitatory dendritic activity can influence several parameters of the single synaptic output such as its amplitude, its time course, the NMDA-component activation and consequently phenomena like STP and LTP.

Resonance Analysis as a Tool for Characterizing Functional Division of Layer 5 Pyramidal Neurons

Evidence suggests that layer 5 pyramidal neurons can be divided into functional zones with unique afferent connectivity and membrane characteristics that allow for post-synaptic integration of feedforward and feedback inputs. To assess the existence of these zones and their interaction, we characterized the resonance properties of a biophysically-realistic compartmental model of a neocortical layer 5 pyramidal neuron. Consistent with recently published theoretical and empirical findings, our model was configured to have a “hot zone” in distal apical dendrite and apical tuft where both high- and low-threshold Ca²⁺ ionic conductances had densities 1–2 orders of magnitude higher than anywhere else in the apical dendrite. We simulated injection of broad spectrum sinusoidal currents with linearly increasing frequency to calculate the input impedance of individual compartments, the transfer impedance between the soma and key compartments within the dendritic tree, and a dimensionless term we introduce called resonance quality. We show that input resonance analysis distinguished at least four distinct zones within the model based on properties of their frequency preferences: basal dendrite which displayed little resonance; soma/proximal apical dendrite which displayed resonance at 5–23 Hz, strongest at 5–10 Hz and hyperpolarized/resting membrane potentials; distal apical dendrite which displayed resonance at 8–19 Hz, strongest at 10 Hz and depolarized membrane potentials; and apical tuft which displayed a weak resonance largely between 8 and 10 Hz across a wide range of membrane potentials. Transfer resonance analysis revealed that changes in subthreshold electrical coupling were found to modulate the transfer resonant frequency of signals transmitted from distal apical dendrite and apical tuft to the soma, which would impact the frequencies that individual neurons are expected to respond to and reinforce. Furthermore, eliminating the hot zone was found to reduce amplification of resonance within the model, which contributes to reduced excitability when perisomatic and distal apical regions receive coincident stimulating current injections. These results indicate that the interactions between different functional zones should be considered in a more complete understanding of neuronal integration. Resonance analysis may therefore be a useful tool for assessing the integration of inputs across the entire neuronal membrane.

Application of Higuchi's fractal dimension from basic to clinical neurophysiology: A review

BACKGROUND AND OBJECTIVE

For more than 20 years, Higuchi's fractal dimension (HFD), as a nonlinear method, has occupied an important place in the analysis of biological signals. The use of HFD has evolved from EEG and single neuron activity analysis to the most recent application in automated assessments of different clinical conditions. Our objective is to provide an updated review of the HFD method applied in basic and clinical neurophysiological research.

METHODS

This article summarizes and critically reviews a broad literature and major findings concerning the applications of HFD for measuring the complexity of neuronal activity during different neurophysiological conditions. The source of information used in this review comes from the PubMed, Scopus, Google Scholar and IEEE Xplore Digital Library databases.

RESULTS

The review process substantiated the significance, advantages and shortcomings of HFD application within all key areas of basic and clinical neurophysiology. Therefore, the paper discusses HFD application alone, combined with other linear or nonlinear measures, or as a part of automated methods for analyzing neurophysiological signals.

CONCLUSIONS

The speed, accuracy and cost of applying the HFD method for research and medical diagnosis make it stand out from the widely used linear methods. However, only a combination of HFD with other nonlinear methods ensures reliable and accurate analysis of a wide range of neurophysiological signals.

Active subthreshold dendritic conductances shape the local field potential

Key points

The local field potential (LFP), the low‐frequency part of extracellular potentials recorded in neural tissue, is often used for probing neural circuit activity. Interpreting the LFP signal is difficult, however.

While the cortical LFP is thought mainly to reflect synaptic inputs onto pyramidal neurons, little is known about the role of the various subthreshold active conductances in shaping the LFP.

By means of biophysical modelling we obtain a comprehensive qualitative understanding of how the LFP generated by a single pyramidal neuron depends on the type and spatial distribution of active subthreshold currents.

For pyramidal neurons, the h‐type channels probably play a key role and can cause a distinct resonance in the LFP power spectrum.

Our results show that the LFP signal can give information about the active properties of neurons and imply that preferred frequencies in the LFP can result from those cellular properties instead of, for example, network dynamics.

Abstract

The main contribution to the local field potential (LFP) is thought to stem from synaptic input to neurons and the ensuing subthreshold dendritic processing. The role of active dendritic conductances in shaping the LFP has received little attention, even though such ion channels are known to affect the subthreshold neuron dynamics. Here we used a modelling approach to investigate the effects of subthreshold dendritic conductances on the LFP. Using a biophysically detailed, experimentally constrained model of a cortical pyramidal neuron, we identified conditions under which subthreshold active conductances are a major factor in shaping the LFP. We found that, in particular, the hyperpolarization‐activated inward current, I_h, can have a sizable effect and cause a resonance in the LFP power spectral density. To get a general, qualitative understanding of how any subthreshold active dendritic conductance and its cellular distribution can affect the LFP, we next performed a systematic study with a simplified model. We found that the effect on the LFP is most pronounced when (1) the synaptic drive to the cell is asymmetrically distributed (i.e. either basal or apical), (2) the active conductances are distributed non‐uniformly with the highest channel densities near the synaptic input and (3) when the LFP is measured at the opposite pole of the cell relative to the synaptic input. In summary, we show that subthreshold active conductances can be strongly reflected in LFP signals, opening up the possibility that the LFP can be used to characterize the properties and cellular distributions of active conductances.

The local field potential (LFP), the low‐frequency part of extracellular potentials recorded in neural tissue, is often used for probing neural circuit activity. Interpreting the LFP signal is difficult, however.

While the cortical LFP is thought mainly to reflect synaptic inputs onto pyramidal neurons, little is known about the role of the various subthreshold active conductances in shaping the LFP.

By means of biophysical modelling we obtain a comprehensive qualitative understanding of how the LFP generated by a single pyramidal neuron depends on the type and spatial distribution of active subthreshold currents.

For pyramidal neurons, the h‐type channels probably play a key role and can cause a distinct resonance in the LFP power spectrum.

Our results show that the LFP signal can give information about the active properties of neurons and imply that preferred frequencies in the LFP can result from those cellular properties instead of, for example, network dynamics.

Grids from bands, or bands from grids? An examination of the effects of single unit contamination on grid cell firing fields

Neural recording technology is improving rapidly, allowing for the detection of spikes from hundreds of cells simultaneously. The limiting step in multielectrode electrophysiology continues to be single cell isolation. However, this step is crucial to the interpretation of data from putative single neurons. We present here, in simulation, an illustration of possibly erroneous conclusions that may be reached when poorly isolated single cell data are analyzed. Grid cells are neurons recorded in rodents, and bats, that spike in equally spaced locations in a hexagonal pattern. One theory states that grid firing patterns arise from a combination of band firing patterns. However, we show here that summing the grid firing patterns of two poorly resolved neurons can result in spurious band-like patterns. Thus, evidence of neurons spiking in band patterns must undergo extreme scrutiny before it is accepted. Toward this aim, we discuss single cell isolation methods and metrics.

Contribution of sublinear and supralinear dendritic integration to neuronal computations

Nonlinear dendritic integration is thought to increase the computational ability of neurons. Most studies focus on how supralinear summation of excitatory synaptic responses arising from clustered inputs within single dendrites result in the enhancement of neuronal firing, enabling simple computations such as feature detection. Recent reports have shown that sublinear summation is also a prominent dendritic operation, extending the range of subthreshold input-output (sI/O) transformations conferred by dendrites. Like supralinear operations, sublinear dendritic operations also increase the repertoire of neuronal computations, but feature extraction requires different synaptic connectivity strategies for each of these operations. In this article we will review the experimental and theoretical findings describing the biophysical determinants of the three primary classes of dendritic operations: linear, sublinear, and supralinear. We then review a Boolean algebra-based analysis of simplified neuron models, which provides insight into how dendritic operations influence neuronal computations. We highlight how neuronal computations are critically dependent on the interplay of dendritic properties (morphology and voltage-gated channel expression), spiking threshold and distribution of synaptic inputs carrying particular sensory features. Finally, we describe how global (scattered) and local (clustered) integration strategies permit the implementation of similar classes of computations, one example being the object feature binding problem.

Grid cell firing patterns may arise from feedback interaction between intrinsic rebound spiking and transverse traveling waves with multiple heading angles

This article presents a model using cellular resonance and rebound properties to model grid cells in medial entorhinal cortex. The model simulates the intrinsic resonance properties of single layer II stellate cells with different frequencies due to the hyperpolarization activated cation current (h current). The stellate cells generate rebound spikes after a delay interval that differs for neurons with different resonance frequency. Stellate cells drive inhibitory interneurons to cause rebound from inhibition in an alternate set of stellate cells that drive interneurons to activate the first set of cells. This allows maintenance of activity with cycle skipping of the spiking of cells that matches recent physiological data on theta cycle skipping. The rebound spiking interacts with subthreshold oscillatory input to stellate cells or interneurons regulated by medial septal input and defined relative to the spatial location coded by neurons. The timing of rebound determines whether the network maintains the activity for the same location or shifts to phases of activity representing a different location. Simulations show that spatial firing patterns similar to grid cells can be generated with a range of different resonance frequencies, indicating how grid cells could be generated with low frequencies present in bats and in mice with knockout of the HCN1 subunit of the h current.

Mechanisms of zero-lag synchronization in cortical motifs

Zero-lag synchronization between distant cortical areas has been observed in a diversity of experimental data sets and between many different regions of the brain. Several computational mechanisms have been proposed to account for such isochronous synchronization in the presence of long conduction delays: Of these, the phenomenon of “dynamical relaying” – a mechanism that relies on a specific network motif – has proven to be the most robust with respect to parameter mismatch and system noise. Surprisingly, despite a contrary belief in the community, the common driving motif is an unreliable means of establishing zero-lag synchrony. Although dynamical relaying has been validated in empirical and computational studies, the deeper dynamical mechanisms and comparison to dynamics on other motifs is lacking. By systematically comparing synchronization on a variety of small motifs, we establish that the presence of a single reciprocally connected pair – a “resonance pair” – plays a crucial role in disambiguating those motifs that foster zero-lag synchrony in the presence of conduction delays (such as dynamical relaying) from those that do not (such as the common driving triad). Remarkably, minor structural changes to the common driving motif that incorporate a reciprocal pair recover robust zero-lag synchrony. The findings are observed in computational models of spiking neurons, populations of spiking neurons and neural mass models, and arise whether the oscillatory systems are periodic, chaotic, noise-free or driven by stochastic inputs. The influence of the resonance pair is also robust to parameter mismatch and asymmetrical time delays amongst the elements of the motif. We call this manner of facilitating zero-lag synchrony resonance-induced synchronization, outline the conditions for its occurrence, and propose that it may be a general mechanism to promote zero-lag synchrony in the brain.

Understanding large-scale neuronal dynamics – and how they relate to the cortical anatomy – is one of the key areas of neuroscience research. Despite a wealth of recent research, the key principles of this relationship have yet to be established. Here we employ computational modeling to study neuronal dynamics on small subgraphs – or motifs – across a hierarchy of spatial scales. We establish a novel organizing principle that we term a “resonance pair” (two mutually coupled nodes), which promotes stable, zero-lag synchrony amongst motif nodes. The bidirectional coupling between a resonance pair acts to mutually adjust their dynamics onto a common and relatively stable synchronized regime, which then propagates and stabilizes the synchronization of other nodes within the motif. Remarkably, we find that this effect can propagate along chains of coupled nodes and hence holds the potential to promote stable zero-lag synchrony in larger sub-networks of cortical systems. Our findings hence suggest a potential unifying account of the existence of zero-lag synchrony, an important phenomenon that may underlie crucial cognitive processes in the brain. Moreover, such pairs of mutually coupled oscillators are found in a wide variety of physical and biological systems suggesting a new, broadly relevant and unifying principle.

Tensegrity, cellular biophysics, and the mechanics of living systems

The recent convergence between physics and biology has led many physicists to enter the fields of cell and developmental biology. One of the most exciting areas of interest has been the emerging field of mechanobiology that centers on how cells control their mechanical properties, and how physical forces regulate cellular biochemical responses, a process that is known as mechanotransduction. In this article, we review the central role that tensegrity (tensional integrity) architecture, which depends on tensile prestress for its mechanical stability, plays in biology. We describe how tensional prestress is a critical governor of cell mechanics and function, and how use of tensegrity by cells contributes to mechanotransduction. Theoretical tensegrity models are also described that predict both quantitative and qualitative behaviors of living cells, and these theoretical descriptions are placed in context of other physical models of the cell. In addition, we describe how tensegrity is used at multiple size scales in the hierarchy of life—from individual molecules to whole living organisms—to both stabilize three-dimensional form and to channel forces from the macroscale to the nanoscale, thereby facilitating mechanochemical conversion at the molecular level.

Neuronal rebound spiking, resonance frequency and theta cycle skipping may contribute to grid cell firing in medial entorhinal cortex

Data show a relationship of cellular resonance and network oscillations in the entorhinal cortex to the spatial periodicity of grid cells. This paper presents a model that simulates the resonance and rebound spiking properties of entorhinal neurons to generate spatial periodicity dependent upon phasic input from medial septum. The model shows that a difference in spatial periodicity can result from a difference in neuronal resonance frequency that replicates data from several experiments. The model also demonstrates a functional role for the phenomenon of theta cycle skipping in the medial entorhinal cortex.

Effect of phase response curve skew on synchronization with and without conduction delays

A central problem in cortical processing including sensory binding and attentional gating is how neurons can synchronize their responses with zero or near-zero time lag. For a spontaneously firing neuron, an input from another neuron can delay or advance the next spike by different amounts depending upon the timing of the input relative to the previous spike. This information constitutes the phase response curve (PRC). We present a simple graphical method for determining the effect of PRC shape on synchronization tendencies and illustrate it using type 1 PRCs, which consist entirely of advances (delays) in response to excitation (inhibition). We obtained the following generic solutions for type 1 PRCs, which include the pulse-coupled leaky integrate and fire model. For pairs with mutual excitation, exact synchrony can be stable for strong coupling because of the stabilizing effect of the causal limit region of the PRC in which an input triggers a spike immediately upon arrival. However, synchrony is unstable for short delays, because delayed inputs arrive during a refractory period and cannot trigger an immediate spike. Right skew destabilizes antiphase and enables modes with time lags that grow as the conduction delay is increased. Therefore, right skew favors near synchrony at short conduction delays and a gradual transition between synchrony and antiphase for pairs coupled by mutual excitation. For pairs with mutual inhibition, zero time lag synchrony is stable for conduction delays ranging from zero to a substantial fraction of the period for pairs. However, for right skew there is a preferred antiphase mode at short delays. In contrast to mutual excitation, left skew destabilizes antiphase for mutual inhibition so that synchrony dominates at short delays as well. These pairwise synchronization tendencies constrain the synchronization properties of neurons embedded in larger networks.

Dendritic signal transmission induced by intracellular charge inhomogeneities

Signal propagation in neuronal dendrites represents the basis for interneuron communication and information processing in the brain. Here we take into account charge inhomogeneities arising in the vicinity of ion channels in cytoplasm and obtain a modified cable equation. We show that charge inhomogeneities acting on a millisecond time scale can lead to the appearance of propagating waves with wavelengths of hundreds of micrometers. They correspond to a certain frequency band predicting the appearance of resonant properties in brain neuron signaling. We also show that membrane potential in spiny dendrites obeys the modified cable equation suggesting a crucial role of the spines in dendritic subthreshold resonance.

Single-neuron criticality optimizes analog dendritic computation

Active dendritic branchlets enable the propagation of dendritic spikes, whose computational functions remain an open question. Here we propose a concrete function to the active channels in large dendritic trees. Modelling the input-output response of large active dendritic arbors subjected to complex spatio-temporal inputs and exhibiting non-stereotyped dendritic spikes, we find that the dendritic arbor can undergo a continuous phase transition from a quiescent to an active state, thereby exhibiting spontaneous and self-sustained localized activity as suggested by experiments. Analogously to the critical brain hypothesis, which states that neuronal networks self-organize near criticality to take advantage of its specific properties, here we propose that neurons with large dendritic arbors optimize their capacity to distinguish incoming stimuli at the critical state. We suggest that "computation at the edge of a phase transition" is more compatible with the view that dendritic arbors perform an analog rather than a digital dendritic computation.

Passive dendrites enable single neurons to compute linearly non-separable functions

Local supra-linear summation of excitatory inputs occurring in pyramidal cell dendrites, the so-called dendritic spikes, results in independent spiking dendritic sub-units, which turn pyramidal neurons into two-layer neural networks capable of computing linearly non-separable functions, such as the exclusive OR. Other neuron classes, such as interneurons, may possess only a few independent dendritic sub-units, or only passive dendrites where input summation is purely sub-linear, and where dendritic sub-units are only saturating. To determine if such neurons can also compute linearly non-separable functions, we enumerate, for a given parameter range, the Boolean functions implementable by a binary neuron model with a linear sub-unit and either a single spiking or a saturating dendritic sub-unit. We then analytically generalize these numerical results to an arbitrary number of non-linear sub-units. First, we show that a single non-linear dendritic sub-unit, in addition to the somatic non-linearity, is sufficient to compute linearly non-separable functions. Second, we analytically prove that, with a sufficient number of saturating dendritic sub-units, a neuron can compute all functions computable with purely excitatory inputs. Third, we show that these linearly non-separable functions can be implemented with at least two strategies: one where a dendritic sub-unit is sufficient to trigger a somatic spike; another where somatic spiking requires the cooperation of multiple dendritic sub-units. We formally prove that implementing the latter architecture is possible with both types of dendritic sub-units whereas the former is only possible with spiking dendrites. Finally, we show how linearly non-separable functions can be computed by a generic two-compartment biophysical model and a realistic neuron model of the cerebellar stellate cell interneuron. Taken together our results demonstrate that passive dendrites are sufficient to enable neurons to compute linearly non-separable functions.

Classical views on single neuron computation treat dendrites as mere collectors of inputs, that is forwarded to the soma for linear summation and causes a spike output if it is sufficiently large. Such a single neuron model can only compute linearly separable input-output functions, representing a small fraction of all possible functions. Recent experimental findings show that in certain pyramidal cells excitatory inputs can be supra-linearly integrated within a dendritic branch, turning this branch into a spiking dendritic sub-unit. Neurons containing many of these dendritic sub-units can compute both linearly separable and linearly non-separable functions. Nevertheless, other neuron types have dendrites which do not spike because the required voltage gated channels are absent. However, these dendrites sub-linearly sum excitatory inputs turning branches into saturating sub-units. We wanted to test if this last type of non-linear summation is sufficient for a single neuron to compute linearly non-separable functions. Using a combination of Boolean algebra and biophysical modeling, we show that a neuron with a single non-linear dendritic sub-unit whether spiking or saturating is able to compute linearly non-separable functions. Thus, in principle, any neuron with a dendritic tree, even passive, can compute linearly non-separable functions.

The robustness of phase-locking in neurons with dendro-dendritic electrical coupling

We examine the effects of dendritic filtering on the existence, stability, and robustness of phase-locked states to heterogeneity and noise in a pair of electrically coupled ball-and-stick neurons with passive dendrites. We use the theory of weakly coupled oscillators and analytically derived filtering properties of the dendritic coupling to systematically explore how the electrotonic length and diameter of dendrites can alter phase-locking. In the case of a fixed value of the coupling conductance (gc) taken from the literature, we find that repeated exchanges in stability between the synchronous and anti-phase states can occur as the electrical coupling becomes more distally located on the dendrites. However, the robustness of the phase-locked states in this case decreases rapidly towards zero as the distance between the electrical coupling and the somata increases. Published estimates of gc are calculated from the experimentally measured coupling coefficient (CC) based on a single-compartment description of a neuron, and therefore may be severe underestimates of gc. With this in mind, we re-examine the stability and robustness of phase-locking using a fixed value of CC, which imposes a limit on the maximum distance the electrical coupling can be located away from the somata. In this case, although the phase-locked states remain robust over the entire range of possible coupling locations, no exchanges in stability with changing coupling position are observed except for a single exchange that occurs in the case of a high somatic firing frequency and a large dendritic radius. Thus, our analysis suggests that multiple exchanges in stability with changing coupling location are unlikely to be observed in real neural systems.

Emergence in the central nervous system

"Emergence" is an idea that has received much attention in consciousness literature, but it is difficult to find characterizations of that concept which are both specific and useful. I will precisely define and characterize a type of epistemic ("weak") emergence and show that it is a property of some neural circuits throughout the CNS, on micro-, meso- and macroscopic levels. I will argue that possession of this property can result in profoundly altered neural dynamics on multiple levels in cortex and other systems. I will first describe emergent neural entities (ENEs) abstractly. I will then show how ENEs function specifically and concretely, and demonstrate some implications of this type of emergence for the CNS.

Voltage dependence of subthreshold resonance frequency in layer II of medial entorhinal cortex

The resonance properties of individual neurons in entorhinal cortex (EC) may contribute to their functional properties in awake, behaving rats. Models propose that entorhinal grid cells could arise from shifts in the intrinsic frequency of neurons caused by changes in membrane potential owing to depolarizing input from neurons coding velocity. To test for potential changes in intrinsic frequency, we measured the resonance properties of neurons at different membrane potentials in neurons in medial and lateral EC. In medial entorhinal neurons, the resonant frequency of individual neurons decreased in a linear manner as the membrane potential was depolarized between -70 and -55 mV. At more hyperpolarized membrane potentials, cells asymptotically approached a maximum resonance frequency. Consistent with the previous studies, near resting potential, the cells of the medial EC possessed a decreasing gradient of resonance frequency along the dorsal to ventral axis, and cells of the lateral EC lacked resonant properties, regardless of membrane potential or position along the medial to lateral axis within lateral EC. Application of 10 μM ZD7288, the H-channel blocker, abolished all resonant properties in MEC cells, and resulted in physiological properties very similar to lateral EC cells. These results on resonant properties show a clear change in frequency response with depolarization that could contribute to the generation of grid cell firing properties in the medial EC.

A homeostatic model of neuronal firing governed by feedback signals from the extracellular matrix

Molecules of the extracellular matrix (ECM) can modulate the efficacy of synaptic transmission and neuronal excitability. These mechanisms are crucial for the homeostatic regulation of neuronal firing over extended timescales. In this study, we introduce a simple mathematical model of neuronal spiking balanced by the influence of the ECM. We consider a neuron receiving random synaptic input in the form of Poisson spike trains and the ECM, which is modeled by a phenomenological variable involved in two feedback mechanisms. One feedback mechanism scales the values of the input synaptic conductance to compensate for changes in firing rate. The second feedback accounts for slow fluctuations of the excitation threshold and depends on the ECM concentration. We show that the ECM-mediated feedback acts as a robust mechanism to provide a homeostatic adjustment of the average firing rate. Interestingly, the activation of feedback mechanisms may lead to a bistability in which two different stable levels of average firing rates can coexist in a spiking network. We discuss the mechanisms of the bistability and how they may be related to memory function.

A model combining oscillations and attractor dynamics for generation of grid cell firing

Different models have been able to account for different features of the data on grid cell firing properties, including the relationship of grid cells to cellular properties and network oscillations. This paper describes a model that combines elements of two major classes of models of grid cells: models using interactions of oscillations and models using attractor dynamics. This model includes a population of units with oscillatory input representing input from the medial septum. These units are termed heading angle cells because their connectivity depends upon heading angle in the environment as well as the spatial phase coded by the cell. These cells project to a population of grid cells. The sum of the heading angle input results in standing waves of circularly symmetric input to the grid cell population. Feedback from the grid cell population increases the activity of subsets of the heading angle cells, resulting in the network settling into activity patterns that resemble the patterns of firing fields in a population of grid cells. The properties of heading angle cells firing as conjunctive grid-by-head-direction cells can shift the grid cell firing according to movement velocity. The pattern of interaction of oscillations requires use of separate populations that fire on alternate cycles of the net theta rhythmic input to grid cells.

Models of grid cell spatial firing published 2005-2011

Since the discovery of grid cells in rat entorhinal cortex, many models of their hexagonally arrayed spatial firing fields have been suggested. We review the models and organize them according to the mechanisms they use to encode position, update the positional code, read it out in the spatial grid pattern, and learn any patterned synaptic connections needed. We mention biological implementations of the models, but focus on the models on Marr’s algorithmic level, where they are not things to individually prove or disprove, but rather are a valuable collection of metaphors of the grid cell system for guiding research that are all likely true to some degree, with each simply emphasizing different aspects of the system. For the convenience of interested researchers, MATLAB implementations of the discussed grid cell models are provided at ModelDB accession 144006 or http://people.bu.edu/zilli/gridmodels.html.

Computational models of grid cells

Grid cells are space-modulated neurons with periodic firing fields. In moving animals, the multiple firing fields of an individual grid cell form a triangular pattern tiling the entire space available to the animal. Collectively, grid cells are thought to provide a context-independent metric representation of the local environment. Since the discovery of grid cells in 2005, a number of models have been proposed to explain the formation of spatially repetitive firing patterns as well as the conversion of these signals to place signals one synapse downstream in the hippocampus. The present article reviews the most recent developments in our understanding of how grid patterns are generated, maintained, and transformed, with particular emphasis on second-generation computational models that have emerged during the past 2-3 years in response to criticism and new data.

Frequency of subthreshold oscillations at different membrane potential voltages in neurons at different anatomical positions on the dorsoventral axis in the rat medial entorhinal cortex

Neurons from layer II of the medial entorhinal cortex show subthreshold membrane potential oscillations (SMPOs) which could contribute to theta-rhythm generation in the entorhinal cortex and to generation of grid cell firing patterns. However, it is unclear whether single neurons have a fixed unique oscillation frequency or whether their frequency varies depending on the mean membrane potential in a cell. We therefore examined the frequency of SMPOs at different membrane potentials in layer II stellate-like cells of the rat medial entorhinal cortex in vitro. Using whole-cell patch recordings, we found that the fluctuations in membrane potential show a broad band of low power frequencies near resting potential that transition to more narrowband oscillation frequencies with depolarization. The transition from broadband to narrowband frequencies depends on the location of the neuron along the dorsoventral axis in the entorhinal cortex, with dorsal neurons transitioning to higher-frequency oscillations relative to ventral neurons transitioning to lower-frequency oscillations. Once SMPOs showed a narrowband frequency, systematic frequency changes were not observed with further depolarization. Using a Hodgkin-Huxley-style model of membrane currents, we show that differences in the influence of depolarization on the frequency of SMPOs at different dorsal to ventral positions could arise from differences in the properties of the h current. The properties of frequency changes in this data are important for evaluating models of the generation of grid cell firing fields with different spacings along the dorsal-to-ventral axis of medial entorhinal cortex.

Effects of conduction delays on the existence and stability of one to one phase locking between two pulse-coupled oscillators

Gamma oscillations can synchronize with near zero phase lag over multiple cortical regions and between hemispheres, and between two distal sites in hippocampal slices. How synchronization can take place over long distances in a stable manner is considered an open question. The phase resetting curve (PRC) keeps track of how much an input advances or delays the next spike, depending upon where in the cycle it is received. We use PRCs under the assumption of pulsatile coupling to derive existence and stability criteria for 1:1 phase-locking that arises via bidirectional pulse coupling of two limit cycle oscillators with a conduction delay of any duration for any 1:1 firing pattern. The coupling can be strong as long as the effect of one input dissipates before the next input is received. We show the form that the generic synchronous and anti-phase solutions take in a system of two identical, identically pulse-coupled oscillators with identical delays. The stability criterion has a simple form that depends only on the slopes of the PRCs at the phases at which inputs are received and on the number of cycles required to complete the delayed feedback loop. The number of cycles required to complete the delayed feedback loop depends upon both the value of the delay and the firing pattern. We successfully tested the predictions of our methods on networks of model neurons. The criteria can easily be extended to include the effect of an input on the cycle after the one in which it is received.

Dendritic slow dynamics enables localized cortical activity to switch between mobile and immobile modes with noisy background input

Mounting lines of evidence suggest the significant computational ability of a single neuron empowered by active dendritic dynamics. This motivates us to study what functionality can be acquired by a network of such neurons. The present paper studies how such rich single-neuron dendritic dynamics affects the network dynamics, a question which has scarcely been specifically studied to date. We simulate neurons with active dendrites networked locally like cortical pyramidal neurons, and find that naturally arising localized activity--called a bump--can be in two distinct modes, mobile or immobile. The mode can be switched back and forth by transient input to the cortical network. Interestingly, this functionality arises only if each neuron is equipped with the observed slow dendritic dynamics and with in vivo-like noisy background input. If the bump activity is considered to indicate a point of attention in the sensory areas or to indicate a representation of memory in the storage areas of the cortex, this would imply that the flexible mode switching would be of great potential use for the brain as an information processing device. We derive these conclusions using a natural extension of the conventional field model, which is defined by combining two distinct fields, one representing the somatic population and the other representing the dendritic population. With this tool, we analyze the spatial distribution of the degree of after-spike adaptation and explain how we can understand the presence of the two distinct modes and switching between the modes. We also discuss the possible functional impact of this mode-switching ability.

Bistability induces episodic spike communication by inhibitory neurons in neuronal networks

Bistability is one of the important features of nonlinear dynamical systems. In neurodynamics, bistability has been found in basic Hodgkin-Huxley equations describing the cell membrane dynamics. When the neuron is clamped near its threshold, the stable rest potential may coexist with the stable limit cycle describing periodic spiking. However, this effect is often neglected in network computations where the neurons are typically reduced to threshold firing units (e.g., integrate-and-fire models). We found that the bistability may induce spike communication by inhibitory coupled neurons in the spiking network. The communication is realized in the form of episodic discharges with synchronous (correlated) spikes during the episodes. A spiking phase map is constructed to describe the synchronization and to estimate basic spike phase locking modes.

Models of neocortical layer 5b pyramidal cells capturing a wide range of dendritic and perisomatic active properties

The thick-tufted layer 5b pyramidal cell extends its dendritic tree to all six layers of the mammalian neocortex and serves as a major building block for the cortical column. L5b pyramidal cells have been the subject of extensive experimental and modeling studies, yet conductance-based models of these cells that faithfully reproduce both their perisomatic Na(+)-spiking behavior as well as key dendritic active properties, including Ca(2+) spikes and back-propagating action potentials, are still lacking. Based on a large body of experimental recordings from both the soma and dendrites of L5b pyramidal cells in adult rats, we characterized key features of the somatic and dendritic firing and quantified their statistics. We used these features to constrain the density of a set of ion channels over the soma and dendritic surface via multi-objective optimization with an evolutionary algorithm, thus generating a set of detailed conductance-based models that faithfully replicate the back-propagating action potential activated Ca(2+) spike firing and the perisomatic firing response to current steps, as well as the experimental variability of the properties. Furthermore, we show a useful way to analyze model parameters with our sets of models, which enabled us to identify some of the mechanisms responsible for the dynamic properties of L5b pyramidal cells as well as mechanisms that are sensitive to morphological changes. This automated framework can be used to develop a database of faithful models for other neuron types. The models we present provide several experimentally-testable predictions and can serve as a powerful tool for theoretical investigations of the contribution of single-cell dynamics to network activity and its computational capabilities.

Phase precession and variable spatial scaling in a periodic attractor map model of medial entorhinal grid cells with realistic after-spike dynamics

We present a model that describes the generation of the spatial (grid fields) and temporal (phase precession) properties of medial entorhinal cortical (MEC) neurons by combining network and intrinsic cellular properties. The model incorporates network architecture derived from earlier attractor map models, and is implemented in 1D for simplicity. Periodic driving of conjunctive (position × head-direction) layer-III MEC cells at theta frequency with intensity proportional to the rat's speed, moves an 'activity bump' forward in network space at a corresponding speed. The addition of prolonged excitatory currents and simple after-spike dynamics resembling those observed in MEC stellate cells (for which new data are presented) accounts for both phase precession and the change in scale of grid fields along the dorso-ventral axis of MEC. Phase precession in the model depends on both synaptic connectivity and intrinsic currents, each of which drive neural spiking either during entry into, or during exit out of a grid field. Thus, the model predicts that the slope of phase precession changes between entry into and exit out of the field. The model also exhibits independent variation in grid spatial period and grid field size, which suggests possible experimental tests of the model.

Role of active dendritic conductances in subthreshold input integration

Dendrites of many types of neurons contain voltage-dependent conductances that are active at subthreshold membrane potentials. To understand the computations neurons perform it is key to understand the role of active dendrites in the subthreshold processing of synaptic inputs. We examine systematically how active dendritic conductances affect the time course of postsynaptic potentials propagating along dendrites, and how they affect the interaction between such signals. Voltage-dependent currents can be classified into two types that have qualitatively different effects on subthreshold input responses: regenerative dendritic currents boost and broaden EPSPs, while restorative currents attenuate and narrow EPSPs. Importantly, the effects of active dendritic currents on EPSP shape increase as the EPSP travels along the dendrite. The effectiveness of active currents in modulating the EPSP shape is determined by their activation time constant: the faster it is, the stronger the effect on EPSP amplitude, while the largest effects on EPSP width occur when it is comparable to the membrane time constant. We finally demonstrate that the two current types can differentially improve precision and robustness of neural computations: restorative currents enhance coincidence detection of dendritic inputs, whereas direction selectivity to sequences of dendritic inputs is enhanced by regenerative dendritic currents.

Coupled noisy spiking neurons as velocity-controlled oscillators in a model of grid cell spatial firing

One of the two primary classes of models of grid cell spatial firing uses interference between oscillators at dynamically modulated frequencies. Generally, these models are presented in terms of idealized oscillators (modeled as sinusoids), which differ from biological oscillators in multiple important ways. Here we show that two more realistic, noisy neural models (Izhikevich's simple model and a biophysical model of an entorhinal cortex stellate cell) can be successfully used as oscillators in a model of this type. When additive noise is included in the models such that uncoupled or sparsely coupled cells show realistic interspike interval variance, both synaptic and gap-junction coupling can synchronize networks of cells to produce comparatively less variable network-level oscillations. We show that the frequency of these oscillatory networks can be controlled sufficiently well to produce stable grid cell spatial firing on the order of at least 2-5 min, despite the high noise level. Our results suggest that the basic principles of oscillatory interference models work with more realistic models of noisy neurons. Nevertheless, a number of simplifications were still made and future work should examine increasingly realistic models.

Neurophysiological and computational principles of cortical rhythms in cognition

Synchronous rhythms represent a core mechanism for sculpting temporal coordination of neural activity in the brain-wide network. This review focuses on oscillations in the cerebral cortex that occur during cognition, in alert behaving conditions. Over the last two decades, experimental and modeling work has made great strides in elucidating the detailed cellular and circuit basis of these rhythms, particularly gamma and theta rhythms. The underlying physiological mechanisms are diverse (ranging from resonance and pacemaker properties of single cells to multiple scenarios for population synchronization and wave propagation), but also exhibit unifying principles. A major conceptual advance was the realization that synaptic inhibition plays a fundamental role in rhythmogenesis, either in an interneuronal network or in a reciprocal excitatory-inhibitory loop. Computational functions of synchronous oscillations in cognition are still a matter of debate among systems neuroscientists, in part because the notion of regular oscillation seems to contradict the common observation that spiking discharges of individual neurons in the cortex are highly stochastic and far from being clocklike. However, recent findings have led to a framework that goes beyond the conventional theory of coupled oscillators and reconciles the apparent dichotomy between irregular single neuron activity and field potential oscillations. From this perspective, a plethora of studies will be reviewed on the involvement of long-distance neuronal coherence in cognitive functions such as multisensory integration, working memory, and selective attention. Finally, implications of abnormal neural synchronization are discussed as they relate to mental disorders like schizophrenia and autism.

Losing phase

In this issue of Neuron, Remme and colleagues examine the biophysics of synchronization between oscillating dendrites and soma. Their findings suggest that oscillators will quickly phase-lock when weakly coupled. These findings are at odds with assumptions of an influential model of grid cell response generation and have implications for grid cell response mechanisms.

Democracy-independence trade-off in oscillating dendrites and its implications for grid cells

Dendritic democracy and independence have been characterized for near-instantaneous processing of synaptic inputs. However, a wide class of neuronal computations requires input integration on long timescales. As a paradigmatic example, entorhinal grid fields have been thought to be generated by the democratic summation of independent dendritic oscillations performing direction-selective path integration. We analyzed how multiple dendritic oscillators embedded in the same neuron integrate inputs separately and determine somatic membrane voltage jointly. We found that the interaction of dendritic oscillations leads to phase locking, which sets an upper limit on the timescale for independent input integration. Factors that increase this timescale also decrease the influence that the dendritic oscillations exert on somatic voltage. In entorhinal stellate cells, interdendritic coupling dominates and causes these cells to act as single oscillators. Our results suggest a fundamental trade-off between local and global processing in dendritic trees integrating ongoing signals.