Probing the dynamics of identified neurons with a data-driven modeling approach

In controlling animal behavior the nervous system has to perform within the operational limits set by the requirements of each specific behavior. The implications for the corresponding range of suitable network, single neuron, and ion channel properties have remained elusive. In this article we approach the question of how well-constrained properties of neuronal systems may be on the neuronal level. We used large data sets of the activity of isolated invertebrate identified cells and built an accurate conductance-based model for this cell type using customized automated parameter estimation techniques. By direct inspection of the data we found that the variability of the neurons is larger when they are isolated from the circuit than when in the intact system. Furthermore, the responses of the neurons to perturbations appear to be more consistent than their autonomous behavior under stationary conditions. In the developed model, the constraints on different parameters that enforce appropriate model dynamics vary widely from some very tightly controlled parameters to others that are almost arbitrary. The model also allows predictions for the effect of blocking selected ionic currents and to prove that the origin of irregular dynamics in the neuron model is proper chaoticity and that this chaoticity is typical in an appropriate sense. Our results indicate that data driven models are useful tools for the in-depth analysis of neuronal dynamics. The better consistency of responses to perturbations, in the real neurons as well as in the model, suggests a paradigm shift away from measuring autonomous dynamics alone towards protocols of controlled perturbations. Our predictions for the impact of channel blockers on the neuronal dynamics and the proof of chaoticity underscore the wide scope of our approach

A Codimension-2 Bifurcation Controlling Endogenous Bursting Activity and Pulse-Triggered Responses of a Neuron Model

The dynamics of individual neurons are crucial for producing functional activity in neuronal networks. An open question is how temporal characteristics can be controlled in bursting activity and in transient neuronal responses to synaptic input. Bifurcation theory provides a framework to discover generic mechanisms addressing this question. We present a family of mechanisms organized around a global codimension-2 bifurcation. The cornerstone bifurcation is located at the intersection of the border between bursting and spiking and the border between bursting and silence. These borders correspond to the blue sky catastrophe bifurcation and the saddle-node bifurcation on an invariant circle (SNIC) curves, respectively. The cornerstone bifurcation satisfies the conditions for both the blue sky catastrophe and SNIC. The burst duration and interburst interval increase as the inverse of the square root of the difference between the corresponding bifurcation parameter and its bifurcation value. For a given set of burst duration and interburst interval, one can find the parameter values supporting these temporal characteristics. The cornerstone bifurcation also determines the responses of silent and spiking neurons. In a silent neuron with parameters close to the SNIC, a pulse of current triggers a single burst. In a spiking neuron with parameters close to the blue sky catastrophe, a pulse of current temporarily silences the neuron. These responses are stereotypical: the durations of the transient intervals–the duration of the burst and the duration of latency to spiking–are governed by the inverse-square-root laws. The mechanisms described here could be used to coordinate neuromuscular control in central pattern generators. As proof of principle, we construct small networks that control metachronal-wave motor pattern exhibited in locomotion. This pattern is determined by the phase relations of bursting neurons in a simple central pattern generator modeled by a chain of oscillators

Modeling the Influence of Ion Channels on Neuron Dynamics in Drosophila

abstract: Voltage gated ion channels play a major role in determining a neuron's firing behavior, resulting in the specific processing of synaptic input patterns. Drosophila and other invertebrates provide valuable model systems for investigating ion channel kinetics and their impact on firing properties. Despite the increasing importance of Drosophila as a model system, few computational models of its ion channel kinetics have been developed. In this study, experimentally observed biophysical properties of voltage gated ion channels from the fruitfly Drosophila melanogaster are used to develop a minimal, conductance based neuron model. We investigate the impact of the densities of these channels on the excitability of the model neuron. Changing the channel densities reproduces different in situ observed firing patterns and induces a switch from integrator to resonator properties. Further, we analyze the preference to input frequency and how it depends on the channel densities and the resulting bifurcation type the system undergoes. An extension to a three dimensional model demonstrates that the inactivation kinetics of the sodium channels play an important role, allowing for firing patterns with a delayed first spike and subsequent high frequency firing as often observed in invertebrates, without altering the kinetics of the delayed rectifier current.View the article as published at http://journal.frontiersin.org/article/10.3389/fncom.2015.00139/ful

Sparse Gamma Rhythms Arising through Clustering in Adapting Neuronal Networks

Gamma rhythms (30–100 Hz) are an extensively studied synchronous brain state responsible for a number of sensory, memory, and motor processes. Experimental evidence suggests that fast-spiking interneurons are responsible for carrying the high frequency components of the rhythm, while regular-spiking pyramidal neurons fire sparsely. We propose that a combination of spike frequency adaptation and global inhibition may be responsible for this behavior. Excitatory neurons form several clusters that fire every few cycles of the fast oscillation. This is first shown in a detailed biophysical network model and then analyzed thoroughly in an idealized model. We exploit the fact that the timescale of adaptation is much slower than that of the other variables. Singular perturbation theory is used to derive an approximate periodic solution for a single spiking unit. This is then used to predict the relationship between the number of clusters arising spontaneously in the network as it relates to the adaptation time constant. We compare this to a complementary analysis that employs a weak coupling assumption to predict the first Fourier mode to destabilize from the incoherent state of an associated phase model as the external noise is reduced. Both approaches predict the same scaling of cluster number with respect to the adaptation time constant, which is corroborated in numerical simulations of the full system. Thus, we develop several testable predictions regarding the formation and characteristics of gamma rhythms with sparsely firing excitatory neurons

Mechanisms of the Coregulation of Multiple Ionic Currents for the Control of Neuronal Activity

An open question in contemporary neuroscience is how neuromodulators coregulate multiple conductances to maintain functional neuronal activity. Neuromodulators enact changes to properties of biophysical characteristics, such as the maximal conductance or voltage of half-activation of an ionic current, which determine the type and properties of neuronal activity. We apply dynamical systems theory to study the changes to neuronal activity that arise from neuromodulation. Neuromulators can act on multiple targets within a cell. The coregulation of mulitple ionic currents extends the scope of dynamic control on neuronal activity. Different aspects of neuronal activity can be independently controlled by different currents. The coregulation of multiple ionic currents provides precise control over the temporal characteristics of neuronal activity. Compensatory changes in multiple ionic currents could be used to avoid dangerous dynamics or maintain some aspect of neuronal activity. The coregulation of multiple ionic currents can be used as bifurcation control to ensure robust dynamics or expand the range of coexisting regimes. Multiple ionic currents could be involved in increasing the range of dynamic control over neuronal activity. The coregulation of multiple ionic currents in neuromodulation expands the range over which biophysical parameters support functional activity

Homologous Neurons and their Locomotor Functions in Nudibranch Molluscs

These studies compare neurotransmitter localization and the behavioral functions of homologous neurons in nudibranch molluscs to determine the types of changes that might underlie the evolution of species-specific behaviors. Serotonin (5-HT) immunohistochemistry in eleven nudibranch species indicated that certain groups of 5 HT-immunoreactive neurons, such as the Cerebral Serotonergic Posterior (CeSP) cluster, are present in all species. However, the locations and numbers of many other 5 HT-immunoreactive neurons were variable. Thus, particular parts of the serotonergic system have changed during the evolution of nudibranchs. To test whether the functions of homologous neurons are phylogenetically variable, comparisons were made in species with divergent behaviors. In Tritonia diomedea, which crawls and also swims via dorsal-ventral body flexions, the CeSP cluster includes the Dorsal Swim Interneurons (DSIs). It was previously shown that the DSIs are members of the swim central pattern generator (CPG); they are rhythmically active during swimming and, along with their neurotransmitter 5-HT, are necessary and sufficient for swimming. It was also known that the DSIs excite efferent neurons used in crawling. DSI homologues, the CeSP-A neurons, were identified in six species that do not exhibit dorsal-ventral swimming. Many physiological characteristics, including excitation of putative crawling neurons were conserved, but the CeSP A neurons were not rhythmically active in any of the six species. In the lateral flexion swimmer, Melibe leonina, the CeSP-A neurons and 5-HT, were sufficient, but not necessary, for swimming. Thus, homologous neurons, and their neurotransmitter, have functionally diverged in species with different behaviors. Homologous neurons in species with similar behaviors also exhibited functional divergence. Like Melibe, Dendronotus iris is a lateral flexion swimmer. Swim interneuron 1 (Si1) is in the Melibe swim CPG. However, its putative homologue in Dendronotus, the Cerebral Posterior ipsilateral Pedal (CPiP) neuron, was not rhythmically active during swim-like motor patterns, but could initiate such a motor pattern. Together, these studies suggest that neurons have changed their functional relationships to neural circuits during the evolution of species-specific behaviors and have functionally diverged even in species that exhibit similar behaviors