Search

The sense of touch is represented by neural activity patterns evoked by mechanosensory input forces. The rodent whisker system is exceptional for studying the neurophysiology of touch in part because these forces can be precisely computed from video of whisker deformation. We evaluate the accuracy of a standard model of whisker bending, which assumes quasi-static dynamics and a linearly tapered conical profile, using controlled whisker deflections. We find significant discrepancies between model and experiment: real whiskers bend more than predicted upon contact at locations in the middle of the whisker and less at distal locations. Thus whiskers behave as if their stiffness near the base and near the tip is larger than expected for a homogeneous cone. We assess whether contact direction, friction, inhomogeneous elasticity, whisker orientation, or nonconical shape could explain these deviations. We show that a thin-middle taper of mouse whisker shape accounts for the majority of this behavior. This taper is conserved across rows and columns of the whisker array. The taper has a large effect on the touch-evoked forces and the ease with which whiskers slip past objects, which are key drivers of neural activity in tactile object localization and identification. This holds for orientations with intrinsic whisker curvature pointed toward, away from, or down from objects, validating two-dimensional models of simple whisker-object interactions. The precision of computational models relating sensory input forces to neural activity patterns can be quantitatively enhanced by taking thin-middle taper into account with a simple corrective function that we provide., (Copyright © 2016 the American Physiological Society.)

Complex olfactory-discrimination (OD) learning results in a series of intrinsic and excitatory synaptic modifications in piriform cortex pyramidal neurons that enhance the circuit excitability. Such overexcitation must be balanced to prevent runway activity while maintaining the efficient ability to store memories. We showed previously that OD learning is accompanied by enhancement of the GABAA-mediated inhibition. Here we show that GABAB-mediated inhibition is also enhanced after learning and study the mechanism underlying such enhancement and explore its functional role. We show that presynaptic, GABAB-mediated synaptic inhibition is enhanced after learning. In contrast, the population-average postsynaptic GABAB-mediated synaptic inhibition is unchanged, but its standard deviation is enhanced. Learning-induced reduction in paired pulse facilitation in the glutamatergic synapses interconnecting pyramidal neurons was abolished by application of the GABAB antagonist CGP55845 but not by blocking G protein-gated inwardly rectifying potassium channels only, indicating enhanced suppression of excitatory synaptic release via presynaptic GABAB-receptor activation. In addition, the correlation between the strengths of the early (GABAA-mediated) and late (GABAB-mediated) synaptic inhibition was much stronger for each particular neuron after learning. Consequently, GABAB-mediated inhibition was also more efficient in controlling epileptic-like activity induced by blocking GABAA receptors. We suggest that complex OD learning is accompanied by enhancement of the GABAB-mediated inhibition that enables the cortical network to store memories, while preventing uncontrolled activity., (Copyright © 2014 the American Physiological Society.)

Vibrissa motoneurons (vMNs) generate rhythmic firing that controls whisker movements, even without cortical, cerebellar, or sensory inputs. vMNs receive serotonergic modulation from brain stem areas, which mainly increases their persistent sodium conductance (g(NaP)) and, possibly, phasic input from a putative central pattern generator (CPG). In response to serotonergic modulation or just-suprathreshold current steps, vMNs fire at low rates, below the firing frequency of exploratory whisking. In response to periodic inputs, vMNs exhibit nonlinear suprathreshold resonance in frequency ranges of exploratory whisking. To determine how firing patterns of vMNs are determined by their 1) intrinsic ionic conductances and 2) responses to periodic input from a putative CPG and to serotonergic modulation, we construct and analyze a single-compartment, conductance-based model of vMNs. Low firing rates are supported in extended regimes by adaptation currents and the minimal firing rate decreases with g(NaP) and increases with M-potassium and h-cation conductances. Suprathreshold resonance results from the locking properties of vMN firing to stimuli and from reduction of firing rates at low frequencies by slow M and afterhyperpolarization potassium conductances. h conductance only slightly affects the suprathreshold resonance. When a vMN is subjected to a small periodic CPG input, serotonergically induced g(NaP) elevation may transfer the system from quiescence to a firing state that is highly locked to the CPG input. Thus we conclude that for vMNs, the CPG controls firing frequency and phase and enables bursting, whereas serotonergic modulation controls transitions from quiescence to firing unless the CPG input is sufficiently strong.

The intrinsic firing modes of adult CA1 pyramidal cells vary along a continuum of "burstiness" from regular firing to rhythmic bursting, depending on the ionic composition of the extracellular milieu. Burstiness is low in neurons exposed to a normal extracellular Ca(2+) concentration ([Ca(2+)](o)), but is markedly enhanced by lowering [Ca(2+)](o), although not by blocking Ca(2+) and Ca(2+)-activated K(+) currents. We show, using intracellular recordings, that burstiness in low [Ca(2+)](o) persists even after truncating the apical dendrites, suggesting that bursts are generated by an interplay of membrane currents at or near the soma. To study the mechanisms of bursting, we have constructed a conductance-based, one-compartment model of CA1 pyramidal neurons. In this neuron model, reduced [Ca(2+)](o) is simulated by negatively shifting the activation curve of the persistent Na(+) current (I(NaP)) as indicated by recent experimental results. The neuron model accounts, with different parameter sets, for the diversity of firing patterns observed experimentally in both zero and normal [Ca(2+)](o). Increasing I(NaP) in the neuron model induces bursting and increases the number of spikes within a burst but is neither necessary nor sufficient for bursting. We show, using fast-slow analysis and bifurcation theory, that the M-type K(+) current (I(M)) allows bursting by shifting neuronal behavior between a silent and a tonically active state provided the kinetics of the spike generating currents are sufficiently, although not extremely, fast. We suggest that bursting in CA1 pyramidal cells can be explained by a single compartment "square bursting" mechanism with one slow variable, the activation of I(M).

A temporal sensory code occurs in posterior medial (POm) thalamus of the rat vibrissa system, where the latency for the spike rate to peak is observed to increase with increasing frequency of stimulation between 2 and 11 Hz. In contrast, the latency of the spike rate in the ventroposterior medial (VPm) thalamus is constant in this frequency range. We consider the hypothesis that two factors are essential for latency coding in the POm. The first is GABAB-mediated feedback inhibition from the reticular thalamic (Rt) nucleus, which provides delayed and prolonged input to thalamic structures. The second is sensory input that leads to an accelerating spike rate in brain stem nuclei. Essential aspects of the experimental observations are replicated by the analytical solution of a rate-based model with a minimal architecture that includes only the POm and Rt nuclei, i.e., an increase in stimulus frequency will increase the level of inhibitory output from Rt thalamus and lead to a longer latency in the activation of POm thalamus. This architecture, however, admits period-doubling at high levels of GABAB-mediated conductance. A full architecture that incorporates the VPm nucleus suppresses period-doubling. A clear match between the experimentally measured spike rates and the numerically calculated rates for the full model occurs when VPm thalamus receives stronger brain stem input and weaker GABAB-mediated inhibition than POm thalamus. Our analysis leads to the prediction that the latency code will disappear if GABAB-mediated transmission is blocked in POm thalamus or if the onset of sensory input is too abrupt. We suggest that GABAB-mediated inhibition is a substrate of temporal coding in normal brain function.

We explore the mechanism of synchronized bursting activity with frequency of approximately 10 Hz that appears in cortical tissues at low extracellular magnesium concentration [Mg2+]o. We hypothesize that this activity is persistent, namely coexists with the quiescent state and depends on slow N-methyl-D-aspartate (NMDA) conductances. To explore this hypothesis, we construct and investigate a conductance-based model of excitatory cortical networks. Population bursting activity can persist for physiological values of the NMDA decay time constant (approximately 100 ms). Neurons are synchronized at the time scale of bursts but not of single spikes. A reduced model of a cell coupled to itself can encompass most of this highly synchronized network behavior and is analyzed using the fast-slow method. Synchronized bursts appear for intermediate values of the NMDA conductance g(NMDA) if NMDA conductances are not too fast. Regular spiking activity appears for larger g(NMDA). If the single cell is a conditional burster, persistent synchronized bursts become more robust. Weakly synchronized states appear for zero AMPA conductance g(AMPA). Enhancing g(AMPA) increases both synchrony and the number of spikes within bursts and decreases the bursting frequency. Too strong g(AMPA), however, prevents the activity because it enhances neuronal intrinsic adaptation. When [Mg2+]o is increased, higher g(NMDA) values are needed to maintain bursting activity. Bursting frequency decreases with [Mg2+]o, and the network is silent with physiological [Mg2+]o. Inhibition weakly decreases the bursting frequency if inhibitory cells receive enough NMDA-mediated excitation. This study explains the importance of conditional bursters in layer V in supporting epileptiform activity at low [Mg2+]o.

Stereotypic paroxysmal discharges that propagate in neocortical tissues after electrical stimulations are used as a probe for studying cortical circuitry. I use modeling to investigate the effects of sparse connectivity, heterogeneity of intrinsic neuronal properties, and synaptic noise on synchronization of evoked propagating neuronal discharges in a network of excitatory, regular spiking neurons with spatially decaying connectivity. The global coherence of the traveling discharge is characterized by the correlation function between spike trains of neurons, averaged over all the pairs of neurons in the system at the same distance. Local coherence of two neurons is characterized by their correlation function averaged over many trials or, for persistent activity, over a long time interval. Spike synchronization between neurons emerges as a result of the transient activity; if activity is persistent, there is no synchrony, and cross-correlation functions are flat. During discharge propagation, system-average cross-correlation between neurons does not depend on their mutual distance except for a time shift. Spike synchronization occurs only when the average number of synapses M a cell receives is large enough. As M increases, there is a cross-over from an asynchronized to a synchronized discharge. Synaptic depression appears to help synchrony; it reduces the M value at the cross-over. The strengths of alpha-amino-3-hydroxy-5-methyl-4-isoxazolepropionic acid (AMPA) and N-methyl-D-aspartate (NMDA) conductances affect synchrony only weakly. Spike synchronization is robust even with large levels of heterogeneity. Synaptic noise reduces synchrony, but strong synchrony is observed at a noise level that cannot evoke spontaneous discharges. System-average spike synchronization is determined by the levels of sparseness, heterogeneity, and noise, whereas trial-average spike synchronization is determined only by the noise level. Therefore, I predict that experiments will reveal local, two-cell spike synchrony, but not global synchrony.

We studied the propagation of paroxysmal discharges in disinhibited neocortical slices by developing and analyzing a model of excitatory regular-spiking neocortical cells with spatially decaying synaptic efficacies and by field potential recording in rat slices. Evoked discharges may propagate both in the model and in the experiment. The model discharge propagates as a traveling pulse with constant velocity and shape. The discharge shape is determined by an interplay between the synaptic driving force and the neuron's intrinsic currents, in particular the slow potassium current. In the model, N-methyl-D-aspartate (NMDA) conductance contributes much less to the discharge velocity than amino-3-hydroxy-5-methyl-4-isoxazolepropionic acid (AMPA) conductance. Blocking NMDA receptors experimentally with 2-amino-5-phosphonovaleric acid (APV) has no significant effect on the discharge velocity. In both model and experiments, propagation occurs for AMPA synaptic coupling gAMPA above a certain threshold, at which the velocity is finite (non-zero). The discharge velocity grows linearly with the gAMPA for gAMPA much above the threshold. In the experiments, blocking AMPA receptors gradually by increasing concentrations of 6-cyano-7-nitroquinoxaline-2,3-dione (CNQX) in the perfusing solution results in a gradual reduction of the discharge velocity until propagation stops altogether, thus confirming the model prediction. When discharges are terminated in the model by the slow potassium current, a network with the same parameter set may display discharges with several forms, which have different velocities and numbers of spikes; initial conditions select the exhibited pattern. When the discharge is also terminated by strong synaptic depression, there is only one discharge form for a particular parameter set; the velocity grows continuously with increased synaptic conductances. No indication for more than one discharge velocity was observed experimentally. If the AMPA decay rate increases while the maximal excitatory postsynaptic conductance (EPSC) a cell receives is kept fixed, the velocity increases by approximately 20% until it reaches a saturated value. Therefore the discharge velocity is determined mainly by the cells' integration time of input EPSCs. We conclude, on the basis of both the experiments and the model, that the total amount of excitatory conductance a typical cell receives in a control slice exhibiting paroxysmal discharges is only approximately 5 times larger than the excitatory conductance needed for raising the potential of a resting cell above its action potential threshold.