On the Firing Rate Dependency of the Phase Response Curve of Rat Purkinje Neurons In Vitro

Synchronous spiking during cerebellar tasks has been observed across Purkinje cells: however, little is known about the intrinsic cellular mechanisms responsible for its initiation, cessation and stability. The Phase Response Curve (PRC), a simple input-output characterization of single cells, can provide insights into individual and collective properties of neurons and networks, by quantifying the impact of an infinitesimal depolarizing current pulse on the time of occurrence of subsequent action potentials, while a neuron is firing tonically. Recently, the PRC theory applied to cerebellar Purkinje cells revealed that these behave as phase-independent integrators at low firing rates, and switch to a phase-dependent mode at high rates. Given the implications for computation and information processing in the cerebellum and the possible role of synchrony in the communication with its post-synaptic targets, we further explored the firing rate dependency of the PRC in Purkinje cells. We isolated key factors for the experimental estimation of the PRC and developed a closed-loop approach to reliably compute the PRC across diverse firing rates in the same cell. Our results show unambiguously that the PRC of individual Purkinje cells is firing rate dependent and that it smoothly transitions from phase independent integrator to a phase dependent mode. Using computational models we show that neither channel noise nor a realistic cell morphology are responsible for the rate dependent shift in the phase response curve

Human Pluripotent Stem-Cell-Derived Cortical Neurons Integrate Functionally into the Lesioned Adult Murine Visual Cortex in an Area-Specific Way

The transplantation of pluripotent stem-cell-derived neurons constitutes a promising avenue for the treatment of several brain diseases. However, their potential for the repair of the cerebral cortex remains unclear, given its complexity and neuronal diversity. Here, we show that human visual cortical cells differentiated from embryonic stem cells can be transplanted and can integrate successfully into the lesioned mouse adult visual cortex. The transplanted human neurons expressed the appropriate repertoire of markers of six cortical layers, projected axons to specific visual cortical targets, and were synaptically active within the adult brain. Moreover, transplant maturation and integration were much less efficient following transplantation into the lesioned motor cortex, as previously observed for transplanted mouse cortical neurons. These data constitute an important milestone for the potential use of human PSC-derived cortical cells for the reassembly of cortical circuits and emphasize the importance of cortical areal identity for successful transplantation. Espuny-Camacho et al. show that transplanted ESC-derived human cortical neurons integrate functionally into the lesioned adult mouse brain. Transplanted neurons display visual cortical identity and show specific restoration of damaged cortical pathways following transplantation into the visual but not the motor cortex, suggesting the importance of areal-identity match for successful cortical repair

Simple, Fast and Accurate Implementation of the Diffusion Approximation Algorithm for Stochastic Ion Channels with Multiple States

The phenomena that emerge from the interaction of the stochastic opening and closing of ion channels (channel noise) with the non-linear neural dynamics are essential to our understanding of the operation of the nervous system. The effects that channel noise can have on neural dynamics are generally studied using numerical simulations of stochastic models. Algorithms based on discrete Markov Chains (MC) seem to be the most reliable and trustworthy, but even optimized algorithms come with a non-negligible computational cost. Diffusion Approximation (DA) methods use Stochastic Differential Equations (SDE) to approximate the behavior of a number of MCs, considerably speeding up simulation times. However, model comparisons have suggested that DA methods did not lead to the same results as in MC modeling in terms of channel noise statistics and effects on excitability. Recently, it was shown that the difference arose because MCs were modeled with coupled activation subunits, while the DA was modeled using uncoupled activation subunits. Implementations of DA with coupled subunits, in the context of a specific kinetic scheme, yielded similar results to MC. However, it remained unclear how to generalize these implementations to different kinetic schemes, or whether they were faster than MC algorithms. Additionally, a steady state approximation was used for the stochastic terms, which, as we show here, can introduce significant inaccuracies. We derived the SDE explicitly for any given ion channel kinetic scheme. The resulting generic equations were surprisingly simple and interpretable - allowing an easy and efficient DA implementation. The algorithm was tested in a voltage clamp simulation and in two different current clamp simulations, yielding the same results as MC modeling. Also, the simulation efficiency of this DA method demonstrated considerable superiority over MC methods.Comment: 32 text pages, 10 figures, 1 supplementary text + figur

The what and where of adding channel noise to the Hodgkin-Huxley equations

One of the most celebrated successes in computational biology is the Hodgkin-Huxley framework for modeling electrically active cells. This framework, expressed through a set of differential equations, synthesizes the impact of ionic currents on a cell's voltage -- and the highly nonlinear impact of that voltage back on the currents themselves -- into the rapid push and pull of the action potential. Latter studies confirmed that these cellular dynamics are orchestrated by individual ion channels, whose conformational changes regulate the conductance of each ionic current. Thus, kinetic equations familiar from physical chemistry are the natural setting for describing conductances; for small-to-moderate numbers of channels, these will predict fluctuations in conductances and stochasticity in the resulting action potentials. At first glance, the kinetic equations provide a far more complex (and higher-dimensional) description than the original Hodgkin-Huxley equations. This has prompted more than a decade of efforts to capture channel fluctuations with noise terms added to the Hodgkin-Huxley equations. Many of these approaches, while intuitively appealing, produce quantitative errors when compared to kinetic equations; others, as only very recently demonstrated, are both accurate and relatively simple. We review what works, what doesn't, and why, seeking to build a bridge to well-established results for the deterministic Hodgkin-Huxley equations. As such, we hope that this review will speed emerging studies of how channel noise modulates electrophysiological dynamics and function. We supply user-friendly Matlab simulation code of these stochastic versions of the Hodgkin-Huxley equations on the ModelDB website (accession number 138950) and http://www.amath.washington.edu/~etsb/tutorials.html.Comment: 14 pages, 3 figures, review articl

The location of the axon initial segment affects the bandwidth of spike initiation dynamics

The dynamics and the sharp onset of action potential (AP) generation have recently been the subject of intense experimental and theoretical investigations. According to the resistive coupling theory, an electrotonic interplay between the site of AP initiation in the axon and the somato-dendritic load determines the AP waveform. This phenomenon not only alters the shape of AP recorded at the soma, but also determines the dynamics of excitability across a variety of time scales. Supporting this statement, here we generalize a previous numerical study and extend it to the quantification of the input-output gain of the neuronal dynamical response. We consider three classes of multicompartmental mathematical models, ranging from ball-and-stick simplified descriptions of neuronal excitability to 3D-reconstructed biophysical models of excitatory neurons of rodent and human cortical tissue. For each model, we demonstrate that increasing the distance between the axonal site of AP initiation and the soma markedly increases the bandwidth of neuronal response properties. We finally consider the Liquid State Machine paradigm, exploring the impact of altering the site of AP initiation at the level of a neuronal population, and demonstrate that an optimal distance exists to boost the computational performance of the network in a simple classification task. Copyright

High Bandwidth Synaptic Communication and Frequency Tracking in Human Neocortex

Neuronal firing, synaptic transmission, and its plasticity form the building blocks for processing and storage of information in the brain. It is unknown whether adult human synapses are more efficient in transferring information between neurons than rodent synapses. To test this, we recorded from connected pairs of pyramidal neurons in acute brain slices of adult human and mouse temporal cortex and probed the dynamical properties of use-dependent plasticity. We found that human synaptic connections were purely depressing and that they recovered three to four times more swiftly from depression than synapses in rodent neocortex. Thereby, during realistic spike trains, the temporal resolution of synaptic information exchange in human synapses substantially surpasses that in mice. Using information theory, we calculate that information transfer between human pyramidal neurons exceeds that of mouse pyramidal neurons by four to nine times, well into the beta and gamma frequency range. In addition, we found that human principal cells tracked fine temporal features, conveyed in received synaptic inputs, at a wider bandwidth than for rodents. Action potential firing probability was reliably phase-locked to input transients up to 1,000 cycles/s because of a steep onset of action potentials in human pyramidal neurons during spike trains, unlike in rodent neurons. Our data show that, in contrast to the widely held views of limited information transfer in rodent depressing synapses, fast recovering synapses of human neurons can actually transfer substantial amounts of information during spike trains. In addition, human pyramidal neurons are equipped to encode high synaptic information content. Thus, adult human cortical microcircuits relay information at a wider bandwidth than rodent microcircuits

Accurate and Fast Simulation of Channel Noise in Conductance-Based Model Neurons by Diffusion Approximation

Stochastic channel gating is the major source of intrinsic neuronal noise whose functional consequences at the microcircuit- and network-levels have been only partly explored. A systematic study of this channel noise in large ensembles of biophysically detailed model neurons calls for the availability of fast numerical methods. In fact, exact techniques employ the microscopic simulation of the random opening and closing of individual ion channels, usually based on Markov models, whose computational loads are prohibitive for next generation massive computer models of the brain. In this work, we operatively define a procedure for translating any Markov model describing voltage- or ligand-gated membrane ion-conductances into an effective stochastic version, whose computer simulation is efficient, without compromising accuracy. Our approximation is based on an improved Langevin-like approach, which employs stochastic differential equations and no Montecarlo methods. As opposed to an earlier proposal recently debated in the literature, our approximation reproduces accurately the statistical properties of the exact microscopic simulations, under a variety of conditions, from spontaneous to evoked response features. In addition, our method is not restricted to the Hodgkin-Huxley sodium and potassium currents and is general for a variety of voltage- and ligand-gated ion currents. As a by-product, the analysis of the properties emerging in exact Markov schemes by standard probability calculus enables us for the first time to analytically identify the sources of inaccuracy of the previous proposal, while providing solid ground for its modification and improvement we present here

Cell type-specific mechanisms of information transfer in data-driven biophysical models of hippocampal CA3 principal neurons

The transformation of synaptic input into action potential output is a fundamental single-cell computation resulting from the complex interaction of distinct cellular morphology and the unique expression profile of ion channels that define the cellular phenotype. Experimental studies aimed at uncovering the mechanisms of the transfer function have led to important insights, yet are limited in scope by technical feasibility, making biophysical simulations an attractive complementary approach to push the boundaries in our understanding of cellular computation. Here we take a data-driven approach by utilizing high-resolution morphological reconstructions and patch-clamp electrophysiology data together with a multi-objective optimization algorithm to build two populations of biophysically detailed models of murine hippocampal CA3 pyramidal neurons based on the two principal cell types that comprise this region. We evaluated the performance of these models and find that our approach quantitatively matches the cell type-specific firing phenotypes and recapitulate the intrinsic population-level variability in the data. Moreover, we confirm that the conductance values found by the optimization algorithm are consistent with differentially expressed ion channel genes in single-cell transcriptomic data for the two cell types. We then use these models to investigate the cell type-specific biophysical properties involved in the generation of complex-spiking output driven by synaptic input through an information-theoretic treatment of their respective transfer functions. Our simulations identify a host of cell type-specific biophysical mechanisms that define the morpho-functional phenotype to shape the cellular transfer function and place these findings in the context of a role for bursting in CA3 recurrent network synchronization dynamics

A method based on a genetic algorithm to find PWL approximations of multivariate nonlinear functions

In this paper we present a systematic approach to find piecewise-linear approximations of multivariate continuous nonlinear functions, by ensuring a good trade-off between approximation accuracy and model complexity. The proposed (suboptimal) method is based on genetic programming and takes into account the circuit constraints concerning the lower bounds for the size of each domain region (called simplex) where a given nonlinear function is approximated linearly. As a benchmark example, we approximate the well-known Hodgkin- Huxley neuron model

On a piecewise linear approximation of the Hindmarsh-Rose neuron model via genetic algorithms suitable for hardware implementation

This paper is concerned with the approximation of the Hindmarsh and Rose neuron model, which is able to reproduce the main neuronal behaviours, in view of its circuit implementation. The method is based on two main tools: a piecewise-linear approximation technique and bifurcation analysis. The piecewise-linear approximation of the Hindmarsh and Rose model is obtained by solving a mixed-integer optimization problem by a genetic algorithm. The result obtained exhibits a good degree of similarity to the original model, both from a qualitative and a quantitative standpoint