Action-driven remapping of hippocampal neuronal populations in jumping rats.

The current dominant view of the hippocampus is that it is a navigation "device" guided by environmental inputs. Yet, a critical aspect of navigation is a sequence of planned, coordinated actions. We examined the role of action in the neuronal organization of the hippocampus by training rats to jump a gap on a linear track. Recording local field potentials and ensembles of single units in the hippocampus, we found that jumping produced a stereotypic behavior associated with consistent electrophysiological patterns, including phase reset of theta oscillations, predictable global firing-rate changes, and population vector shifts of hippocampal neurons. A subset of neurons ("jump cells") were systematically affected by the gap but only in one direction of travel. Novel place fields emerged and others were either boosted or attenuated by jumping, yet the theta spike phase versus animal position relationship remained unaltered. Thus, jumping involves an action plan for the animal to traverse the same route as without jumping, which is faithfully tracked by hippocampal neuronal activity.

Classification of bursting patterns: A tale of two ducks.

Bursting is one of the fundamental rhythms that excitable cells can generate either in response to incoming stimuli or intrinsically. It has been a topic of intense research in computational biology for several decades. The classification of bursting oscillations in excitable systems has been the subject of active research since the early 1980s and is still ongoing. As a by-product, it establishes analytical and numerical foundations for studying complex temporal behaviors in multiple timescale models of cellular activity. In this review, we first present the seminal works of Rinzel and Izhikevich in classifying bursting patterns of excitable systems. We recall a complementary mathematical classification approach by Bertram and colleagues, and then by Golubitsky and colleagues, which, together with the Rinzel-Izhikevich proposals, provide the state-of-the-art foundations to these classifications. Beyond classical approaches, we review a recent bursting example that falls outside the previous classification systems. Generalizing this example leads us to propose an extended classification, which requires the analysis of both fast and slow subsystems of an underlying slow-fast model and allows the dissection of a larger class of bursters. Namely, we provide a general framework for bursting systems with both subthreshold and superthreshold oscillations. A new class of bursters with at least 2 slow variables is then added, which we denote folded-node bursters, to convey the idea that the bursts are initiated or annihilated via a folded-node singularity. Key to this mechanism are so-called canard or duck orbits, organizing the underpinning excitability structure. We describe the 2 main families of folded-node bursters, depending upon the phase (active/spiking or silent/nonspiking) of the bursting cycle during which folded-node dynamics occurs. We classify both families and give examples of minimal systems displaying these novel bursting patterns. Finally, we provide a biophysical example by reinterpreting a generic conductance-based episodic burster as a folded-node burster, showing that the associated framework can explain its subthreshold oscillations over a larger parameter region than the fast subsystem approach.

Bistability at the onset of neuronal oscillations.

The Hodgkin-Huxley (HH) model and squid axon (bathed in reduced Ca2+) fire repetitively for steady current injection. Moreover, for a current-range just suprathreshold, repetitive firing coexists with a stable steady state. Neuronal excitability, as such, shows bistability and hysteresis providing the opportunity for the system to perform as switchable between firing and non-firing states with transient input and providing the backbone as a dynamical mechanism for bursting oscillations. Some conditions for bistability can be derived by intricate analysis (bifurcation theory) and characterized by simulation, but conditions for emergence and robustness of such bistability do not typically follow from intuition. Here, we demonstrate with a semi-quantitative two-variable, V-w, reduction of the HH model features that promote/reduce bistability. Visualization of flow and trajectories in the V-w phase plane provides an intuitive grasp for bistability. The geometry of action potential recovery involves a late phase during which the dynamic negative feedback of [Formula: see text] inactivation and [Formula: see text] activation over/undershoot, respectively, their resting values, thereby leading to hyperexcitabilty and an intrinsically generated opportunity to by-pass the spiral-like stable rest state and initiate the next spike upstroke. We illustrate control of bistability and dependence of the degree of hysteresis on recovery timescales and gating properties. Our dynamical dissection reveals the strongly attracting depolarized phase of the spike, enabling approximations like the resetting feature of adapting integrate-and-fire models. We extend our insights and show that the Morris-Lecar model can also exhibit robust bistability.

A biophysical counting mechanism for keeping time.

The ability to estimate and produce appropriately timed responses is central to many behaviors including speaking, dancing, and playing a musical instrument. A classical framework for estimating or producing a time interval is the pacemaker-accumulator model in which pulses of a pacemaker are counted and compared to a stored representation. However, the neural mechanisms for how these pulses are counted remain an open question. The presence of noise and stochasticity further complicates the picture. We present a biophysical model of how to keep count of a pacemaker in the presence of various forms of stochasticity using a system of bistable Wilson-Cowan units asymmetrically connected in a one-dimensional array; all units receive the same input pulses from a central clock but only one unit is active at any point in time. With each pulse from the clock, the position of the activated unit changes thereby encoding the total number of pulses emitted by the clock. This neural architecture maps the counting problem into the spatial domain, which in turn translates count to a time estimate. We further extend the model to a hierarchical structure to be able to robustly achieve higher counts.

Stimulus-Driven and Spontaneous Dynamics in Excitatory-Inhibitory Recurrent Neural Networks for Sequence Representation.

Recurrent neural networks (RNNs) have been widely used to model sequential neural dynamics ("neural sequences") of cortical circuits in cognitive and motor tasks. Efforts to incorporate biological constraints and Dale's principle will help elucidate the neural representations and mechanisms of underlying circuits. We trained an excitatory-inhibitory RNN to learn neural sequences in a supervised manner and studied the representations and dynamic attractors of the trained network. The trained RNN was robust to trigger the sequence in response to various input signals and interpolated a time-warped input for sequence representation. Interestingly, a learned sequence can repeat periodically when the RNN evolved beyond the duration of a single sequence. The eigenspectrum of the learned recurrent connectivity matrix with growing or damping modes, together with the RNN's nonlinearity, were adequate to generate a limit cycle attractor. We further examined the stability of dynamic attractors while training the RNN to learn two sequences. Together, our results provide a general framework for understanding neural sequence representation in the excitatory-inhibitory RNN.

Buildup and bistability in auditory streaming as an evidence accumulation process with saturation.

A repeating triplet-sequence ABA- of non-overlapping brief tones, A and B, is a valued paradigm for studying auditory stream formation and the cocktail party problem. The stimulus is "heard" either as a galloping pattern (integration) or as two interleaved streams (segregation); the initial percept is typically integration then followed by spontaneous alternations between segregation and integration, each being dominant for a few seconds. The probability of segregation grows over seconds, from near-zero to a steady value, defining the buildup function, BUF. Its stationary level increases with the difference in tone frequencies, DF, and the BUF rises faster. Percept durations have DF-dependent means and are gamma-like distributed. Behavioral and computational studies usually characterize triplet streaming either during alternations or during buildup. Here, our experimental design and modeling encompass both. We propose a pseudo-neuromechanistic model that incorporates spiking activity in primary auditory cortex, A1, as input and resolves perception along two network-layers downstream of A1. Our model is straightforward and intuitive. It describes the noisy accumulation of evidence against the current percept which generates switches when reaching a threshold. Accumulation can saturate either above or below threshold; if below, the switching dynamics resemble noise-induced transitions from an attractor state. Our model accounts quantitatively for three key features of data: the BUFs, mean durations, and normalized dominance duration distributions, at various DF values. It describes perceptual alternations without competition per se, and underscores that treating triplets in the sequence independently and averaging across trials, as implemented in earlier widely cited studies, is inadequate.

A neuromechanistic model for rhythmic beat generation.

When listening to music, humans can easily identify and move to the beat. Numerous experimental studies have identified brain regions that may be involved with beat perception and representation. Several theoretical and algorithmic approaches have been proposed to account for this ability. Related to, but different from the issue of how we perceive a beat, is the question of how we learn to generate and hold a beat. In this paper, we introduce a neuronal framework for a beat generator that is capable of learning isochronous rhythms over a range of frequencies that are relevant to music and speech. Our approach combines ideas from error-correction and entrainment models to investigate the dynamics of how a biophysically-based neuronal network model synchronizes its period and phase to match that of an external stimulus. The model makes novel use of on-going faster gamma rhythms to form a set of discrete clocks that provide estimates, but not exact information, of how well the beat generator spike times match those of a stimulus sequence. The beat generator is endowed with plasticity allowing it to quickly learn and thereby adjust its spike times to achieve synchronization. Our model makes generalizable predictions about the existence of asymmetries in the synchronization process, as well as specific predictions about resynchronization times after changes in stimulus tempo or phase. Analysis of the model demonstrates that accurate rhythmic time keeping can be achieved over a range of frequencies relevant to music, in a manner that is robust to changes in parameters and to the presence of noise.

Soma-axon coupling configurations that enhance neuronal coincidence detection.

Coincidence detector neurons transmit timing information by responding preferentially to concurrent synaptic inputs. Principal cells of the medial superior olive (MSO) in the mammalian auditory brainstem are superb coincidence detectors. They encode sound source location with high temporal precision, distinguishing submillisecond timing differences among inputs. We investigate computationally how dynamic coupling between the input region (soma and dendrite) and the spike-generating output region (axon and axon initial segment) can enhance coincidence detection in MSO neurons. To do this, we formulate a two-compartment neuron model and characterize extensively coincidence detection sensitivity throughout a parameter space of coupling configurations. We focus on the interaction between coupling configuration and two currents that provide dynamic, voltage-gated, negative feedback in subthreshold voltage range: sodium current with rapid inactivation and low-threshold potassium current, IKLT. These currents reduce synaptic summation and can prevent spike generation unless inputs arrive with near simultaneity. We show that strong soma-to-axon coupling promotes the negative feedback effects of sodium inactivation and is, therefore, advantageous for coincidence detection. Furthermore, the feedforward combination of strong soma-to-axon coupling and weak axon-to-soma coupling enables spikes to be generated efficiently (few sodium channels needed) and with rapid recovery that enhances high-frequency coincidence detection. These observations detail the functional benefit of the strongly feedforward configuration that has been observed in physiological studies of MSO neurons. We find that IKLT further enhances coincidence detection sensitivity, but with effects that depend on coupling configuration. For instance, in models with weak soma-to-axon and weak axon-to-soma coupling, IKLT in the axon enhances coincidence detection more effectively than IKLT in the soma. By using a minimal model of soma-to-axon coupling, we connect structure, dynamics, and computation. Although we consider the particular case of MSO coincidence detectors, our method for creating and exploring a parameter space of two-compartment models can be applied to other neurons.

Canard-induced complex oscillations in an excitatory network.

In Neuroscience, mathematical modelling involving multiple spatial and temporal scales can unveil complex oscillatory activity such as excitable responses to an input current, subthreshold oscillations, spiking or bursting. While the number of slow and fast variables and the geometry of the system determine the type of the complex oscillations, canard structures define boundaries between them. In this study, we use geometric singular perturbation theory to identify and characterise boundaries between different dynamical regimes in multiple-timescale firing rate models of the developing spinal cord. These rate models are either three or four dimensional with state variables chosen within an overall group of two slow and two fast variables. The fast subsystem corresponds to a recurrent excitatory network with fast activity-dependent synaptic depression, and the slow variables represent the cell firing threshold and slow activity-dependent synaptic depression, respectively. We start by demonstrating canard-induced bursting and mixed-mode oscillations in two different three-dimensional rate models. Then, in the full four-dimensional model we show that a canard-mediated slow passage creates dynamics that combine these complex oscillations and give rise to mixed-mode bursting oscillations (MMBOs). We unveil complicated isolas along which MMBOs exist in parameter space. The profile of solutions along each isola undergoes canard-mediated transitions between the sub-threshold regime and the bursting regime; these explosive transitions change the number of oscillations in each regime. Finally, we relate the MMBO dynamics to experimental recordings and discuss their effects on the silent phases of bursting patterns as well as their potential role in creating subthreshold fluctuations that are often interpreted as noise. The mathematical framework used in this paper is relevant for modelling multiple timescale dynamics in excitable systems.

Attention model of binocular rivalry.

When the corresponding retinal locations in the two eyes are presented with incompatible images, a stable percept gives way to perceptual alternations in which the two images compete for perceptual dominance. As perceptual experience evolves dynamically under constant external inputs, binocular rivalry has been used for studying intrinsic cortical computations and for understanding how the brain regulates competing inputs. Converging behavioral and EEG results have shown that binocular rivalry and attention are intertwined: binocular rivalry ceases when attention is diverted away from the rivalry stimuli. In addition, the competing image in one eye suppresses the target in the other eye through a pattern of gain changes similar to those induced by attention. These results require a revision of the current computational theories of binocular rivalry, in which the role of attention is ignored. Here, we provide a computational model of binocular rivalry. In the model, competition between two images in rivalry is driven by both attentional modulation and mutual inhibition, which have distinct selectivity (feature vs. eye of origin) and dynamics (relatively slow vs. relatively fast). The proposed model explains a wide range of phenomena reported in rivalry, including the three hallmarks: (i) binocular rivalry requires attention; (ii) various perceptual states emerge when the two images are swapped between the eyes multiple times per second; (iii) the dominance duration as a function of input strength follows Levelt's propositions. With a bifurcation analysis, we identified the parameter space in which the model's behavior was consistent with experimental results.

Order-indeterminant event-based maps for learning a beat.

The process by which humans synchronize to a musical beat is believed to occur through error-correction where an individual's estimates of the period and phase of the beat time are iteratively adjusted to align with an external stimuli. Mathematically, error-correction can be described using a two-dimensional map where convergence to a fixed point corresponds to synchronizing to the beat. In this paper, we show how a neural system, called a beat generator, learns to adapt its oscillatory behavior through error-correction to synchronize to an external periodic signal. We construct a two-dimensional event-based map, which iteratively adjusts an internal parameter of the beat generator to speed up or slow down its oscillatory behavior to bring it into synchrony with the periodic stimulus. The map is novel in that the order of events defining the map are not a priori known. Instead, the type of error-correction adjustment made at each iterate of the map is determined by a sequence of expected events. The map possesses a rich repertoire of dynamics, including periodic solutions and chaotic orbits.

Firing rate models for gamma oscillations.

Gamma oscillations are readily observed in a variety of brain regions during both waking and sleeping states. Computational models of gamma oscillations typically involve simulations of large networks of synaptically coupled spiking units. These networks can exhibit strongly synchronized gamma behavior, whereby neurons fire in near synchrony on every cycle, or weakly modulated gamma behavior, corresponding to stochastic, sparse firing of the individual units on each cycle of the population gamma rhythm. These spiking models offer valuable biophysical descriptions of gamma oscillations; however, because they involve many individual neuronal units they are limited in their ability to communicate general network-level dynamics. Here we demonstrate that few-variable firing rate models with established synaptic timescales can account for both strongly synchronized and weakly modulated gamma oscillations. These models go beyond the classical formulations of rate models by including at least two dynamic variables per population: firing rate and synaptic activation. The models' flexibility to capture the broad range of gamma behavior depends directly on the timescales that represent recruitment of the excitatory and inhibitory firing rates. In particular, we find that weakly modulated gamma oscillations occur robustly when the recruitment timescale of inhibition is faster than that of excitation. We present our findings by using an extended Wilson-Cowan model and a rate model derived from a network of quadratic integrate-and-fire neurons. These biophysical rate models capture the range of weakly modulated and coherent gamma oscillations observed in spiking network models, while additionally allowing for greater tractability and systems analysis. NEW & NOTEWORTHY Here we develop simple and tractable models of gamma oscillations, a dynamic feature observed throughout much of the brain with significant correlates to behavior and cognitive performance in a variety of experimental contexts. Our models depend on only a few dynamic variables per population, but despite this they qualitatively capture features observed in previous biophysical models of gamma oscillations that involve many individual spiking units.

Function and energy consumption constrain neuronal biophysics in a canonical computation: Coincidence detection.

Neural morphology and membrane properties vary greatly between cell types in the nervous system. The computations and local circuit connectivity that neurons support are thought to be the key factors constraining the cells' biophysical properties. Nevertheless, additional constraints can be expected to further shape neuronal design. Here, we focus on a particularly energy-intense system (as indicated by metabolic markers): principal neurons in the medial superior olive (MSO) nucleus of the auditory brainstem. Based on a modeling approach, we show that a trade-off between the level of performance of a functionally relevant computation and energy consumption predicts optimal ranges for cell morphology and membrane properties. The biophysical parameters appear most strongly constrained by functional needs, while energy use is minimized as long as function can be maintained. The key factors that determine model performance and energy consumption are 1) the saturation of the synaptic conductance input and 2) the temporal resolution of the postsynaptic signals as they reach the soma, which is largely determined by active membrane properties. MSO cells seem to operate close to pareto optimality, i.e., the trade-off boundary between performance and energy consumption that is formed by the set of optimal models. Good performance for drastically lower costs could in theory be achieved by small neurons without dendrites, as seen in the avian auditory system, pointing to additional constraints for mammalian MSO cells, including their circuit connectivity.

Signatures of Somatic Inhibition and Dendritic Excitation in Auditory Brainstem Field Potentials.

Extracellular voltage recordings (Ve ; field potentials) provide an accessible view of in vivo neural activity, but proper interpretation of field potentials is a long-standing challenge. Computational modeling can aid in identifying neural generators of field potentials. In the auditory brainstem of cats, spatial patterns of sound-evoked Ve can resemble, strikingly, Ve generated by current dipoles. Previously, we developed a biophysically-based model of a binaural brainstem nucleus, the medial superior olive (MSO), that accounts qualitatively for observed dipole-like Ve patterns in sustained responses to monaural tones with frequencies >∼1000 Hz (Goldwyn et al., 2014). We have observed, however, that Ve patterns in cats of both sexes appear more monopole-like for lower-frequency tones. Here, we enhance our theory to accurately reproduce dipole and non-dipole features of Ve responses to monaural tones with frequencies ranging from 600 to 1800 Hz. By applying our model to data, we estimate time courses of paired input currents to MSO neurons. We interpret these inputs as dendrite-targeting excitation and soma-targeting inhibition (the latter contributes non-dipole-like features to Ve responses). Aspects of inferred inputs are consistent with synaptic inputs to MSO neurons including the tendencies of inhibitory inputs to attenuate in response to high-frequency tones and to precede excitatory inputs. Importantly, our updated theory can be tested experimentally by blocking synaptic inputs. MSO neurons perform a critical role in sound localization and binaural hearing. By solving an inverse problem to uncover synaptic inputs from Ve patterns we provide a new perspective on MSO physiology.SIGNIFICANCE STATEMENT Extracellular voltages (field potentials) are a common measure of brain activity. Ideally, one could infer from these data the activity of neurons and synapses that generate field potentials, but this "inverse problem" is not easily solved. We study brainstem field potentials in the region of the medial superior olive (MSO); a critical center in the auditory pathway. These field potentials exhibit distinctive spatial and temporal patterns in response to pure tone sounds. We use mathematical modeling in combination with physiological and anatomical knowledge of MSO neurons to plausibly explain how dendrite-targeting excitation and soma-targeting inhibition generate these field potentials. Inferring putative synaptic currents from field potentials advances our ability to study neural processing of sound in the MSO.

Modeling fast and slow gamma oscillations with interneurons of different subtype.

Experimental and theoretical studies demonstrate that neuronal gamma oscillations crucially depend on interneurons, but current models do not consider the diversity of known interneuron subtypes. Moreover, in CA1 of the hippocampus, experimental evidence indicates the presence of multiple gamma oscillators, two of which may be coordinated by differing interneuron populations. In this article, we show that models of networks with competing interneuron populations with different postsynaptic effects are sufficient to generate, within CA1, distinct oscillatory regimes. We find that strong mutual inhibition between the interneuron populations permits distinct fast and slow gamma states, whereas weak mutual inhibition generates mixed gamma states. We develop idealized firing rate models to illuminate dynamic properties of these competitive gamma networks, and reinforce these concepts with basic spiking models. The models make several explicit predictions about gamma oscillators in CA1. Specifically, interneurons of different subtype phase-lock to different gamma states, and one population of interneurons is silenced and the other active during fast and slow gamma events. Finally, mutual inhibition between interneuron populations is necessary to generate distinct gamma states. Previous experimental studies indicate that fast and slow gamma oscillations reflect different information processing modes, although it is unclear whether these rhythms are intrinsic or imposed. The models outlined demonstrate that basic architectures can locally generate these oscillations, as well as capture other features of fast and slow gamma, including theta-phase preference and spontaneous transitions between gamma states. These models may extend to describe general dynamics in networks with diverse interneuron populations.NEW & NOTEWORTHY The oscillatory coordination of neural signals is crucial to healthy brain function. We have developed an idealized neuronal model that generates distinct fast and slow gamma oscillations, a known feature of the rodent hippocampus. Our work provides a mechanism of this phenomenon, as well as a theoretical framework for future experiments concerning hippocampal gamma. It moreover offers a tractable model of competitive gamma oscillations that is generalizable across the nervous system.

High-Frequency Resonance in the Gerbil Medial Superior Olive.

A high-frequency, subthreshold resonance in the guinea pig medial superior olive (MSO) was recently linked to the efficient extraction of spatial cues from the fine structure of acoustic stimuli. We report here that MSO neurons in gerbil also have resonant properties and, based on our whole-cell recordings and computational modeling, that a low-voltage-gated potassium current, IKLT, underlies the resonance. We show that resonance was lost following dynamic clamp replacement of IKLT with a leak conductance and in the model when voltage-gating of IKLT was suppressed. Resonance was characterized using small amplitude sinusoidal stimuli to generate impedance curves as typically done for linear systems analysis. Extending our study into the nonlinear, voltage-dependent regime, we increased stimulus amplitude and found, experimentally and in simulations, that the subthreshold resonant frequency (242Hz for weak stimuli) increased continuously to the resonant frequency for spiking (285Hz). The spike resonance of these phasic-firing (type III excitable) MSO neurons and of the model is of particular interest also because previous studies of resonance typically involved neurons/models (type II excitable, such as the standard Hodgkin-Huxley model) that can fire tonically for steady inputs. To probe more directly how these resonances relate to MSO neurons as slope-detectors, we presented periodic trains of brief, fast-rising excitatory post-synaptic potentials (EPSCs) to the model. While weak subthreshold EPSC trains were essentially low-pass filtered, resonance emerged as EPSC amplitude increased. Interestingly, for spike-evoking EPSC trains, the threshold amplitude at spike resonant frequency (317Hz) was lower than the single ESPC threshold. Our finding of a frequency-dependent threshold for repetitive brief EPSC stimuli and preferred frequency for spiking calls for further consideration of both subthreshold and suprathreshold resonance to fast and precise temporal processing in the MSO.

Subthreshold resonance properties contribute to the efficient coding of auditory spatial cues.

Neurons in the medial superior olive (MSO) and lateral superior olive (LSO) of the auditory brainstem code for sound-source location in the horizontal plane, extracting interaural time differences (ITDs) from the stimulus fine structure and interaural level differences (ILDs) from the stimulus envelope. Here, we demonstrate a postsynaptic gradient in temporal processing properties across the presumed tonotopic axis; neurons in the MSO and the low-frequency limb of the LSO exhibit fast intrinsic electrical resonances and low input impedances, consistent with their processing of ITDs in the temporal fine structure. Neurons in the high-frequency limb of the LSO show low-pass electrical properties, indicating they are better suited to extracting information from the slower, modulated envelopes of sounds. Using a modeling approach, we assess ITD and ILD sensitivity of the neural filters to natural sounds, demonstrating that the transformation in temporal processing along the tonotopic axis contributes to efficient extraction of auditory spatial cues.

Neuronal coupling by endogenous electric fields: cable theory and applications to coincidence detector neurons in the auditory brain stem.

The ongoing activity of neurons generates a spatially and time-varying field of extracellular voltage (Ve). This Ve field reflects population-level neural activity, but does it modulate neural dynamics and the function of neural circuits? We provide a cable theory framework to study how a bundle of model neurons generates Ve and how this Ve feeds back and influences membrane potential (Vm). We find that these "ephaptic interactions" are small but not negligible. The model neural population can generate Ve with millivolt-scale amplitude, and this Ve perturbs the Vm of "nearby" cables and effectively increases their electrotonic length. After using passive cable theory to systematically study ephaptic coupling, we explore a test case: the medial superior olive (MSO) in the auditory brain stem. The MSO is a possible locus of ephaptic interactions: sounds evoke large (millivolt scale)Vein vivo in this nucleus. The Ve response is thought to be generated by MSO neurons that perform a known neuronal computation with submillisecond temporal precision (coincidence detection to encode sound source location). Using a biophysically based model of MSO neurons, we find millivolt-scale ephaptic interactions consistent with the passive cable theory results. These subtle membrane potential perturbations induce changes in spike initiation threshold, spike time synchrony, and time difference sensitivity. These results suggest that ephaptic coupling may influence MSO function.

Neuromechanistic Model of Auditory Bistability.

Sequences of higher frequency A and lower frequency B tones repeating in an ABA- triplet pattern are widely used to study auditory streaming. One may experience either an integrated percept, a single ABA-ABA- stream, or a segregated percept, separate but simultaneous streams A-A-A-A- and -B---B--. During minutes-long presentations, subjects may report irregular alternations between these interpretations. We combine neuromechanistic modeling and psychoacoustic experiments to study these persistent alternations and to characterize the effects of manipulating stimulus parameters. Unlike many phenomenological models with abstract, percept-specific competition and fixed inputs, our network model comprises neuronal units with sensory feature dependent inputs that mimic the pulsatile-like A1 responses to tones in the ABA- triplets. It embodies a neuronal computation for percept competition thought to occur beyond primary auditory cortex (A1). Mutual inhibition, adaptation and noise are implemented. We include slow NDMA recurrent excitation for local temporal memory that enables linkage across sound gaps from one triplet to the next. Percepts in our model are identified in the firing patterns of the neuronal units. We predict with the model that manipulations of the frequency difference between tones A and B should affect the dominance durations of the stronger percept, the one dominant a larger fraction of time, more than those of the weaker percept-a property that has been previously established and generalized across several visual bistable paradigms. We confirm the qualitative prediction with our psychoacoustic experiments and use the behavioral data to further constrain and improve the model, achieving quantitative agreement between experimental and modeling results. Our work and model provide a platform that can be extended to consider other stimulus conditions, including the effects of context and volition.

A model of the medial superior olive explains spatiotemporal features of local field potentials.

Local field potentials are important indicators of in vivo neural activity. Sustained, phase-locked, sound-evoked extracellular fields in the mammalian auditory brainstem, known as the auditory neurophonic, reflect the activity of neurons in the medial superior olive (MSO). We develop a biophysically based model of the neurophonic that accounts for features of in vivo extracellular recordings in the cat auditory brainstem. By making plausible idealizations regarding the spatial symmetry of MSO neurons and the temporal synchrony of their afferent inputs, we reduce the challenging problem of computing extracellular potentials in a 3D volume conductor to a one-dimensional problem. We find that postsynaptic currents in bipolar MSO neuron models generate extracellular voltage responses that strikingly resemble in vivo recordings. Simulations reproduce distinctive spatiotemporal features of the in vivo neurophonic response to monaural pure tones: large oscillations (hundreds of microvolts to millivolts), broad spatial reach (millimeter scale), and a dipole-like spatial profile. We also explain how somatic inhibition and the relative timing of bilateral excitation may shape the spatial profile of the neurophonic. We observe in simulations, and find supporting evidence in in vivo data, that coincident excitatory inputs on both dendrites lead to a drastically reduced spatial reach of the neurophonic. This outcome surprises because coincident inputs are thought to evoke maximal firing rates in MSO neurons, and it reconciles previously puzzling evoked potential results in humans and animals. The success of our model, which has no axon or spike-generating sodium currents, suggests that MSO spikes do not contribute appreciably to the neurophonic.