Experimentally guided modelling of dendritic excitability in rat neocortical pyramidal neurones

A biophysical and statistical modeling paradigm for connecting neural physiology and function

To understand single neuron computation, it is necessary to know how specific physiological parameters affect neural spiking patterns that emerge in response to specific stimuli. Here we present a computational pipeline combining biophysical and statistical models that provides a link between variation in functional ion channel expression and changes in single neuron stimulus encoding. More specifically, we create a mapping from biophysical model parameters to stimulus encoding statistical model parameters. Biophysical models provide mechanistic insight, whereas statistical models can identify associations between spiking patterns and the stimuli they encode. We used public biophysical models of two morphologically and functionally distinct projection neuron cell types: mitral cells (MCs) of the main olfactory bulb, and layer V cortical pyramidal cells (PCs). We first simulated sequences of action potentials according to certain stimuli while scaling individual ion channel conductances. We then fitted point process generalized linear models (PP-GLMs), and we constructed a mapping between the parameters in the two types of models. This framework lets us detect effects on stimulus encoding of changing an ion channel conductance. The computational pipeline combines models across scales and can be applied as a screen of channels, in any cell type of interest, to identify ways that channel properties influence single neuron computation.

MMM - The molecular model of memory

Identifying mechanisms underlying neurons ability to process information including acquisition, storage, and retrieval plays an important role in the understanding of the different types of memory, pathogenesis of many neurological diseases affecting memory and therapeutic target discovery. However, the traditional understanding of the mechanisms of memory associated with the electrical signals having a unique combination of frequency and amplitude does not answer the question how the memories can survive for life-long periods of time, while exposed to synaptic noise. Recent evidence suggests that, apart from neuronal circuits, a diversity of the molecular memory (MM) carriers, are essential for memory performance. The molecular model of memory (MMM) is proposed, according to which each item of incoming information (the elementary memory item - eMI) is encoded by both circuitries, with the unique for a given MI electrical parameters, and also the MM carriers, unique by its molecular composition. While operating as the carriers of incoming information, the MMs, are functioning within the neuron plasma membrane. Inactive (latent) initially, during acquisition each of the eMIs is activated to become a virtual copy of some real fact or events bygone. This activation is accompanied by the considerable remodeling of the MM molecule associated with the resonance effect.

Kinetic and thermodynamic modeling of a voltage-gated sodium channel

Like all biological and chemical reactions, ion channel kinetics are highly sensitive to changes in temperature. Therefore, it is prudent to investigate channel dynamics at physiological temperatures. However, most ion channel investigations are performed at room temperature due to practical considerations, such as recording stability and technical limitations. This problem is especially severe for the fast voltage-gated sodium channel, whose activation kinetics are faster than the time constant of the standard patch-clamp amplifier at physiological temperatures. Thus, biologically detailed simulations of the action potential generation evenly scale the kinetic models of voltage-gated channels acquired at room temperature. To quantitatively study voltage-gated sodium channels' temperature sensitivity, we recorded sodium currents from nucleated patches extracted from the rat's layer five neocortical pyramidal neurons at several temperatures from 13.5 to 30 °C. We use these recordings to model the kinetics of the voltage-gated sodium channel as a function of temperature. We show that the temperature dependence of activation differs from that of inactivation. Furthermore, the data indicate that the sustained current has a different temperature dependence than the fast current. Our kinetic and thermodynamic analysis of the current provided a numerical model spanning the entire temperature range. This model reproduced vital features of channel activation and inactivation. Furthermore, the model also reproduced action potential dependence on temperature. Thus, we provide an essential building block for the generation of biologically detailed models of cortical neurons.

NeuroGPU: Accelerating multi-compartment, biophysically detailed neuron simulations on GPUs

BACKGROUND

The membrane potential of individual neurons depends on a large number of interacting biophysical processes operating on spatial-temporal scales spanning several orders of magnitude. The multi-scale nature of these processes dictates that accurate prediction of membrane potentials in specific neurons requires the utilization of detailed simulations. Unfortunately, constraining parameters within biologically detailed neuron models can be difficult, leading to poor model fits. This obstacle can be overcome partially by numerical optimization or detailed exploration of parameter space. However, these processes, which currently rely on central processing unit (CPU) computation, often incur orders of magnitude increases in computing time for marginal improvements in model behavior. As a result, model quality is often compromised to accommodate compute resources.

NEW METHOD

Here, we present a simulation environment, NeuroGPU, that takes advantage of the inherent parallelized structure of the graphics processing unit (GPU) to accelerate neuronal simulation.

RESULTS & COMPARISON WITH EXISTING METHODS

NeuroGPU can simulate most biologically detailed models 10-200 times faster than NEURON simulation running on a single core and 5 times faster than GPU simulators (CoreNEURON). NeuroGPU is designed for model parameter tuning and best performs when the GPU is fully utilized by running multiple (> 100) instances of the same model with different parameters. When using multiple GPUs, NeuroGPU can reach to a speed-up of 800 fold compared to single core simulations, especially when simulating the same model morphology with different parameters. We demonstrate the power of NeuoGPU through large-scale parameter exploration to reveal the response landscape of a neuron. Finally, we accelerate numerical optimization of biophysically detailed neuron models to achieve highly accurate fitting of models to simulation and experimental data.

CONCLUSIONS

Thus, NeuroGPU is the fastest available platform that enables rapid simulation of multi-compartment, biophysically detailed neuron models on commonly used computing systems accessible by many scientists.

Integrating single-cell transcriptomics and microcircuit computer modeling

Biophysically realistic computer modeling of neuronal microcircuitry has served as a testing ground for hypotheses related to the structure and function of different brain microcircuits. Recent advances in single-cell transcriptomics provide snapshots of a neuron's molecular state and have demonstrated that cell-specific genetic markers engineer the electrophysiological properties of a neuron. Integrating these molecular details with biophysical modeling can allow unprecedented mechanistic insights. In this opinion review, we consider systems biology-based strategies involving statistical deconvolution and gene ontology to integrate the two approaches. We foresee that this integration will infer the nonlinear interactions between the transcriptomically detailed neurons in different brain states. For an initial assessment of these integrative strategies, we recommend testing them on a penetrant phenotype such as epilepsy or a basic organism model such as Caenorhabditis elegans.

Modeling Reveals Human-Rodent Differences in H-Current Kinetics Influencing Resonance in Cortical Layer 5 Neurons

While our understanding of human neurons is often inferred from rodent data, inter-species differences between neurons can be captured by building cellular models specifically from human data. This includes understanding differences at the level of ion channels and their implications for human brain function. Thus, we here present a full spiking, biophysically detailed multi-compartment model of a human layer 5 (L5) cortical pyramidal cell. Model development was primarily based on morphological and electrophysiological data from the same human L5 neuron, avoiding confounds of experimental variability. Focus was placed on describing the behavior of the hyperpolarization-activated cation (h-) channel, given increasing interest in this channel due to its role in pacemaking and differentiating cell types. We ensured that the model exhibited post-inhibitory rebound spiking considering its relationship with the h-current, along with other general spiking characteristics. The model was validated against data not used in its development, which highlighted distinctly slower kinetics of the human h-current relative to the rodent setting. We linked the lack of subthreshold resonance observed in human L5 neurons to these human-specific h-current kinetics. This work shows that it is possible and necessary to build human-specific biophysical neuron models in order to understand human brain dynamics.

Global and Multiplexed Dendritic Computations under In Vivo-like Conditions

Dendrites integrate inputs nonlinearly, but it is unclear how these nonlinearities contribute to the overall input-output transformation of single neurons. We developed statistically principled methods using a hierarchical cascade of linear-nonlinear subunits (hLN) to model the dynamically evolving somatic response of neurons receiving complex, in vivo-like spatiotemporal synaptic input patterns. We used the hLN to predict the somatic membrane potential of an in vivo-validated detailed biophysical model of a L2/3 pyramidal cell. Linear input integration with a single global dendritic nonlinearity achieved above 90% prediction accuracy. A novel hLN motif, input multiplexing into parallel processing channels, could improve predictions as much as conventionally used additional layers of local nonlinearities. We obtained similar results in two other cell types. This approach provides a data-driven characterization of a key component of cortical circuit computations: the input-output transformation of neurons during in vivo-like conditions.

Biophysical Psychiatry-How Computational Neuroscience Can Help to Understand the Complex Mechanisms of Mental Disorders

The brain is the most complex of human organs, and the pathophysiology underlying abnormal brain function in psychiatric disorders is largely unknown. Despite the rapid development of diagnostic tools and treatments in most areas of medicine, our understanding of mental disorders and their treatment has made limited progress during the last decades. While recent advances in genetics and neuroscience have a large potential, the complexity and multidimensionality of the brain processes hinder the discovery of disease mechanisms that would link genetic findings to clinical symptoms and behavior. This applies also to schizophrenia, for which genome-wide association studies have identified a large number of genetic risk loci, spanning hundreds of genes with diverse functionalities. Importantly, the multitude of the associated variants and their prevalence in the healthy population limit the potential of a reductionist functional genetics approach as a stand-alone solution to discover the disease pathology. In this review, we outline the key concepts of a "biophysical psychiatry," an approach that employs large-scale mechanistic, biophysics-founded computational modelling to increase transdisciplinary understanding of the pathophysiology and strive toward robust predictions. We discuss recent scientific advances that allow a synthesis of previously disparate fields of psychiatry, neurophysiology, functional genomics, and computational modelling to tackle open questions regarding the pathophysiology of heritable mental disorders. We argue that the complexity of the increasing amount of genetic data exceeds the capabilities of classical experimental assays and requires computational approaches. Biophysical psychiatry, based on modelling diseased brain networks using existing and future knowledge of basic genetic, biochemical, and functional properties on a single neuron to a microcircuit level, may allow a leap forward in deriving interpretable biomarkers and move the field toward novel treatment options.

Parameter Optimization Using Covariance Matrix Adaptation-Evolutionary Strategy (CMA-ES), an Approach to Investigate Differences in Channel Properties Between Neuron Subtypes

Computational models in neuroscience can be used to predict causal relationships between biological mechanisms in neurons and networks, such as the effect of blocking an ion channel or synaptic connection on neuron activity. Since developing a biophysically realistic, single neuron model is exceedingly difficult, software has been developed for automatically adjusting parameters of computational neuronal models. The ideal optimization software should work with commonly used neural simulation software; thus, we present software which works with models specified in declarative format for the MOOSE simulator. Experimental data can be specified using one of two different file formats. The fitness function is customizable as a weighted combination of feature differences. The optimization itself uses the covariance matrix adaptation-evolutionary strategy, because it is robust in the face of local fluctuations of the fitness function, and deals well with a high-dimensional and discontinuous fitness landscape. We demonstrate the versatility of the software by creating several model examples of each of four types of neurons (two subtypes of spiny projection neurons and two subtypes of globus pallidus neurons) by tuning to current clamp data. Optimizations reached convergence within 1,600-4,000 model evaluations (200-500 generations × population size of 8). Analysis of the parameters of the best fitting models revealed differences between neuron subtypes, which are consistent with prior experimental results. Overall our results suggest that this easy-to-use, automatic approach for finding neuron channel parameters may be applied to current clamp recordings from neurons exhibiting different biochemical markers to help characterize ionic differences between other neuron subtypes.

A stepwise neuron model fitting procedure designed for recordings with high spatial resolution: Application to layer 5 pyramidal cells

BACKGROUND

Recent progress in electrophysiological and optical methods for neuronal recordings provides vast amounts of high-resolution data. In parallel, the development of computer technology has allowed simulation of ever-larger neuronal circuits. A challenge in taking advantage of these developments is the construction of single-cell and network models in a way that faithfully reproduces neuronal biophysics with subcellular level of details while keeping the simulation costs at an acceptable level.

NEW METHOD

In this work, we develop and apply an automated, stepwise method for fitting a neuron model to data with fine spatial resolution, such as that achievable with voltage sensitive dyes (VSDs) and Ca²⁺ imaging.

RESULT

We apply our method to simulated data from layer 5 pyramidal cells (L5PCs) and construct a model with reduced neuronal morphology. We connect the reduced-morphology neurons into a network and validate against simulated data from a high-resolution L5PC network model.

COMPARISON WITH EXISTING METHODS

Our approach combines features from several previously applied model-fitting strategies. The reduced-morphology neuron model obtained using our approach reliably reproduces the membrane-potential dynamics across the dendrites as predicted by the full-morphology model.

CONCLUSIONS

The network models produced using our method are cost-efficient and predict that interconnected L5PCs are able to amplify delta-range oscillatory inputs across a large range of network sizes and topologies, largely due to the medium after hyperpolarization mediated by the Ca²⁺-activated SK current.

Neuroelectric Tuning of Cortical Oscillations by Apical Dendrites in Loop Circuits

Bundles of relatively long apical dendrites dominate the neurons that make up the thickness of the cerebral cortex. It is proposed that a major function of the apical dendrite is to produce sustained oscillations at a specific frequency that can serve as a common timing unit for the processing of information in circuits connected to that apical dendrite. Many layer 5 and 6 pyramidal neurons are connected to thalamic neurons in loop circuits. A model of the apical dendrites of these pyramidal neurons has been used to simulate the electric activity of the apical dendrite. The results of that simulation demonstrated that subthreshold electric pulses in these apical dendrites can be tuned to specific frequencies and also can be fine-tuned to narrow bandwidths of less than one Hertz (1 Hz). Synchronous pulse outputs from the circuit loops containing apical dendrites can tune subthreshold membrane oscillations of neurons they contact. When the pulse outputs are finely tuned, they function as a local “clock,” which enables the contacted neurons to synchronously communicate with each other. Thus, a shared tuning frequency can select neurons for membership in a circuit. Unlike layer 6 apical dendrites, layer 5 apical dendrites can produce burst firing in many of their neurons, which increases the amplitude of signals in the neurons they contact. This difference in amplitude of signals serves as basis of selecting a sub-circuit for specialized processing (e.g., sustained attention) within the typically larger layer 6-based circuit. After examining the sustaining of oscillations in loop circuits and the processing of spikes in network circuits, we propose that cortical functioning can be globally viewed as two systems: a loop system and a network system. The loop system oscillations influence the network system’s timing and amplitude of pulse signals, both of which can select circuits that are momentarily dominant in cortical activity.

Is realistic neuronal modeling realistic?

Scientific models are abstractions that aim to explain natural phenomena. A successful model shows how a complex phenomenon arises from relatively simple principles while preserving major physical or biological rules and predicting novel experiments. A model should not be a facsimile of reality; it is an aid for understanding it. Contrary to this basic premise, with the 21st century has come a surge in computational efforts to model biological processes in great detail. Here we discuss the oxymoronic, realistic modeling of single neurons. This rapidly advancing field is driven by the discovery that some neurons don't merely sum their inputs and fire if the sum exceeds some threshold. Thus researchers have asked what are the computational abilities of single neurons and attempted to give answers using realistic models. We briefly review the state of the art of compartmental modeling highlighting recent progress and intrinsic flaws. We then attempt to address two fundamental questions. Practically, can we realistically model single neurons? Philosophically, should we realistically model single neurons? We use layer 5 neocortical pyramidal neurons as a test case to examine these issues. We subject three publically available models of layer 5 pyramidal neurons to three simple computational challenges. Based on their performance and a partial survey of published models, we conclude that current compartmental models are ad hoc, unrealistic models functioning poorly once they are stretched beyond the specific problems for which they were designed. We then attempt to plot possible paths for generating realistic single neuron models.

Identification and modelling of fast and slow Ih current components in vestibular ganglion neurons

Previous experimental data indicates the hyperpolarization‐activated cation (I_h) current, in the inner ear, consists of two components [different hyperpolarization‐activated cyclic nucleotide‐gated (HCN) subunits] which are impossible to pharmacologically isolate. To confirm the presence of these two components in vestibular ganglion neurons we have applied a parameter identification algorithm which is able to discriminate the parameters of the two components from experimental data. Using simulated data we have shown that this algorithm is able to identify the parameters of two populations of non‐inactivated ionic channels more accurately than a classical method. Moreover, the algorithm was demonstrated to be insensitive to the key parameter variations. We then applied this algorithm to I_h current recordings from mouse vestibular ganglion neurons. The algorithm revealed the presence of a high‐voltage‐activated slow component and a low‐voltage‐activated fast component. Finally, the electrophysiological significance of these two I_h components was tested individually in computational vestibular ganglion neuron models (sustained and transient), in the control case and in the presence of cAMP, an intracellular cyclic nucleotide that modulates HCN channel activity. The results suggest that, first, the fast and slow components modulate differently the action potential excitability and the excitatory postsynaptic potentials in both sustained and transient vestibular neurons and, second, the fast and slow components, in the control case, provide different information about characteristics of the stimulation and this information is significantly modified after modulation by cAMP.

A quantitative description of dendritic conductances and its application to dendritic excitation in layer 5 pyramidal neurons

Postsynaptic integration is a complex function of passive membrane properties and nonlinear activation of voltage-gated channels. Some cortical neurons express many voltage-gated channels, with each displaying heterogeneous dendritic conductance gradients. This complexity has hindered the construction of experimentally based mechanistic models of cortical neurons. Here we show that it is possible to overcome this obstacle. We recorded the membrane potential from the soma and apical dendrite of layer 5 (L5) pyramidal neurons of the rat somatosensory cortex. A combined experimental and numerical parameter peeling procedure was implemented to optimize a detailed ionic mechanism for the generation and propagation of dendritic spikes in neocortical L5 pyramidal neurons. In the optimized model, the density of voltage-gated Ca(2+) channels decreased linearly from the soma, and leveled at the distal apical dendrite. The density of the small-conductance Ca(2+)-activated channel decreased along the apical dendrite, whereas the density of the large-conductance Ca(2+)-gated K(+) channel was uniform throughout the apical dendrite. The model predicted an ionic mechanism for the generation of a dendritic spike, the interaction of this spike with the backpropagating action potential, the mechanism responsible for the ability of the proximal apical dendrite to control the coupling between the axon and the dendrite, and the generation of NMDA spikes in the distal apical tuft. Moreover, in addition to faithfully predicting many experimental results recorded from the apical dendrite of L5 pyramidal neurons, the model validates a new methodology for mechanistic modeling of neurons in the CNS.

Accelerating compartmental modeling on a graphical processing unit

Compartmental modeling is a widely used tool in neurophysiology but the detail and scope of such models is frequently limited by lack of computational resources. Here we implement compartmental modeling on low cost Graphical Processing Units (GPUs), which significantly increases simulation speed compared to NEURON. Testing two methods for solving the current diffusion equation system revealed which method is more useful for specific neuron morphologies. Regions of applicability were investigated using a range of simulations from a single membrane potential trace simulated in a simple fork morphology to multiple traces on multiple realistic cells. A runtime peak 150-fold faster than the CPU was achieved. This application can be used for statistical analysis and data fitting optimizations of compartmental models and may be used for simultaneously simulating large populations of neurons. Since GPUs are forging ahead and proving to be more cost-effective than CPUs, this may significantly decrease the cost of computation power and open new computational possibilities for laboratories with limited budgets.

Strength-duration relationship for intra- versus extracellular stimulation with microelectrodes

Chronaxie, a historically introduced excitability time parameter for electrical stimulation, has been assumed to be closely related to the time constant of the cell membrane. Therefore, it is perplexing that significantly larger chronaxies have been found for intracellular than for extracellular stimulation. Using compartmental model analysis, this controversy is explained on the basis that extracellular stimulation also generates hyperpolarized regions of the cell membrane hindering a steady excitation as seen in the intracellular case. The largest inside/outside chronaxie ratio for microelectrode stimulation is found in close vicinity of the cell. In the case of monophasic cathodic stimulation, the length of the primarily excited zone which is situated between the hyperpolarized regions increases with electrode-cell distance. For distant electrodes this results in an excitation process comparable to the temporal behavior of intracellular stimulation. Chronaxie also varies along the neural axis, being small for electrode positions at the nodes of Ranvier and axon initial segment and larger at the soma and dendrites. As spike initiation site can change for short and long pulses, in some cases strength-duration curves have a bimodal shape, and thus, they deviate from a classical monotonic curve as described by the formulas of Lapicque or Weiss.

Optimizing ion channel models using a parallel genetic algorithm on graphical processors

We have recently shown that we can semi-automatically constrain models of voltage-gated ion channels by combining a stochastic search algorithm with ionic currents measured using multiple voltage-clamp protocols. Although numerically successful, this approach is highly demanding computationally, with optimization on a high performance Linux cluster typically lasting several days. To solve this computational bottleneck we converted our optimization algorithm for work on a graphical processing unit (GPU) using NVIDIA's CUDA. Parallelizing the process on a Fermi graphic computing engine from NVIDIA increased the speed ∼180 times over an application running on an 80 node Linux cluster, considerably reducing simulation times. This application allows users to optimize models for ion channel kinetics on a single, inexpensive, desktop "super computer," greatly reducing the time and cost of building models relevant to neuronal physiology. We also demonstrate that the point of algorithm parallelization is crucial to its performance. We substantially reduced computing time by solving the ODEs (Ordinary Differential Equations) so as to massively reduce memory transfers to and from the GPU. This approach may be applied to speed up other data intensive applications requiring iterative solutions of ODEs.

A self-organizing state-space-model approach for parameter estimation in hodgkin-huxley-type models of single neurons

Traditional approaches to the problem of parameter estimation in biophysical models of neurons and neural networks usually adopt a global search algorithm (for example, an evolutionary algorithm), often in combination with a local search method (such as gradient descent) in order to minimize the value of a cost function, which measures the discrepancy between various features of the available experimental data and model output. In this study, we approach the problem of parameter estimation in conductance-based models of single neurons from a different perspective. By adopting a hidden-dynamical-systems formalism, we expressed parameter estimation as an inference problem in these systems, which can then be tackled using a range of well-established statistical inference methods. The particular method we used was Kitagawa's self-organizing state-space model, which was applied on a number of Hodgkin-Huxley-type models using simulated or actual electrophysiological data. We showed that the algorithm can be used to estimate a large number of parameters, including maximal conductances, reversal potentials, kinetics of ionic currents, measurement and intrinsic noise, based on low-dimensional experimental data and sufficiently informative priors in the form of pre-defined constraints imposed on model parameters. The algorithm remained operational even when very noisy experimental data were used. Importantly, by combining the self-organizing state-space model with an adaptive sampling algorithm akin to the Covariance Matrix Adaptation Evolution Strategy, we achieved a significant reduction in the variance of parameter estimates. The algorithm did not require the explicit formulation of a cost function and it was straightforward to apply on compartmental models and multiple data sets. Overall, the proposed methodology is particularly suitable for resolving high-dimensional inference problems based on noisy electrophysiological data and, therefore, a potentially useful tool in the construction of biophysical neuron models.

Parameter estimation is a problem of central importance and, perhaps, the most laborious task in biophysical modeling of neurons and neural networks. An emerging trend is to treat parameter estimation in this context as yet another statistical inference problem, which can be tackled using well-established methods from Computational Statistics. Inspired by these recent advances, we adopted a self-organizing state-space-model approach augmented with an adaptive sampling algorithm akin to the Covariance Matrix Adaptation Evolution Strategy in order to estimate a large number of parameters in a number of Hodgkin-Huxley-type models of single neurons. Parameter estimation was based on noisy electrophysiological data and involved the maximal conductances, reversal potentials, levels of noise and, unlike most mainstream work, the kinetics of ionic currents in the examined models. Our main conclusion was that parameters in complex, conductance-based neuron models can be inferred using the aforementioned methodology, if sufficiently informative priors regarding the unknown model parameters are available. Importantly, the use of an adaptive algorithm for sampling new parameter vectors significantly reduced the variance of parameter estimates. Flexibility and scalability are additional advantages of the proposed method, which is particularly suited to resolve high-dimensional inference problems.

Effective stimuli for constructing reliable neuron models

The rich dynamical nature of neurons poses major conceptual and technical challenges for unraveling their nonlinear membrane properties. Traditionally, various current waveforms have been injected at the soma to probe neuron dynamics, but the rationale for selecting specific stimuli has never been rigorously justified. The present experimental and theoretical study proposes a novel framework, inspired by learning theory, for objectively selecting the stimuli that best unravel the neuron's dynamics. The efficacy of stimuli is assessed in terms of their ability to constrain the parameter space of biophysically detailed conductance-based models that faithfully replicate the neuron's dynamics as attested by their ability to generalize well to the neuron's response to novel experimental stimuli. We used this framework to evaluate a variety of stimuli in different types of cortical neurons, ages and animals. Despite their simplicity, a set of stimuli consisting of step and ramp current pulses outperforms synaptic-like noisy stimuli in revealing the dynamics of these neurons. The general framework that we propose paves a new way for defining, evaluating and standardizing effective electrical probing of neurons and will thus lay the foundation for a much deeper understanding of the electrical nature of these highly sophisticated and non-linear devices and of the neuronal networks that they compose.

Theory of electric resonance in the neocortical apical dendrite

Pyramidal neurons of the neocortex display a wide range of synchronous EEG rhythms, which arise from electric activity along the apical dendrites of neocortical pyramidal neurons. Here we present a theoretical description of oscillation frequency profiles along apical dendrites which exhibit resonance frequencies in the range of 10 to 100 Hz. The apical dendrite is modeled as a leaky coaxial cable coated with a dielectric, in which a series of compartments act as coupled electric circuits that gradually narrow the resonance profile. The tuning of the peak frequency is assumed to be controlled by the average amplitude of voltage-gated outward currents, which in turn are regulated by the subthreshold noise in the thousands of synaptic spines that are continuously bombarded by local circuits. The results of simulations confirmed the ability of the model both to tune the peak frequency in the 10–100 Hz range and to gradually narrow the resonance profile. Considerable additional narrowing of the resonance profile is provided by repeated looping through the apical dendrite via the corticothalamocortical circuit, which reduced the width of each resonance curve (at half-maximum) to approximately 1 Hz. Synaptic noise in the neural circuit is discussed in relation to the ways it can influence the narrowing process.

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.

Mechanisms of magnetic stimulation of central nervous system neurons

Transcranial magnetic stimulation (TMS) is a stimulation method in which a magnetic coil generates a magnetic field in an area of interest in the brain. This magnetic field induces an electric field that modulates neuronal activity. The spatial distribution of the induced electric field is determined by the geometry and location of the coil relative to the brain. Although TMS has been used for several decades, the biophysical basis underlying the stimulation of neurons in the central nervous system (CNS) is still unknown. To address this problem we developed a numerical scheme enabling us to combine realistic magnetic stimulation (MS) with compartmental modeling of neurons with arbitrary morphology. The induced electric field for each location in space was combined with standard compartmental modeling software to calculate the membrane current generated by the electromagnetic field for each segment of the neuron. In agreement with previous studies, the simulations suggested that peripheral axons were excited by the spatial gradients of the induced electric field. In both peripheral and central neurons, MS amplitude required for action potential generation was inversely proportional to the square of the diameter of the stimulated compartment. Due to the importance of the fiber's diameter, magnetic stimulation of CNS neurons depolarized the soma followed by initiation of an action potential in the initial segment of the axon. Passive dendrites affect this process primarily as current sinks, not sources. The simulations predict that neurons with low current threshold are more susceptible to magnetic stimulation. Moreover, they suggest that MS does not directly trigger dendritic regenerative mechanisms. These insights into the mechanism of MS may be relevant for the design of multi-intensity TMS protocols, may facilitate the construction of magnetic stimulators, and may aid the interpretation of results of TMS of the CNS.

Transcranial magnetic stimulation (TMS) is a widely applied tool for probing cognitive function in humans and is one of the best tools for clinical treatments and interfering with cognitive tasks. Surprisingly, while TMS has been commercially available for decades, the cellular mechanisms underlying magnetic stimulation remain unclear. Here we investigate these mechanisms using compartmental modeling. We generated a numerical scheme allowing simulation of the physiological response to magnetic stimulation of neurons with arbitrary morphologies and active properties. Computational experiments using this scheme suggested that TMS affects neurons in the central nervous system (CNS) primarily by somatic stimulation. Since magnetic stimulation appears to cause somatic depolarization, its effects are highly correlated with the neuron's current threshold. Our simulations therefore predict that subpopulations of CNS neurons with different firing thresholds will respond differently to magnetic stimulation. For example, low-intensity TMS may be used to stimulate low-threshold cortical inhibitory interneurons. At higher intensities we predict that both inhibitory and excitatory neurons are activated. These predictions may be tested at the cellular level and may impact cognitive experiments in humans. Furthermore, our simulations may be used to design TMS coils, devices and protocols.

Kinetic modeling of Nav1.7 provides insight into erythromelalgia-associated F1449V mutation

Gain-of-function mutations of the voltage-gated sodium channel (VGSC) Na(v)1.7 have been linked to human pain disorders. The mutation F1449V, located at the intracellular end of transmembrane helix S6 of domain III, induces the inherited pain syndrome erythromelalgia. A kinetic model of wild-type (WT) and F1449V Na(v)1.7 may provide a basis for predicting putative intraprotein interactions. We semiautomatically constrained a Markov model using stochastic search algorithms and whole cell patch-clamp recordings from human embryonic kidney cells transfected with Na(v)1.7 and its F1449V mutation. The best models obtained simulated known differences in action potential thresholds and firing patterns in spinal sensory neurons expressing WT and F1449V. The most suitable Markov model consisted of three closed, one open, and two inactivated states. The model predicted that the F1449V mutation shifts occupancy of the closed states closer to the open state, making it easier for the channel pore to open. It also predicted that F1449V's second inactivated state is more than four times more likely to be occupied than the equivalent state in WT at hyperpolarized potentials, although the mutation still lowered the firing threshold of action potentials. The differences between WT and F1449V were not limited to a single transition. Thus a point mutation in position F1449, while phenotypically most probably affecting the activation gate, may also modify channel functions mediated by structures in more distant areas of the channel protein.

Which elements of the mammalian central nervous system are excited by low current stimulation with microelectrodes?

Low current cortex stimulation produces a sparse and distributed set of activated cells often with distances of several hundred micrometers between cell bodies and the microelectrode. A modeling study based on recently measured densities of high threshold sodium channels Nav1.2 in dendrites and soma and low threshold sodium channels Nav1.6 in the axon shall identify spike initiation sites including a discussion on dendritic spikes. Varying excitability along the neural axis has been observed while studying different electrode positions and configurations. Although the axon initial segment (AIS) and nodes of Ranvier are most excitable, many thin axons and dendrites which are likely to be close to the electrode in the densely packed cortical regions are also proper candidates for spike initiation sites. Cathodic threshold ratio for thin axons and dendrites is about 1:3, whereas 0.2 mum diameter axons passing the electrode tip in 10 mum distance can be activated by 100 mus pulses with 2.6 muA. Direct cathodic excitation of dendrites requires a minimum electrode-fiber distance, which increases with dendrite diameter. Therefore thin dendrites can profit from the stronger electrical field close to the electrode but low current stimulation cannot activate large diameter dendrites, contrary to the inverse recruitment order known from peripheral nerve stimulation. When local depolarization fails to generate a dendritic spike, stimulation is possible via intracellular current flow that initiates an action potential, for example 200 mum distant in the low threshold AIS or in certain cases at the distal dendrite ending. Beside these exceptions, spike initiation site for cathodic low current stimulation appears rather close to the electrode.