Stochastic Differential Equation Model for Cerebellar Granule Cell Excitability

Neurons in the brain express intrinsic dynamic behavior which is known to be stochastic in nature. A crucial question in building models of neuronal excitability is how to be able to mimic the dynamic behavior of the biological counterpart accurately and how to perform simulations in the fastest possible way. The well-established Hodgkin-Huxley formalism has formed to a large extent the basis for building biophysically and anatomically detailed models of neurons. However, the deterministic Hodgkin-Huxley formalism does not take into account the stochastic behavior of voltage-dependent ion channels. Ion channel stochasticity is shown to be important in adjusting the transmembrane voltage dynamics at or close to the threshold of action potential firing, at the very least in small neurons. In order to achieve a better understanding of the dynamic behavior of a neuron, a new modeling and simulation approach based on stochastic differential equations and Brownian motion is developed. The basis of the work is a deterministic one-compartmental multi-conductance model of the cerebellar granule cell. This model includes six different types of voltage-dependent conductances described by Hodgkin-Huxley formalism and simple calcium dynamics. A new model for the granule cell is developed by incorporating stochasticity inherently present in the ion channel function into the gating variables of conductances. With the new stochastic model, the irregular electrophysiological activity of an in vitro granule cell is reproduced accurately, with the same parameter values for which the membrane potential of the original deterministic model exhibits regular behavior. The irregular electrophysiological activity includes experimentally observed random subthreshold oscillations, occasional spontaneous spikes, and clusters of action potentials. As a conclusion, the new stochastic differential equation model of the cerebellar granule cell excitability is found to expand the range of dynamics in comparison to the original deterministic model. Inclusion of stochastic elements in the operation of voltage-dependent conductances should thus be emphasized more in modeling the dynamic behavior of small neurons. Furthermore, the presented approach is valuable in providing faster computation times compared to the Markov chain type of modeling approaches and more sophisticated theoretical analysis tools compared to previously presented stochastic modeling approaches

Validating Stochastic Models: Invariance Criteria for Systems of Stochastic Differential Equations and the Selection of a Stochastic Hodgkin-Huxley Type Model

In recent years, many difficulties appeared when taking into account the inherent stochastic behavior of neurons and voltage-dependent ion channels in Hodgking-Huxley type models. In particular, an open problem for a stochastic model of cerebellar granule cell excitability was to ensure that the values of the gating variables remain within the unit interval. In this paper, we provide an answer to this modeling issue and obtain a class of viable stochastic models. We select the stochastic models thanks to a general criterion for the flow invariance of rectangular subsets under systems of stochastic differential equations. We formulate explicit necessary and sufficient conditions, that are valid for both, It\^o's and Stratonovich's interpretation of stochastic differential equations, improving a previous result obtained by A. Milian [A.Milian, Coll. Math. 1995] in the It\^o case. These invariance criteria allow to validate stochastic models in many applications. To illustrate our results we present numerical simulations for a stochastic Hodgkin-Huxley model

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 effect of desflurane on neuronal communication at a central synapse

Although general anesthetics are thought to modify critical neuronal functions, their impact on neuronal communication has been poorly examined. We have investigated the effect induced by desflurane, a clinically used general anesthetic, on information transfer at the synapse between mossy fibers and granule cells of cerebellum, where this analysis can be carried out extensively. Mutual information values were assessed by measuring the variability of postsynaptic output in relationship to the variability of a given set of presynaptic inputs. Desflurane synchronized granule cell firing and reduced mutual information in response to physiologically relevant mossy fibers patterns. The decrease in spike variability was due to an increased postsynaptic membrane excitability, which made granule cells more prone to elicit action potentials, and to a strengthened synaptic inhibition, which markedly hampered membrane depolarization. These concomitant actions on granule cells firing indicate that desflurane re-shapes the transfer of information between neurons by providing a less informative neurotransmission rather than completely silencing neuronal activity

Computational Reconstruction of Pacemaking and Intrinsic Electroresponsiveness in Cerebellar Golgi Cells

The Golgi cells have been recently shown to beat regularly in vitro (Forti et al., 2006. J. Physiol. 574, 711–729). Four main currents were shown to be involved, namely a persistent sodium current (I Na-p), an h current (I h), an SK-type calcium-dependent potassium current (I K-AHP), and a slow M-like potassium current (I K-slow). These ionic currents could take part, together with others, also to different aspects of neuronal excitability like responses to depolarizing and hyperpolarizing current injection. However, the ionic mechanisms and their interactions remained largely hypothetical. In this work, we have investigated the mechanisms of Golgi cell excitability by developing a computational model. The model predicts that pacemaking is sustained by subthreshold oscillations tightly coupled to spikes. I Na-p and I K-slow emerged as the critical determinants of oscillations. I h also played a role by setting the oscillatory mechanism into the appropriate membrane potential range. I K-AHP, though taking part to the oscillation, appeared primarily involved in regulating the ISI following spikes. The combination with other currents, in particular a resurgent sodium current (I Na-r) and an A-current (I K-A), allowed a precise regulation of response frequency and delay. These results provide a coherent reconstruction of the ionic mechanisms determining Golgi cell intrinsic electroresponsiveness and suggests important implications for cerebellar signal processing, which will be fully developed in a companion paper (Solinas et al., 2008. Front. Neurosci. 2:4)

Functional properties and pharmacology of extrasynaptic GABA-A receptors

The “ambient” GABA that is present in the extracellular space surrounding all neurons of the brain is believed to be capable of persistently activating high‐affinity extrasynaptic GABA-A receptors to generate a tonic membrane conductance. This generates a form of shunting inhibition that is capable of influencing cellular and network excitability. Extrasynaptic δ subunit-containing GABA-A receptors are known to generate this form of tonic inhibition in a number of defined brain regions and these are emerging as important clinical drug targets for the treatment of a number of neurological conditions. This thesis examines the functional and pharmacological properties of recombinant and native GABA-A receptors that allow them to function as ambient GABA detectors. Surprisingly, the data presented in this Thesis shows that the behaviour of these extrasynaptic GABA-A receptor populations is dramatically influenced by the steady-state GABA concentration they experience. For example, recombinant α4βδ and α6βδ receptor populations are shown to exhibit profound levels of desensitization in the presence of low ambient GABA levels that will limit their ability to respond to changes in ambient GABA. We also find that the action of certain sedative/hypnotic drugs on extrasynaptic GABA-A receptors expressed in cerebellar granule neurons is critically dependent upon the concentration of ambient GABA. For example, we show that the intravenous anaesthetic propofol will only enhance tonic inhibition when ambient GABA levels are below 100 nM. Similarly, we show that the GABA-A receptor agonist Gaboxadol is not capable of enhancing tonic inhibition when ambient GABA levels are high. In contrast to the behaviour of drugs like propofol and Gaboxadol, we find that neurosteroid enhancement of tonic inhibition will occur regardless of ambient GABA levels. This issue will be important when considering therapeutic strategies to target tonic inhibition in the treatment of neurological disorders. Furthermore, we show for the first time, that copper ions can potently block extrasynaptic GABA-A receptors, suggesting that copper may provide a means to selectively block tonic inhibition in the brain, and may even represent a novel source of extrasynaptic GABA-A receptor modulation in vivo

Computational reconstruction of pacemaking and intrinsic electroresponsiveness in cerebellar Golgi cells.

The Golgi cells have been recently shown to beat regularly in vitro (Forti et al., 2006. J. Physiol. 574, 711-729). Four main currents were shown to be involved, namely a persistent sodium current (I(Na-p)), an h current (I(h)), an SK-type calcium-dependent potassium current (I(K-AHP)), and a slow M-like potassium current (I(K-slow)). These ionic currents could take part, together with others, also to different aspects of neuronal excitability like responses to depolarizing and hyperpolarizing current injection. However, the ionic mechanisms and their interactions remained largely hypothetical. In this work, we have investigated the mechanisms of Golgi cell excitability by developing a computational model. The model predicts that pacemaking is sustained by subthreshold oscillations tightly coupled to spikes. I(Na-p) and I(K-slow) emerged as the critical determinants of oscillations. I(h) also played a role by setting the oscillatory mechanism into the appropriate membrane potential range. I(K-AHP), though taking part to the oscillation, appeared primarily involved in regulating the ISI following spikes. The combination with other currents, in particular a resurgent sodium current (I(Na-r)) and an A-current (I(K-A)), allowed a precise regulation of response frequency and delay. These results provide a coherent reconstruction of the ionic mechanisms determining Golgi cell intrinsic electroresponsiveness and suggests important implications for cerebellar signal processing, which will be fully developed in a companion paper (Solinas et al., 2008. Front. Neurosci. 2:4)

Computational Properties of Cerebellar Nucleus Neurons: Effects of Stochastic Ion Channel Gating and Input Location

The function of the nervous system is shaped by the refined integration of synaptic inputs taking place at the single neuron level. Gain modulation is a computational principle that is widely used across the brain, in which the response of a neuronal unit to a set of inputs is affected in a multiplicative fashion by a second set of inputs, but without any effect on its selectivity. The arithmetic operations performed by pyramidal cells in cortical brain areas have been well characterised, along with the underlying mechanisms at the level of networks and cells, for instance background synaptic noise and dendritic saturation. However, in spite of the vast amount of research on the cerebellum and its function, little is known about neuronal computations carried out by its cellular components. A particular area of interest are the cerebellar nuclei, the main output gate of the cerebellum to the brain stem and cortical areas. The aim of this thesis is to contribute to an understanding of the arithmetic operations performed by neurons in the cerebellar nuclei. Focus is placed on two putative determinants, the location of the synaptic input and the presence of channel noise. To analyse the effect of channel noise, the known voltage-gated ion channels of a cerebellar nucleus neuron model are translated to stochastic Markov formalisms and their electrophysiologial behaviour is compared to their deterministic Hodgkin-Huxley counterparts. The findings demonstrate that in most cases, the behaviour of stochastic channels matches the reference deterministic models, with the notable exception of voltage-gated channels with fast kinetics. Two potential explanations are suggested for this discrepancy. Firstly, channels with fast kinetics are strongly affected by the artefactual loss of gating events in the simulation that is caused by the use of a finite-length time step. While this effect can be mitigated, in part, by using very small time steps, the second source of simulation artefacts is the rectification of the distribution of open channels, when channel kinetics characteristics allow the generation of a window current, with an temporal-averaged equilibrium close to zero. Further, stochastic gating is implemented in a realistic cerebellar nucleus neuronal model. The resulting stochastic model exhibits probabilistic spiking and a similar output rate as the corresponding deterministic cerebellar nucleus neuronal model. However, the outcomes of this thesis indicate the computational properties of the cerebellar nucleus neuronal model are independent of the presence of ion channel noise. The main result of this thesis is that the synaptic input location determines the single neuron computational properties, both in the cerebellar nucleus and layer Vb pyramidal neuronal models. The extent of multiplication increases systematically with the distance from the soma, for the cerebellar nucleus, but not for the layer Vb pyramidal neuron, where it is smaller than it would be expected for the distance from the soma. For both neurons, the underlying mechanism is related to the combined effect of nonlinearities introduced by dendritic saturation and the synaptic input noise. However, while excitatory inputs in the perisomatic areas in the cerebellar nucleus undergo additive operations and the distal areas multiplicative, in the layer Vb pyramidal neuron the integration of the excitatory driving input is always multiplicative. In addition, the change in gain is sensitive to the synchronicity of the excitatory synaptic input in the layer Vb pyramidal neuron, but not in the cerebellar nucleus neuron. These observations indicate that the same gain control mechanism might be utilized in distinct ways, in different computational contexts and across different areas, based on the neuronal type and its function

DEVELOPMENT OF A CEREBELLAR MEAN FIELD MODEL: THE THEORETICAL FRAMEWORK, THE IMPLEMENTATION AND THE FIRST APPLICATION

Brain modeling constantly evolves to improve the accuracy of the simulated brain dynamics with the ambitious aim to build a digital twin of the brain. Specific models tuned on brain regions specific features empower the brain simulations introducing bottom-up physiology properties into data-driven simulators. Despite the cerebellum contains 80 % of the neurons and is deeply involved in a wide range of functions, from sensorimotor to cognitive ones, a specific cerebellar model is still missing. Furthermore, its quasi-crystalline multi-layer circuitry deeply differs from the cerebral cortical one, therefore is hard to imagine a unique general model suitable for the realistic simulation of both cerebellar and cerebral cortex. The present thesis tackles the challenge of developing a specific model for the cerebellum. Specifically, multi-neuron multi-layer mean field (MF) model of the cerebellar network, including Granule Cells, Golgi Cells, Molecular Layer Interneurons, and Purkinje Cells, was implemented, and validated against experimental data and the corresponding spiking neural network microcircuit model. The cerebellar MF model was built using a system of interdependent equations, where the single neuronal populations and topological parameters were captured by neuron-specific inter- dependent Transfer Functions. The model time resolution was optimized using Local Field Potentials recorded experimentally with high-density multielectrode array from acute mouse cerebellar slices. The present MF model satisfactorily captured the average discharge of different microcircuit neuronal populations in response to various input patterns and was able to predict the changes in Purkinje Cells firing patterns occurring in specific behavioral conditions: cortical plasticity mapping, which drives learning in associative tasks, and Molecular Layer Interneurons feed-forward inhibition, which controls Purkinje Cells activity patterns. The cerebellar multi-layer MF model thus provides a computationally efficient tool that will allow to investigate the causal relationship between microscopic neuronal properties and ensemble brain activity in health and pathological conditions. Furthermore, preliminary attempts to simulate a pathological cerebellum were done in the perspective of introducing our multi-layer cerebellar MF model in whole-brain simulators to realize patient-specific treatments, moving ahead towards personalized medicine. Two preliminary works assessed the relevant impact of the cerebellum on whole-brain dynamics and its role in modulating complex responses in causal connected cerebral regions, confirming that a specific model is required to further investigate the cerebellum-on- cerebrum influence. The framework presented in this thesis allows to develop a multi-layer MF model depicting the features of a specific brain region (e.g., cerebellum, basal ganglia), in order to define a general strategy to build up a pool of biology grounded MF models for computationally feasible simulations. Interconnected bottom-up MF models integrated in large-scale simulators would capture specific features of different brain regions, while the applications of a virtual brain would have a substantial impact on the reality ranging from the characterization of neurobiological processes, subject-specific preoperative plans, and development of neuro-prosthetic devices