Mechanism of gain modulation at single neuron and network levels

Gain modulation, in which the sensitivity of a neural response to one input is modified by a second input, is studied at single-neuron and network levels. At the single neuron level, gain modulation can arise if the two inputs are subject to a direct multiplicative interaction. Alternatively, these inputs can be summed in a linear manner by the neuron and gain modulation can arise, instead, from a nonlinear input–output relationship. We derive a mathematical constraint that can distinguish these two mechanisms even though they can look very similar, provided sufficient data of the appropriate type are available. Previously, it has been shown in coordinate transformation studies that artificial neurons with sigmoid transfer functions can acquire a nonlinear additive form of gain modulation through learning-driven adjustment of synaptic weights. We use the constraint derived for single-neuron studies to compare responses in this network with those of another network model based on a biologically inspired transfer function that can support approximately multiplicative interactions.     

Dynamical estimation of neuron and network properties II: path integral Monte Carlo methods

Hodgkin–Huxley (HH) models of neuronal membrane dynamics consist of a set of nonlinear differential equations that describe the time-varying conductance of various ion channels. Using observations of voltage alone we show how to estimate the unknown parameters and unobserved state variables of an HH model in the expected circumstance that the measurements are noisy, the model has errors, and the state of the neuron is not known when observations commence. The joint probability distribution of the observed membrane voltage and the unobserved state variables and parameters of these models is a path integral through the model state space. The solution to this integral allows estimation of the parameters and thus a characterization of many biological properties of interest, including channel complement and density, that give rise to a neuron’s electrophysiological behavior. This paper describes a method for directly evaluating the path integral using a Monte Carlo numerical approach. This provides estimates not only of the expected values of model parameters but also of their posterior uncertainty. Using test data simulated from neuronal models comprising several common channels, we show that short (<50 ms) intracellular recordings from neurons stimulated with a complex time-varying current yield accurate and precise estimates of the model parameters as well as accurate predictions of the future behavior of the neuron. We also show that this method is robust to errors in model specification, supporting model development for biological preparations in which the channel expression and other biophysical properties of the neurons are not fully known.     

A hemicord locomotor network of excitatory interneurons: a simulation study

Locomotor burst generation is simulated using a full-scale network model of the unilateral excitatory interneuronal population. Earlier small-scale models predicted that a population of excitatory neurons would be sufficient to produce burst activity, and this has recently been experimentally confirmed. Here we simulate the hemicord activity induced under various experimental conditions, including pharmacological activation by NMDA and AMPA as well as electrical stimulation. The model network comprises a realistic number of cells and synaptic connectivity patterns. Using similar distributions of cellular and synaptic parameters, as have been estimated experimentally, a large variation in dynamic characteristics like firing rates, burst, and cycle durations were seen in single cells. On the network level an overall rhythm was generated because the synaptic interactions cause partial synchronization within the population. This network rhythm not only emerged despite the distributed cellular parameters but relied on this variability, in particular, in reproducing variations of the activity during the cycle and showing recruitment in interneuronal populations. A slow rhythm (0.4–2 Hz) can be induced by tonic activation of NMDA-sensitive channels, which are voltage dependent and generate depolarizing plateaus. The rhythm emerges through a synchronization of bursts of the individual neurons. A fast rhythm (4–12 Hz), induced by AMPA, relies on spike synchronization within the population, and each burst is composed of single spikes produced by different neurons. The dynamic range of the fast rhythm is limited by the ability of the network to synchronize oscillations and depends on the strength of synaptic connections and the duration of the slow after hyperpolarization. The model network also produces prolonged bouts of rhythmic activity in response to brief electrical activations, as seen experimentally. The mutual excitation can sustain long-lasting activity for a realistic set of synaptic parameters. The bout duration depends on the strength of excitatory synaptic connections, the level of persistent depolarization, and the influx of Ca²⁺ ions and activation of Ca²⁺-dependent K⁺ current.     

Dynamical estimation of neuron and network properties I: variational methods

We present a method for using measurements of membrane voltage in individual neurons to estimate the parameters and states of the voltage-gated ion channels underlying the dynamics of the neuron's behavior. Short injections of a complex time-varying current provide sufficient data to determine the reversal potentials, maximal conductances, and kinetic parameters of a diverse range of channels, representing tens of unknown parameters and many gating variables in a model of the neuron's behavior. These estimates are used to predict the response of the model at times beyond the observation window. This method of [Formula: see text] extends to the general problem of determining model parameters and unobserved state variables from a sparse set of observations, and may be applicable to networks of neurons. We describe an exact formulation of the tasks in nonlinear data assimilation when one has noisy data, errors in the models, and incomplete information about the state of the system when observations commence. This is a high dimensional integral along the path of the model state through the observation window. In this article, a stationary path approximation to this integral, using a variational method, is described and tested employing data generated using neuronal models comprising several common channels with Hodgkin-Huxley dynamics. These numerical experiments reveal a number of practical considerations in designing stimulus currents and in determining model consistency. The tools explored here are computationally efficient and have paths to parallelization that should allow large individual neuron and network problems to be addressed.     

Bursting Reverberation as a Multiscale Neuronal Network Process Driven by Synaptic Depression-Facilitation

Neuronal networks can generate complex patterns of activity that depend on membrane properties of individual neurons as well as on functional synapses. To decipher the impact of synaptic properties and connectivity on neuronal network behavior, we investigate the responses of neuronal ensembles from small (5–30 cells in a restricted sphere) and large (acute hippocampal slice) networks to single electrical stimulation: in both cases, a single stimulus generated a synchronous long-lasting bursting activity. While an initial spike triggered a reverberating network activity that lasted 2–5 seconds for small networks, we found here that it lasted only up to 300 milliseconds in slices. To explain this phenomena present at different scales, we generalize the depression-facilitation model and extracted the network time constants. The model predicts that the reverberation time has a bell shaped relation with the synaptic density, revealing that the bursting time cannot exceed a maximum value. Furthermore, before reaching its maximum, the reverberation time increases sub-linearly with the synaptic density of the network. We conclude that synaptic dynamics and connectivity shape the mean burst duration, a property present at various scales of the networks. Thus bursting reverberation is a property of sufficiently connected neural networks, and can be generated by collective depression and facilitation of underlying functional synapses.     

Astrocyte- neuron interaction as a mechanism responsible for generation of neural synchrony: a study based on modeling and experiments

Neural synchronization is considered as an important mechanism for information processing. In addition, based on recent neurophysiologic findings, it is believed that astrocytes regulate the synaptic transmission of neuronal networks. Therefore, the present study focused on determining the functional contribution of astrocytes in neuronal synchrony using both computer simulations and extracellular field potential recordings. For computer simulations, as a first step, a minimal network model is constructed by connecting two Morris-Lecar neuronal models. In this minimal model, astrocyte-neuron interactions are considered in a functional-based procedure. Next, the minimal network is extended and a biologically plausible neuronal population model is developed which considers functional outcome of astrocyte-neuron interactions too. The employed structure is based on the physiological and anatomical network properties of the hippocampal CA1 area. Utilizing these two different levels of modeling, it is demonstrated that astrocytes are able to change the threshold value of transition from synchronous to asynchronous behavior among neurons. In this way, variations in the interaction between astrocytes and neurons lead to the emergence of synchronous/asynchronous patterns in neural responses. Furthermore, population spikes are recorded from CA1 pyramidal neurons in rat hippocampal slices to validate the modeling results. It demonstrates that astrocytes play a primary role in neuronal firing synchronicity and synaptic coordination. These results may offer a new insight into understanding the mechanism by which astrocytes contribute to stabilizing neural activities.     

Clique of Functional Hubs Orchestrates Population Bursts in Developmentally Regulated Neural Networks

It has recently been discovered that single neuron stimulation can impact network dynamics in immature and adult neuronal circuits. Here we report a novel mechanism which can explain in neuronal circuits, at an early stage of development, the peculiar role played by a few specific neurons in promoting/arresting the population activity. For this purpose, we consider a standard neuronal network model, with short-term synaptic plasticity, whose population activity is characterized by bursting behavior. The addition of developmentally inspired constraints and correlations in the distribution of the neuronal connectivities and excitabilities leads to the emergence of functional hub neurons, whose stimulation/deletion is critical for the network activity. Functional hubs form a clique, where a precise sequential activation of the neurons is essential to ignite collective events without any need for a specific topological architecture. Unsupervised time-lagged firings of supra-threshold cells, in connection with coordinated entrainments of near-threshold neurons, are the key ingredients to orchestrate population activity.     

Inferring network activity from synaptic noise

During intense network activity in vivo, cortical neurons are in a high-conductance state, in which the membrane potential (V(m)) is subject to a tremendous fluctuating activity. Clearly, this "synaptic noise" contains information about the activity of the network, but there are presently no methods available to extract this information. We focus here on this problem from a computational neuroscience perspective, with the aim of drawing methods to analyze experimental data. We start from models of cortical neurons, in which high-conductance states stem from the random release of thousands of excitatory and inhibitory synapses. This highly complex system can be simplified by using global synaptic conductances described by effective stochastic processes. The advantage of this approach is that one can derive analytically a number of properties from the statistics of resulting V(m) fluctuations. For example, the global excitatory and inhibitory conductances can be extracted from synaptic noise, and can be related to the mean activity of presynaptic neurons. We show here that extracting the variances of excitatory and inhibitory synaptic conductances can provide estimates of the mean temporal correlation-or level of synchrony-among thousands of neurons in the network. Thus, "probing the network" through intracellular V(m) activity is possible and constitutes a promising approach, but it will require a continuous effort combining theory, computational models and intracellular physiology.     

Consistent dynamics suggests tight regulation of biophysical parameters in a small network of bursting neurons

The neuronal firing patterns in the pyloric network of crustaceans are remarkably consistent among animals. Although this characteristic of the pyloric network is well-known, the biophysical mechanisms underlying the regulation of the systems output are receiving renewed attention. Computer simulations of the pyloric network recently demonstrated that consistent motor output can be achieved from neurons with disparate biophysical parameters among animals. Here we address this hypothesis by pharmacologically manipulating the pyloric network and analyzing the emerging voltage oscillations and firing patterns. Our results show that the pyloric network of the lobster stomatogastric ganglion maintains consistent and regular firing patterns even when entire populations of specific voltage-gated channels and synaptic receptors are blocked. The variations of temporal parameters used to characterize the burst patterns of the neurons as well as their intraburst spike dynamics do not display statistically significant increase after blocking the transient K-currents (with 4-aminopyridine), the glutamatergic inhibitory synapses (with picrotoxin), or the cholinergic synapses (with atropine) in pyloric networks from different animals. These data suggest that in this very compact circuit, the biophysical parameters are cell-specific and tightly regulated.     

Identifying and Tracking Simulated Synaptic Inputs from Neuronal Firing: Insights from In Vitro Experiments

Accurately describing synaptic interactions between neurons and how interactions change over time are key challenges for systems neuroscience. Although intracellular electrophysiology is a powerful tool for studying synaptic integration and plasticity, it is limited by the small number of neurons that can be recorded simultaneously in vitro and by the technical difficulty of intracellular recording in vivo. One way around these difficulties may be to use large-scale extracellular recording of spike trains and apply statistical methods to model and infer functional connections between neurons. These techniques have the potential to reveal large-scale connectivity structure based on the spike timing alone. However, the interpretation of functional connectivity is often approximate, since only a small fraction of presynaptic inputs are typically observed. Here we use in vitro current injection in layer 2/3 pyramidal neurons to validate methods for inferring functional connectivity in a setting where input to the neuron is controlled. In experiments with partially-defined input, we inject a single simulated input with known amplitude on a background of fluctuating noise. In a fully-defined input paradigm, we then control the synaptic weights and timing of many simulated presynaptic neurons. By analyzing the firing of neurons in response to these artificial inputs, we ask 1) How does functional connectivity inferred from spikes relate to simulated synaptic input? and 2) What are the limitations of connectivity inference? We find that individual current-based synaptic inputs are detectable over a broad range of amplitudes and conditions. Detectability depends on input amplitude and output firing rate, and excitatory inputs are detected more readily than inhibitory. Moreover, as we model increasing numbers of presynaptic inputs, we are able to estimate connection strengths more accurately and detect the presence of connections more quickly. These results illustrate the possibilities and outline the limits of inferring synaptic input from spikes.     

Dynamics of spontaneous activity in random networks with multiple neuron subtypes and synaptic noise

During slow-wave sleep, brain electrical activity is dominated by the slow (< 1 Hz) electroencephalogram (EEG) oscillations characterized by the periodic transitions between active (or Up) and silent (or Down) states in the membrane voltage of the cortical and thalamic neurons. Sleep slow oscillation is believed to play critical role in consolidation of recent memories. Past computational studies, based on the Hodgkin-Huxley type neuronal models, revealed possible intracellular and network mechanisms of the neuronal activity during sleep, however, they failed to explore the large-scale cortical network dynamics depending on collective behavior in the large populations of neurons. In this new study, we developed a novel class of reduced discrete time spiking neuron models for large-scale network simulations of wake and sleep dynamics. In addition to the spiking mechanism, the new model implemented nonlinearities capturing effects of the leak current, the Ca²⁺ dependent K⁺ current and the persistent Na⁺ current that were found to be critical for transitions between Up and Down states of the slow oscillation. We applied the new model to study large-scale two-dimensional cortical network activity during slow-wave sleep. Our study explained traveling wave dynamics and characteristic synchronization properties of transitions between Up and Down states of the slow oscillation as observed in vivo in recordings from cats. We further predict a critical role of synaptic noise and slow adaptive currents for spike sequence replay as found during sleep related memory consolidation.     

Recent advances in the use of neurotropic viruses for circuit analysis

In the past decade, the literature describing how viruses can be used to define the synaptic architecture of the brain has increased dramatically. Early studies focused on determining the specificity of viral transport through synaptic connections. Analysis of the assembly and intracellular transport of viruses as well as of the role of the brain response to infection were central to this literature. With the growing acceptance that the transport of viruses is circuit-related, attention has shifted to application of the method to define the functional architecture of neural systems. The development of attenuated recombinant viruses that maintain neuroinvasiveness has been instrumental to the generation of increasingly powerful experimental approaches for the functional dissection of neural circuits. These approaches include the use of recombinant viruses that express unique reporters to address issues of axon collateralization in complex circuits, the use of green-fluorescent-protein-expressing recombinants to characterize the electrophysiological properties of projection-specific neurons in live slices of brain, and the exploitation of the Cre recombinase system for conditional replication of virus in phenotypically defined populations of neurons.     

Mimicking synaptic effects of addictive drugs with selective dopamine neuron stimulation

The synaptic changes induced by initial drug exposure leave a trace on neural systems that can eventually manifest in compulsive drug-seeking behavior. A single injection of cocaine has been shown to induce a change in the AMPA receptor (AMPAR) subunit composition at glutamatergic synapses onto ventral tegmental area (VTA) dopamine (DA) neurons. This change is long-lasting (up to months following self-administration) and represents an important functional change at the synaptic level following cocaine use. We recently published findings that cocaine's action at the DA transporter (DAT) is necessary for the induction of AMPAR redistribution and that this can also be mimicked by selective DA neuron stimulation. The stimulation effect is dependent on D1 receptors within the VTA. Furthermore other addictive drugs, although they act through distinct mechanisms, also induce this synaptic change. Here we discuss literature that expands on these observations in an attempt to further clarify the synaptic changes following early drug use.     

锥体神经元的连接模式对突触后局部电位的依赖

脑皮层的功能连接模式与突触可塑性密切相关,受突触空间分布和刺激模式等多种因素的影响。尽管越来越多的证据表明突触可塑性不仅受突触后动作电位而且还受突触后局部树突电位的影响,但是目前尚不清楚神经元的功能连接模式是否和怎样依赖于突触后局部电位的。为此,本文建立了一个无需硬边界设置的、突触后局部膜电位依赖的可塑性模型。该模型具有突触强度的自平衡能力并且能够再现多种突触可塑性实验结果。基于该模型对两个锥体神经元的功能连接模式进行仿真的结果表明,当突触后局部电位都处于亚阈值时两个神经元无功能连接,如果一个神经元的突触后膜电位高于阈值电位则产生向该神经元的单向连接,当两个神经元的突触后膜电位都超过阈值电位时则产生双向连接,说明突触后局部膜电位分布是神经元功能连接模式形成的关键。研究结果加深了神经网络连接模式形成机制的理解,对学习和记忆的研究具有重要意义。     

Variability v.s. synchronicity of neuronal activity in local cortical network models with different wiring topologies

Dynamical behavior of a biological neuronal network depends significantly on the spatial pattern of synaptic connections among neurons. While neuronal network dynamics has extensively been studied with simple wiring patterns, such as all-to-all or random synaptic connections, not much is known about the activity of networks with more complicated wiring topologies. Here, we examined how different wiring topologies may influence the response properties of neuronal networks, paying attention to irregular spike firing, which is known as a characteristic of in vivo cortical neurons, and spike synchronicity. We constructed a recurrent network model of realistic neurons and systematically rewired the recurrent synapses to change the network topology, from a localized regular and a “small-world” network topology to a distributed random network topology. Regular and small-world wiring patterns greatly increased the irregularity or the coefficient of variation (Cv) of output spike trains, whereas such an increase was small in random connectivity patterns. For given strength of recurrent synapses, the firing irregularity exhibited monotonous decreases from the regular to the random network topology. By contrast, the spike coherence between an arbitrary neuron pair exhibited a non-monotonous dependence on the topological wiring pattern. More precisely, the wiring pattern to maximize the spike coherence varied with the strength of recurrent synapses. In a certain range of the synaptic strength, the spike coherence was maximal in the small-world network topology, and the long-range connections introduced in this wiring changed the dependence of spike synchrony on the synaptic strength moderately. However, the effects of this network topology were not really special in other properties of network activity. Action Editor: Xiao-Jing Wang     

The information theory of developmental pruning: Optimizing global network architectures using local synaptic rules

During development, biological neural networks produce more synapses and neurons than needed. Many of these synapses and neurons are later removed in a process known as neural pruning. Why networks should initially be over-populated, and the processes that determine which synapses and neurons are ultimately pruned, remains unclear. We study the mechanisms and significance of neural pruning in model neural networks. In a deep Boltzmann machine model of sensory encoding, we find that (1) synaptic pruning is necessary to learn efficient network architectures that retain computationally-relevant connections, (2) pruning by synaptic weight alone does not optimize network size and (3) pruning based on a locally-available measure of importance based on Fisher information allows the network to identify structurally important vs. unimportant connections and neurons. This locally-available measure of importance has a biological interpretation in terms of the correlations between presynaptic and postsynaptic neurons, and implies an efficient activity-driven pruning rule. Overall, we show how local activity-dependent synaptic pruning can solve the global problem of optimizing a network architecture. We relate these findings to biology as follows: (I) Synaptic over-production is necessary for activity-dependent connectivity optimization. (II) In networks that have more neurons than needed, cells compete for activity, and only the most important and selective neurons are retained. (III) Cells may also be pruned due to a loss of synapses on their axons. This occurs when the information they convey is not relevant to the target population.     

Asymptotics applied to a neural network

A mathematical model of neural processing is proposed which incorporates a theory for the storage of information. The model consists of a network of neurons that linearly processes incoming neural activity. The network stores the input by modifying the synaptic properties of all of its neurons. The model lends support to a distributive theory of memory using synaptic modification. The dynamics of the processing and storage are represented by a discrete system. Asymptotic analysis is applied to the system to show the learning capabilities of the network under constant input. Results are also given to predict the network's ability to learn periodic input, and input subjected to small random fluctuations.     

The significance of action potential bursting in the brain reward circuit

The brain reward circuit consists of specialized cortical and subcortical structural components that code for various cognitive aspects of goal-directed behavior. These components include the prefrontal cortex (PFC), amygdala (AMY), nucleus accumbens (Nac), subiculum (SUB) of the hippocampal formation, and the dopamine (DA) neurons in the ventral tegmental area (VTA). Both serial and parallel processing in the different components of the circuit code the various aspects of reward-related behavior. Individual neurons within each component have developed specialized intrinsic membrane properties that have led them to be typically defined as either single spiking or high frequency burst-firing neurons. However, a strict definition based on the output mode may not be appropriate. Under the right conditions, neurons can switch between bursting and single-spiking modes, therefore providing a conditional output state. The preferred mode of each individual neuron depends on a combination of different plastic neuronal properties such as, dendritic architecture, neuromodulation, intracellular calcium (Ca(++)) buffering, excitatory and inhibitory synaptic strength, and the spatial distribution and density of voltage and ligand-gated channels. It is likely that, in vivo, most neurons in the circuit, despite variations in intrinsic membrane properties, are conditional output neurons equipped with the versatility of switching between output modes under appropriate conditions. Bursting mode may be used to boost the gain of neural signaling of important or novel events by enhancing transmitter release and enhancing dendritic depolarization, thereby increasing synaptic potentiation. Conversely, single spiking mode may be used to dampen neuronal signaling and may be associated with habituation to unimportant events. Mode switching may provide flexibility to the circuit allowing different sets of neurons to conditionally code for the various aspects of reward-related memory and behavior.     

Estimating conductances of dual-recorded neurons within a network of coupled cells

Simultaneous pre- and postsynaptic cell recordings are used to calculate gap junction conductance based on an equivalent electrical circuit of an electrically coupled pair of cells. This calculation is imprecise when recording from a cell pair that is coupled to neighboring cells providing indirect conductance paths between the recorded cells. Despite this imprecision, junctional conductance has been calculated for coupled cell networks during the past 40 years since a more accurate method was lacking. The present study simulated a three-dimensional network of electrically coupled heterogeneous neurons and used mathematical modeling to reduce the complexity to the simplest equations that could more accurately estimate the electrical properties of dual-recorded cells in the network. Analyses of the simulations showed that knowledge of the number of unrecorded cells directly linked to the recorded cells and of the voltage responses of these recorded cells were largely sufficient to accurately predict the direct junctional resistance linking the recorded cells as well as the input resistance of the recorded cells that would exist in the absence of junctional coupling. All model parameters could be obtained from real dual-intracellular penetrations which allow electrophysiological recordings and intracellular staining.     

Computation of the electroencephalogram (EEG) from network models of point neurons

The electroencephalogram (EEG) is a major tool for non-invasively studying brain function and dysfunction. Comparing experimentally recorded EEGs with neural network models is important to better interpret EEGs in terms of neural mechanisms. Most current neural network models use networks of simple point neurons. They capture important properties of cortical dynamics, and are numerically or analytically tractable. However, point neurons cannot generate an EEG, as EEG generation requires spatially separated transmembrane currents. Here, we explored how to compute an accurate approximation of a rodent’s EEG with quantities defined in point-neuron network models. We constructed different approximations (or proxies) of the EEG signal that can be computed from networks of leaky integrate-and-fire (LIF) point neurons, such as firing rates, membrane potentials, and combinations of synaptic currents. We then evaluated how well each proxy reconstructed a ground-truth EEG obtained when the synaptic currents of the LIF model network were fed into a three-dimensional network model of multicompartmental neurons with realistic morphologies. Proxies based on linear combinations of AMPA and GABA currents performed better than proxies based on firing rates or membrane potentials. A new class of proxies, based on an optimized linear combination of time-shifted AMPA and GABA currents, provided the most accurate estimate of the EEG over a wide range of network states. The new linear proxies explained 85–95% of the variance of the ground-truth EEG for a wide range of network configurations including different cell morphologies, distributions of presynaptic inputs, positions of the recording electrode, and spatial extensions of the network. Non-linear EEG proxies using a convolutional neural network (CNN) on synaptic currents increased proxy performance by a further 2–8%. Our proxies can be used to easily calculate a biologically realistic EEG signal directly from point-neuron simulations thus facilitating a quantitative comparison between computational models and experimental EEG recordings.