Journal of Mathematical Neuroscience最新文献

The modeling of neural fields in the visual cortex involves geometrical structures which describe in mathematical formalism the functional architecture of this cortical area. The case of contour detection and orientation tuning has been extensively studied and has become a paradigm for the mathematical analysis of image processing by the brain. Ten years ago an attempt was made to extend these models by replacing orientation (an angle) with a second-order tensor built from the gradient of the image intensity, and it was named the structure tensor. This assumption does not follow from biological observations (experimental evidence is still lacking) but from the idea that the effectiveness of texture processing with the structure tensor in computer vision may well be exploited by the brain itself. The drawback is that in this case the geometry is not Euclidean but hyperbolic instead, which complicates the analysis substantially. The purpose of this review is to present the methodology that was developed in a series of papers to investigate this quite unusual problem, specifically from the point of view of tuning and pattern formation. These methods, which rely on bifurcation theory with symmetry in the hyperbolic context, might be of interest for the modeling of other features such as color vision or other brain functions.

Neurons are biological cells with uniquely complex dendritic morphologies that are not present in other cell types. Electrical signals in a neuron with branching dendrites can be studied by cable theory which provides a general mathematical modelling framework of spatio-temporal voltage dynamics. Typically such models need to be solved numerically unless the cell membrane is modelled either by passive or quasi-active dynamics, in which cases analytical solutions can be reduced to calculation of the Green's function describing the fundamental input-output relationship in a given morphology. Such analytically tractable models often assume individual dendritic segments to be cylinders. However, it is known that dendritic segments in many types of neurons taper, i.e. their radii decline from proximal to distal ends. Here we consider a generalised form of cable theory which takes into account both branching and tapering structures of dendritic trees. We demonstrate that analytical solutions can be found in compact algebraic forms in an arbitrary branching neuron with a class of tapering dendrites studied earlier in the context of single neuronal cables by Poznanski (Bull. Math. Biol. 53(3):457-467, 1991). We apply this extended framework to a number of simplified neuronal models and contrast their output dynamics in the presence of tapering versus cylindrical segments.

The canonical computational model for the cognitive process underlying two-alternative forced-choice decision making is the so-called drift-diffusion model (DDM). In this model, a decision variable keeps track of the integrated difference in sensory evidence for two competing alternatives. Here I extend the notion of a drift-diffusion process to multiple alternatives. The competition between n alternatives takes place in a linear subspace of [Formula: see text] dimensions; that is, there are [Formula: see text] decision variables, which are coupled through correlated noise sources. I derive the multiple-alternative DDM starting from a system of coupled, linear firing rate equations. I also show that a Bayesian sequential probability ratio test for multiple alternatives is, in fact, equivalent to these same linear DDMs, but with time-varying thresholds. If the original neuronal system is nonlinear, one can once again derive a model describing a lower-dimensional diffusion process. The dynamics of the nonlinear DDM can be recast as the motion of a particle on a potential, the general form of which is given analytically for an arbitrary number of alternatives.

Background: Trigeminal neuralgia (TN) is a severe neuropathic pain, which has an electric shock-like characteristic. There are some common treatments for this pain such as medicine, microvascular decompression or radio frequency. In this regard, transcranial direct current stimulation (tDCS) is another therapeutic method to reduce pain, which has been recently attracting the therapists' attention. The positive effect of tDCS on TN was shown in many previous studies. However, the mechanism of the tDCS effect has remained unclear.

Objective: This study aims to model the neuronal behavior of the main known regions of the brain participating in TN pathways to study the effect of transcranial direct current stimulation.

Method: The proposed model consists of several blocks: (1) trigeminal nerve, (2) trigeminal ganglion, (3) PAG (periaqueductal gray in the brainstem), (4) thalamus, (5) motor cortex (M1) and (6) somatosensory cortex (S1). Each of these components is represented by a modified Hodgkin-Huxley (HH) model. The modification of the HH model was done based on some neurological facts of pain sodium channels. The input of the model involves any stimuli to the 'trigeminal nerve,' which cause the pain, and the output is the activity of the somatosensory cortex. An external current, which is considered as an electrical current, was applied to the motor cortex block of the model.

Result: The results showed that by decreasing the conductivity of the slow sodium channels (pain channels) and applying tDCS over the M1, the activity of the somatosensory cortex would be reduced. This reduction can cause pain relief.

Conclusion: The proposed model provided some possible suggestions about the relationship between the effects of tDCS and associated components in TN, and also the relationship between the pain measurement index, somatosensory cortex activity, and the strength of tDCS.

To establish and exploit novel biomarkers of demyelinating diseases requires a mechanistic understanding of axonal propagation. Here, we present a novel computational framework called the stochastic spike-diffuse-spike (SSDS) model for assessing the effects of demyelination on axonal transmission. It models transmission through nodal and internodal compartments with two types of operations: a stochastic integrate-and-fire operation captures nodal excitability and a linear filtering operation describes internodal propagation. The effects of demyelinated segments on the probability of transmission, transmission delay and spike time jitter are explored. We argue that demyelination-induced impedance mismatch prevents propagation mostly when the action potential leaves a demyelinated region, not when it enters a demyelinated region. In addition, we model sodium channel remodeling as a homeostatic control of nodal excitability. We find that the effects of mild demyelination on transmission probability and delay can be largely counterbalanced by an increase in excitability at the nodes surrounding the demyelination. The spike timing jitter, however, reflects the level of demyelination whether excitability is fixed or is allowed to change in compensation. This jitter can accumulate over long axons and leads to a broadening of the compound action potential, linking microscopic defects to a mesoscopic observable. Our findings articulate why action potential jitter and compound action potential dispersion can serve as potential markers of weak and sporadic demyelination.

Understanding nervous system function requires careful study of transient (non-equilibrium) neural response to rapidly changing, noisy input from the outside world. Such neural response results from dynamic interactions among multiple, heterogeneous brain regions. Realistic modeling of these large networks requires enormous computational resources, especially when high-dimensional parameter spaces are considered. By assuming quasi-steady-state activity, one can neglect the complex temporal dynamics; however, in many cases the quasi-steady-state assumption fails. Here, we develop a new reduction method for a general heterogeneous firing-rate model receiving background correlated noisy inputs that accurately handles highly non-equilibrium statistics and interactions of heterogeneous cells. Our method involves solving an efficient set of nonlinear ODEs, rather than time-consuming Monte Carlo simulations or high-dimensional PDEs, and it captures the entire set of first and second order statistics while allowing significant heterogeneity in all model parameters.

The ongoing acquisition of large and multifaceted data sets in neuroscience requires new mathematical tools for quantitatively grounding these experimental findings. Since 2015, the International Conference on Mathematical Neuroscience (ICMNS) has provided a forum for researchers to discuss current mathematical innovations emerging in neuroscience. This special issue assembles current research and tutorials that were presented at the 2017 ICMNS held in Boulder, Colorado from May 30 to June 2. Topics discussed at the meeting include correlation analysis of network activity, information theory for plastic synapses, combinatorics for attractor neural networks, and novel data assimilation methods for neuroscience-all of which are represented in this special issue.