Synaptic homeostasis transiently leverages Hebbian mechanisms for a multiphasic response to inactivity.

Homeostatic regulation of synapses is vital for nervous system function and key to understanding a range of neurological conditions. Synaptic homeostasis is proposed to operate over hours to counteract the destabilizing influence of long-term potentiation (LTP) and long-term depression (LTD). The prevailing view holds that synaptic scaling is a slow first-order process that regulates postsynaptic glutamate receptors and fundamentally differs from LTP or LTD. Surprisingly, we find that the dynamics of scaling induced by neuronal inactivity are not exponential or monotonic, and the mechanism requires calcineurin and CaMKII, molecules dominant in LTD and LTP. Our quantitative model of these enzymes reconstructs the unexpected dynamics of homeostatic scaling and reveals how synapses can efficiently safeguard future capacity for synaptic plasticity. This mechanism of synaptic adaptation supports a broader set of homeostatic changes, including action potential autoregulation, and invites further inquiry into how such a mechanism varies in health and disease.

Dynamic models for musical rhythm perception and coordination.

Rhythmicity permeates large parts of human experience. Humans generate various motor and brain rhythms spanning a range of frequencies. We also experience and synchronize to externally imposed rhythmicity, for example from music and song or from the 24-h light-dark cycles of the sun. In the context of music, humans have the ability to perceive, generate, and anticipate rhythmic structures, for example, "the beat." Experimental and behavioral studies offer clues about the biophysical and neural mechanisms that underlie our rhythmic abilities, and about different brain areas that are involved but many open questions remain. In this paper, we review several theoretical and computational approaches, each centered at different levels of description, that address specific aspects of musical rhythmic generation, perception, attention, perception-action coordination, and learning. We survey methods and results from applications of dynamical systems theory, neuro-mechanistic modeling, and Bayesian inference. Some frameworks rely on synchronization of intrinsic brain rhythms that span the relevant frequency range; some formulations involve real-time adaptation schemes for error-correction to align the phase and frequency of a dedicated circuit; others involve learning and dynamically adjusting expectations to make rhythm tracking predictions. Each of the approaches, while initially designed to answer specific questions, offers the possibility of being integrated into a larger framework that provides insights into our ability to perceive and generate rhythmic patterns.

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.

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.

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.

Attentional control via synaptic gain mechanisms in auditory streaming.

Attention is a crucial component in sound source segregation allowing auditory objects of interest to be both singled out and held in focus. Our study utilizes a fundamental paradigm for sound source segregation: a sequence of interleaved tones, A and B, of different frequencies that can be heard as a single integrated stream or segregated into two streams (auditory streaming paradigm). We focus on the irregular alternations between integrated and segregated that occur for long presentations, so-called auditory bistability. Psychaoustic experiments demonstrate how attentional control, a listener's intention to experience integrated or segregated, biases perception in favour of different perceptual interpretations. Our data show that this is achieved by prolonging the dominance times of the attended percept and, to a lesser extent, by curtailing the dominance times of the unattended percept, an effect that remains consistent across a range of values for the difference in frequency between A and B. An existing neuromechanistic model describes the neural dynamics of perceptual competition downstream of primary auditory cortex (A1). The model allows us to propose plausible neural mechanisms for attentional control, as linked to different attentional strategies, in a direct comparison with behavioural data. A mechanism based on a percept-specific input gain best accounts for the effects of attentional control.

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.

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.

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.

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.

Auditory streaming and bistability paradigm extended to a dynamic environment.

We explore stream segregation with temporally modulated acoustic features using behavioral experiments and modelling. The auditory streaming paradigm in which alternating high- A and low-frequency tones B appear in a repeating ABA-pattern, has been shown to be perceptually bistable for extended presentations (order of minutes). For a fixed, repeating stimulus, perception spontaneously changes (switches) at random times, every 2-15 s, between an integrated interpretation with a galloping rhythm and segregated streams. Streaming in a natural auditory environment requires segregation of auditory objects with features that evolve over time. With the relatively idealized ABA-triplet paradigm, we explore perceptual switching in a non-static environment by considering slowly and periodically varying stimulus features. Our previously published model captures the dynamics of auditory bistability and predicts here how perceptual switches are entrained, tightly locked to the rising and falling phase of modulation. In psychoacoustic experiments we find that entrainment depends on both the period of modulation and the intrinsic switch characteristics of individual listeners. The extended auditory streaming paradigm with slowly modulated stimulus features presented here will be of significant interest for future imaging and neurophysiology experiments by reducing the need for subjective perceptual reports of ongoing perception.

Computational models of auditory perception from feature extraction to stream segregation and behavior.

Audition is by nature dynamic, from brainstem processing on sub-millisecond time scales, to segregating and tracking sound sources with changing features, to the pleasure of listening to music and the satisfaction of getting the beat. We review recent advances from computational models of sound localization, of auditory stream segregation and of beat perception/generation. A wealth of behavioral, electrophysiological and imaging studies shed light on these processes, typically with synthesized sounds having regular temporal structure. Computational models integrate knowledge from different experimental fields and at different levels of description. We advocate a neuromechanistic modeling approach that incorporates knowledge of the auditory system from various fields, that utilizes plausible neural mechanisms, and that bridges our understanding across disciplines.

NREM sleep in the rodent neocortex and hippocampus reflects excitable dynamics.

During non-rapid eye movement (NREM) sleep, neuronal populations in the mammalian forebrain alternate between periods of spiking and inactivity. Termed the slow oscillation in the neocortex and sharp wave-ripples in the hippocampus, these alternations are often considered separately but are both crucial for NREM functions. By directly comparing experimental observations of naturally-sleeping rats with a mean field model of an adapting, recurrent neuronal population, we find that the neocortical alternations reflect a dynamical regime in which a stable active state is interrupted by transient inactive states (slow waves) while the hippocampal alternations reflect a stable inactive state interrupted by transient active states (sharp waves). We propose that during NREM sleep in the rodent, hippocampal and neocortical populations are excitable: each in a stable state from which internal fluctuations or external perturbation can evoke the stereotyped population events that mediate NREM functions.

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.

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.

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.

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.

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.