Effects of neuromodulation in a cortical network model of object working memory dominated by recurrent inhibition

Constructing neural networks with pre-specified dynamics

A main goal in neuroscience is to understand the computations carried out by neural populations that give animals their cognitive skills. Neural network models allow to formulate explicit hypotheses regarding the algorithms instantiated in the dynamics of a neural population, its firing statistics, and the underlying connectivity. Neural networks can be defined by a small set of parameters, carefully chosen to procure specific capabilities, or by a large set of free parameters, fitted with optimization algorithms that minimize a given loss function. In this work we alternatively propose a method to make a detailed adjustment of the network dynamics and firing statistic to better answer questions that link dynamics, structure, and function. Our algorithm-termed generalised Firing-to-Parameter (gFTP)-provides a way to construct binary recurrent neural networks whose dynamics strictly follows a user pre-specified transition graph that details the transitions between population firing states triggered by stimulus presentations. Our main contribution is a procedure that detects when a transition graph is not realisable in terms of a neural network, and makes the necessary modifications in order to obtain a new transition graph that is realisable and preserves all the information encoded in the transitions of the original graph. With a realisable transition graph, gFTP assigns values to the network firing states associated with each node in the graph, and finds the synaptic weight matrices by solving a set of linear separation problems. We test gFTP performance by constructing networks with random dynamics, continuous attractor-like dynamics that encode position in 2-dimensional space, and discrete attractor dynamics. We then show how gFTP can be employed as a tool to explore the link between structure, function, and the algorithms instantiated in the network dynamics.

Incorporating slow NMDA-type receptors with nonlinear voltage-dependent magnesium block in a next generation neural mass model: derivation and dynamics

We derive a next generation neural mass model of a population of quadratic-integrate-and-fire neurons, with slow adaptation, and conductance-based AMPAR, GABAR and nonlinear NMDAR synapses. We show that the Lorentzian ansatz assumption can be satisfied by introducing a piece-wise polynomial approximation of the nonlinear voltage-dependent magnesium block of NMDAR current. We study the dynamics of the resulting system for two example cases of excitatory cortical neurons and inhibitory striatal neurons. Bifurcation diagrams are presented comparing the different dynamical regimes as compared to the case of linear NMDAR currents, along with sample comparison simulation time series demonstrating different possible oscillatory solutions. The omission of the nonlinearity of NMDAR currents results in a shift in the range (and possible disappearance) of the constant high firing rate regime, along with a modulation in the amplitude and frequency power spectrum of oscillations. Moreover, nonlinear NMDAR action is seen to be state-dependent and can have opposite effects depending on the type of neurons involved and the level of input firing rate received. The presented model can serve as a computationally efficient building block in whole brain network models for investigating the differential modulation of different types of synapses under neuromodulatory influence or receptor specific malfunction.

A Basal Ganglia model for understanding working memory functions in healthy and Parkinson's conditions

Working memory (WM) is considered as the scratchpad for reading, writing, and processing information necessary to perform cognitive tasks. The Basal Ganglia (BG) and Prefrontal Cortex are two important parts of the brain that are involved in WM functions, and both structures receive projections from dopaminergic nuclei. In this modelling study, we specifically focus on modelling the WM functions of the BG, the WM deficits in Parkinson's disease (PD) conditions, and the impact of dopamine deficiency on different kinds of WM functions. Though there are many experimental and modelling studies of WM properties, there is a paucity of models of the BG that provide insights into the contributions of the BG in WM functions. The proposed model of BG uses bistable flip-flop neurons to model striatal up-down neurons, a network of nonlinear oscillators to model the oscillations of the Indirect Pathway of BG and race-model for action selection. Five different WM tasks are used to demonstrate the generalisation ability of the proposed model. Experimental data from the four tasks are compared with model performance in both control and PD conditions. The model is extended to predict the response time of subjects and in the PD version of the model, the effect of dopaminergic medication on WM performance is also simulated. The proposed model of BG is a unified model that can explain the WM functions of the BG over a wide variety of tasks in both normal and PD conditions, and can be used to understand why specific WM functions are impaired whereas others remain intact in PD.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11571-023-10056-y.

Mixed selectivity: Cellular computations for complexity

The property of mixed selectivity has been discussed at a computational level and offers a strategy to maximize computational power by adding versatility to the functional role of each neuron. Here, we offer a biologically grounded implementational-level mechanistic explanation for mixed selectivity in neural circuits. We define pure, linear, and nonlinear mixed selectivity and discuss how these response properties can be obtained in simple neural circuits. Neurons that respond to multiple, statistically independent variables display mixed selectivity. If their activity can be expressed as a weighted sum, then they exhibit linear mixed selectivity; otherwise, they exhibit nonlinear mixed selectivity. Neural representations based on diverse nonlinear mixed selectivity are high dimensional; hence, they confer enormous flexibility to a simple downstream readout neural circuit. However, a simple neural circuit cannot possibly encode all possible mixtures of variables simultaneously, as this would require a combinatorially large number of mixed selectivity neurons. Gating mechanisms like oscillations and neuromodulation can solve this problem by dynamically selecting which variables are mixed and transmitted to the readout.

Pairing cellular and synaptic dynamics into building blocks of rhythmic neural circuits. A tutorial

In this study we focus on two subnetworks common in the circuitry of swim central pattern generators (CPGs) in the sea slugs, Melibe leonina and Dendronotus iris and show that they are independently capable of stably producing emergent network bursting. This observation raises the question of whether the coordination of redundant bursting mechanisms plays a role in the generation of rhythm and its regulation in the given swim CPGs. To address this question, we investigate two pairwise rhythm-generating networks and examine the properties of their fundamental components: cellular and synaptic, which are crucial for proper network assembly and its stable function. We perform a slow-fast decomposition analysis of cellular dynamics and highlight its significant bifurcations occurring in isolated and coupled neurons. A novel model for slow synapses with high filtering efficiency and temporal delay is also introduced and examined. Our findings demonstrate the existence of two modes of oscillation in bicellular rhythm-generating networks with network hysteresis: i) a half-center oscillator and ii) an excitatory-inhibitory pair. These 2-cell networks offer potential as common building blocks combined in modular organization of larger neural circuits preserving robust network hysteresis.

The memory systems of the human brain and generative artificial intelligence

Generative Artificial Intelligence foundation models (for example Generative Pre-trained Transformer - GPT - models) can generate the next token given a sequence of tokens. How can this 'generative AI' be compared with the 'real' intelligence of the human brain, when for example a human generates a whole memory in response to an incomplete retrieval cue, and then generates further prospective thoughts? Here these two types of generative intelligence, artificial in machines and real in the human brain are compared, and it is shown how when whole memories are generated by hippocampal recall in response to an incomplete retrieval cue, what the human brain computes, and how it computes it, are very different from generative AI. Key differences are the use of local associative learning rules in the hippocampal memory system, and of non-local backpropagation of error learning in AI. Indeed, it is argued that the whole operation of the human brain is performed computationally very differently to what is implemented in generative AI. Moreover, it is emphasized that the primate including human hippocampal system includes computations about spatial view and where objects and people are in scenes, whereas in rodents the emphasis is on place cells and path integration by movements between places. This comparison with generative memory and processing in the human brain has interesting implications for the further development of generative AI and for neuroscience research.

In vivo whole-cortex marker of excitation-inhibition ratio indexes cortical maturation and cognitive ability in youth

A balanced excitation-inhibition ratio (E/I ratio) is critical for healthy brain function. Normative development of cortex-wide E/I ratio remains unknown. Here, we noninvasively estimate a putative marker of whole-cortex E/I ratio by fitting a large-scale biophysically plausible circuit model to resting-state functional MRI (fMRI) data. We first confirm that our model generates realistic brain dynamics in the Human Connectome Project. Next, we show that the estimated E/I ratio marker is sensitive to the gamma-aminobutyric acid (GABA) agonist benzodiazepine alprazolam during fMRI. Alprazolam-induced E/I changes are spatially consistent with positron emission tomography measurement of benzodiazepine receptor density. We then investigate the relationship between the E/I ratio marker and neurodevelopment. We find that the E/I ratio marker declines heterogeneously across the cerebral cortex during youth, with the greatest reduction occurring in sensorimotor systems relative to association systems. Importantly, among children with the same chronological age, a lower E/I ratio marker (especially in the association cortex) is linked to better cognitive performance. This result is replicated across North American (8.2 to 23.0 y old) and Asian (7.2 to 7.9 y old) cohorts, suggesting that a more mature E/I ratio indexes improved cognition during normative development. Overall, our findings open the door to studying how disrupted E/I trajectories may lead to cognitive dysfunction in psychopathology that emerges during youth.

Spontaneous persistent activity and inactivity in vivo reveals differential cortico-entorhinal functional connectivity

Understanding the functional connectivity between brain regions and its emergent dynamics is a central challenge. Here we present a theory-experiment hybrid approach involving iteration between a minimal computational model and in vivo electrophysiological measurements. Our model not only predicted spontaneous persistent activity (SPA) during Up-Down-State oscillations, but also inactivity (SPI), which has never been reported. These were confirmed in vivo in the membrane potential of neurons, especially from layer 3 of the medial and lateral entorhinal cortices. The data was then used to constrain two free parameters, yielding a unique, experimentally determined model for each neuron. Analytic and computational analysis of the model generated a dozen quantitative predictions about network dynamics, which were all confirmed in vivo to high accuracy. Our technique predicted functional connectivity; e. g. the recurrent excitation is stronger in the medial than lateral entorhinal cortex. This too was confirmed with connectomics data. This technique uncovers how differential cortico-entorhinal dialogue generates SPA and SPI, which could form an energetically efficient working-memory substrate and influence the consolidation of memories during sleep. More broadly, our procedure can reveal the functional connectivity of large networks and a theory of their emergent dynamics.

Stochastic attractor models of visual working memory

This paper investigates models of working memory in which memory traces evolve according to stochastic attractor dynamics. These models have previously been shown to account for response-biases that are manifest across multiple trials of a visual working memory task. Here we adapt this approach by making the stable fixed points correspond to the multiple items to be remembered within a single-trial, in accordance with standard dynamical perspectives of memory, and find evidence that this multi-item model can provide a better account of behavioural data from continuous-report tasks. Additionally, the multi-item model proposes a simple mechanism by which swap-errors arise: memory traces diffuse away from their initial state and are captured by the attractors of other items. Swap-error curves reveal the evolution of this process as a continuous function of time throughout the maintenance interval and can be inferred from experimental data. Consistent with previous findings, we find that empirical memory performance is not well characterised by a purely-diffusive process but rather by a stochastic process that also embodies error-correcting dynamics.

A general description of criticality in neural network models

Recent experimental observations have supported the hypothesis that the cerebral cortex operates in a dynamical regime near criticality, where the neuronal network exhibits a mixture of ordered and disordered patterns. However, A comprehensive study of how criticality emerges and how to reproduce it is still lacking. In this study, we investigate coupled networks with conductance-based neurons and illustrate the co-existence of different spiking patterns, including asynchronous irregular (AI) firing and synchronous regular (SR) state, along with a scale-invariant neuronal avalanche phenomenon (criticality). We show that fast-acting synaptic coupling can evoke neuronal avalanches in the mean-dominated regime but has little effect in the fluctuation-dominated regime. In a narrow region of parameter space, the network exhibits avalanche dynamics with power-law avalanche size and duration distributions. We conclude that three stages which may be responsible for reproducing the synchronized bursting: mean-dominated subthreshold dynamics, fast-initiating a spike event, and time-delayed inhibitory cancellation. Remarkably, we illustrate the mechanisms underlying critical avalanches in the presence of noise, which can be explained as a stochastic crossing state around the Hopf bifurcation under the mean-dominated regime. Moreover, we apply the ensemble Kalman filter to determine and track effective connections for the neuronal network. The method is validated on noisy synthetic BOLD signals and could exactly reproduce the corresponding critical network activity. Our results provide a special perspective to understand and model the criticality, which can be useful for large-scale modeling and computation of brain dynamics.

On a framework of data assimilation for hyperparameter estimation of spiking neuronal networks

When handling real-world data modeled by a complex network dynamical system, the number of the parameters is often much more than the size of the data. Therefore, in many cases, it is impossible to estimate these parameters and the exact value of each parameter is frequently less interesting than the distribution of the parameters. In this paper, we aim to estimate the distribution of the parameters in the mesoscopic neuronal network model from the macroscopic experimental data, for example, the BOLD (blood oxygen level dependent) signal. Herein, we assume that the parameters of the neurons and synapses are inhomogeneous but independently and identically distributed from certain distributions with unknown hyperparameters. Thus, we estimate these hyperparameters of the distributions of the parameters, instead of estimating the parameters themselves. We formulate this problem under the framework of data assimilation and hierarchical Bayesian method and present an efficient method named Hierarchical Data Assimilation (HDA) to conduct the statistical inference on the neuronal network model with the BOLD signal data simulated by the hemodynamic model. We consider the Leaky Integral-Fire (LIF) neuronal networks with four synapses and show that the proposed algorithm can estimate the BOLD signals and the hyperparameters with high preciseness. In addition, we discuss the influence on the performance of the algorithm configuration and the LIF network model setup.

Exploring the dynamic interplay between learning and working memory within various cognitive contexts

Introduction

The intertwined relationship between reinforcement learning and working memory in the brain is a complex subject, widely studied across various domains in neuroscience. Research efforts have focused on identifying the specific brain areas responsible for these functions, understanding their contributions in accomplishing the related tasks, and exploring their adaptability under conditions such as cognitive impairment or aging.

Methods

Numerous models have been introduced to formulate either these two subsystems of reinforcement learning and working memory separately or their combination and relationship in executing cognitive tasks. This study adopts the RLWM model as a computational framework to analyze the behavioral parameters of subjects with varying cognitive abilities due to age or cognitive status. A related RLWM task is employed to assess a group of subjects across different age groups and cognitive abilities, as measured by the Montreal Cognitive Assessment tool (MoCA).

Results

Analysis reveals a decline in overall performance accuracy and speed with differing age groups (young vs. middle-aged). Significant differences are observed in model parameters such as learning rate, WM decay, and decision noise. Furthermore, among the middle-aged group, distinctions emerge between subjects categorized as normal vs. MCI based on MoCA scores, notably in speed, performance accuracy, and decision noise.

A prefrontal network model operating near steady and oscillatory states links spike desynchronization and synaptic deficits in schizophrenia

Schizophrenia results in part from a failure of prefrontal networks but we lack full understanding of how disruptions at a synaptic level cause failures at the network level. This is a crucial gap in our understanding because it prevents us from discovering how genetic mutations and environmental risks that alter synaptic function cause prefrontal network to fail in schizophrenia. To address that question, we developed a recurrent spiking network model of prefrontal local circuits that can explain the link between NMDAR synaptic and 0-lag spike synchrony deficits we recently observed in a pharmacological monkey model of prefrontal network failure in schizophrenia. We analyze how the balance between AMPA and NMDA components of recurrent excitation and GABA inhibition in the network influence oscillatory spike synchrony to inform the biological data. We show that reducing recurrent NMDAR synaptic currents prevents the network from shifting from a steady to oscillatory state in response to extrinsic inputs such as might occur during behavior. These findings strongly parallel dynamic modulation of 0-lag spike synchrony we observed between neurons in monkey prefrontal cortex during behavior, as well as the suppression of this 0-lag spiking by administration of NMDAR antagonists. As such, our cortical network model provides a plausible mechanism explaining the link between NMDAR synaptic and 0-lag spike synchrony deficits observed in a pharmacological monkey model of prefrontal network failure in schizophrenia.

Mean-field method for generic conductance-based integrate-and-fire neurons with finite timescales

The construction of transfer functions in theoretical neuroscience plays an important role in determining the spiking rate behavior of neurons in networks. These functions can be obtained through various fitting methods, but the biological relevance of the parameters is not always clear. However, for stationary inputs, such functions can be obtained without the adjustment of free parameters by using mean-field methods. In this work, we expand current Fokker-Planck approaches to account for the concurrent influence of colored and multiplicative noise terms on generic conductance-based integrate-and-fire neurons. We reduce the resulting stochastic system through the application of the diffusion approximation to a one-dimensional Langevin equation. An effective Fokker-Planck is then constructed using Fox Theory, which is solved numerically using a newly developed double integration procedure to obtain the transfer function and the membrane potential distribution. The solution is capable of reproducing the transfer function and the stationary voltage distribution of simulated neurons across a wide range of parameters. The method can also be easily extended to account for different sources of noise with various multiplicative terms, and it can be used in other types of problems in principle.

Information representation in an oscillating neural field model modulated by working memory signals

We study how stimulus information can be represented in the dynamical signatures of an oscillatory model of neural activity-a model whose activity can be modulated by input akin to signals involved in working memory (WM). We developed a neural field model, tuned near an oscillatory instability, in which the WM-like input can modulate the frequency and amplitude of the oscillation. Our neural field model has a spatial-like domain in which an input that preferentially targets a point-a stimulus feature-on the domain will induce feature-specific activity changes. These feature-specific activity changes affect both the mean rate of spikes and the relative timing of spiking activity to the global field oscillation-the phase of the spiking activity. From these two dynamical signatures, we define both a spike rate code and an oscillatory phase code. We assess the performance of these two codes to discriminate stimulus features using an information-theoretic analysis. We show that global WM input modulations can enhance phase code discrimination while simultaneously reducing rate code discrimination. Moreover, we find that the phase code performance is roughly two orders of magnitude larger than that of the rate code defined for the same model solutions. The results of our model have applications to sensory areas of the brain, to which prefrontal areas send inputs reflecting the content of WM. These WM inputs to sensory areas have been established to induce oscillatory changes similar to our model. Our model results suggest a mechanism by which WM signals may enhance sensory information represented in oscillatory activity beyond the comparatively weak representations based on the mean rate activity.

Static and dynamic coding in distinct cell types during associative learning in the prefrontal cortex

The prefrontal cortex maintains information in memory through static or dynamic population codes depending on task demands, but whether the population coding schemes used are learning-dependent and differ between cell types is currently unknown. We investigate the population coding properties and temporal stability of neurons recorded from male macaques in two mapping tasks during and after stimulus-response associative learning, and then we use a Strategy task with the same stimuli and responses as control. We identify a heterogeneous population coding for stimuli, responses, and novel associations: static for putative pyramidal cells and dynamic for putative interneurons that show the strongest selectivity for all the variables. The population coding of learned associations shows overall the highest stability driven by cell types, with interneurons changing from dynamic to static coding after successful learning. The results support that prefrontal microcircuitry expresses mixed population coding governed by cell types and changes its stability during associative learning.

Reduced evoked cortical beta and gamma activity and neuronal synchronization in succinic semialdehyde dehydrogenase deficiency, a disorder of γ-aminobutyric acid metabolism

Succinic semialdehyde dehydrogenase deficiency is a rare autosomal recessively inherited metabolic disorder of γ-aminobutyric acid catabolism manifested by intellectual disability, expressive aphasia, movement disorders, psychiatric ailments and epilepsy. Subjects with succinic semialdehyde dehydrogenase deficiency are characterized by elevated γ-aminobutyric acid and related metabolites, such as γ-guanidinobutyric acid, and an age-dependent downregulation of cerebral γ-aminobutyric acid receptors. These findings indicate impaired γ-aminobutyric acid and γ-aminobutyric acid sub-type A (GABAA) receptor signalling as major factors underlying the pathophysiology of this neurometabolic disorder. We studied the cortical oscillation patterns and their relationship with γ-aminobutyric acid metabolism in 18 children affected by this condition and 10 healthy controls. Using high-density EEG, we recorded somatosensory cortical responses and resting-state activity. Using electrical source imaging, we estimated the relative power changes (compared with baseline) in both stimulus-evoked and stimulus-induced responses for physiologically relevant frequency bands and resting-state power. Stimulus-evoked oscillations are phase locked to the stimulus, whereas induced oscillations are not. Power changes for both evoked and induced responses as well as resting-state power were correlated with plasma γ-aminobutyric acid and γ-guanidinobutyric acid concentrations and with cortical γ-aminobutyric acid measured by proton magnetic resonance spectroscopy. Plasma γ-aminobutyric acid, γ-guanidinobutyric acid and cortical γ-aminobutyric acid were higher in patients than in controls (P < 0.001 for both). Beta and gamma relative power were suppressed for evoked responses in patients versus controls (P < 0.01). No group differences were observed for induced activity (P > 0.05). The mean gamma frequency of evoked responses was lower in patients versus controls (P = 0.002). Resting-state activity was suppressed in patients for theta (P = 0.011) and gamma (P < 0.001) bands. Evoked power changes were inversely correlated with plasma γ-aminobutyric acid and with γ-guanidinobutyric acid for beta (P < 0.001) and gamma (P < 0.001) bands. Similar relationships were observed between the evoked power changes and cortical γ-aminobutyric acid for all tested areas in the beta band (P < 0.001) and for the posterior cingulate gyrus in the gamma band (P < 0.001). We also observed a negative correlation between resting-state activity and plasma γ-aminobutyric acid and γ-guanidinobutyric acid for theta (P < 0.001; P = 0.003), alpha (P = 0.003; P = 0.02) and gamma (P = 0.02; P = 0.01) bands. Our findings indicate that increased γ-aminobutyric acid concentration is associated with reduced sensory-evoked beta and gamma activity and impaired neuronal synchronization in patients with succinic semialdehyde dehydrogenase deficiency. This further elucidates the pathophysiology of this neurometabolic disorder and serves as a potential biomarker for therapeutic trials.

Enhanced brain structure-function tethering in transmodal cortex revealed by high-frequency eigenmodes

While the link between brain structure and function remains an ongoing challenge, the prevailing hypothesis is that the structure-function relationship may itself be gradually decoupling from unimodal to transmodal cortex. However, this hypothesis is constrained by the underlying models which may neglect requisite information. Here we relate structural and functional connectivity derived from diffusion and functional MRI through orthogonal eigenmodes governing frequency-specific diffusion patterns. We find that low-frequency eigenmodes contribute little to functional interactions in transmodal cortex, resulting in divergent structure-function relationships. Conversely, high-frequency eigenmodes predominantly support neuronal coactivation patterns in these areas, inducing structure-function convergence along a unimodal-transmodal hierarchy. High-frequency information, although weak and scattered, could enhance the structure-function tethering, especially in transmodal association cortices. Our findings suggest that the structure-function decoupling may not be an intrinsic property of brain organization, but can be narrowed through multiplexed and regionally specialized spatiotemporal propagation regimes.

Memory rescue and learning in synaptic impaired neuronal circuits

Neuronal impairment is a characteristic of Alzheimer's disease (AD), but its effect on neural activity dynamics underlying memory deficits is unclear. Here, we studied the effects of synaptic impairment on neural activities associated with memory recall, memory rescue, and learning a new memory, in an integrate-and-fire neuronal network. Our results showed that reducing connectivity decreases the neuronal synchronization of memory neurons and impairs memory recall performance. Although, slow-gamma stimulation rescued memory recall and slow-gamma oscillations, the rescue caused a side effect of activating mixed memories. During the learning of a new memory, reducing connectivity caused impairment in storing the new memory, but did not affect previously stored memories. We also explored the effects of other types of impairments including neuronal loss and excitation-inhibition imbalance and the rescue by general increase of excitability. Our results reveal potential computational mechanisms underlying the memory deficits caused by impairment in AD.

The h-current controls cortical recurrent network activity through modulation of dendrosomatic communication

A fundamental feature of the cerebral cortex is the ability to rapidly turn on and off maintained activity within ensembles of neurons through recurrent excitation balanced by inhibition. Here we demonstrate that reduction of the h-current, which is especially prominent in pyramidal cell dendrites, strongly increases the ability of local cortical networks to generate maintained recurrent activity. Reduction of the h-current resulted in hyperpolarization and increase in input resistance of both the somata and apical dendrites of layer 5 pyramidal cells, while strongly increasing the dendrosomatic transfer of low (<20 Hz) frequencies, causing an increased responsiveness to dynamic clamp-induced recurrent network-like activity injected into the dendrites and substantially increasing the duration of spontaneous Up states. We propose that modulation of the h-current may strongly control the ability of cortical networks to generate recurrent persistent activity and the formation and dissolution of neuronal ensembles.

Different resting membrane potentials in posterior parietal cortex and prefrontal cortex in the view of recurrent synaptic strengths and neural network dynamics

In this study, we introduce the importance of elevated membrane potentials (MPs) in the prefrontal cortex (PFC) compared to that in the posterior parietal cortex (PPC), based on new observations of different MP levels in these areas. Through experimental data and spiking neural network modeling, we investigated a possible mechanism of the elevated membrane potential in the PFC and how these physiological differences affect neural network dynamics and cognitive functions in the PPC and PFC. Our findings indicate that NMDA receptors may be a main contributor to the elevated MP in the PFC region and highlight the potential of using a modeling toolkit to investigate the means by which changes in synaptic properties can affect neural dynamics and potentiate desirable cognitive functions through population activities in the corresponding brain regions.

Optogenetic manipulation of inhibitory interneurons can be used to validate a model of spatiotemporal sequence learning

The brain uses temporal information to link discrete events into memory structures supporting recognition, prediction, and a wide variety of complex behaviors. It is still an open question how experience-dependent synaptic plasticity creates memories including temporal and ordinal information. Various models have been proposed to explain how this could work, but these are often difficult to validate in a living brain. A recent model developed to explain sequence learning in the visual cortex encodes intervals in recurrent excitatory synapses and uses a learned offset between excitation and inhibition to generate precisely timed "messenger" cells that signal the end of an instance of time. This mechanism suggests that the recall of stored temporal intervals should be particularly sensitive to the activity of inhibitory interneurons that can be easily targeted in vivo with standard optogenetic tools. In this work we examined how simulated optogenetic manipulations of inhibitory cells modifies temporal learning and recall based on these mechanisms. We show that disinhibition and excess inhibition during learning or testing cause characteristic errors in recalled timing that could be used to validate the model in vivo using either physiological or behavioral measurements.

Dopaminergic Modulation of Prefrontal Cortex Inhibition

The prefrontal cortex is the highest stage of integration in the mammalian brain. Its functions vary greatly, from working memory to decision-making, and are primarily related to higher cognitive functions. This explains the considerable effort devoted to investigating this area, revealing the complex molecular, cellular, and network organization, and the essential role of various regulatory controls. In particular, the dopaminergic modulation and the impact of local interneurons activity are critical for prefrontal cortex functioning, controlling the excitatory/inhibitory balance and the overall network processing. Though often studied separately, the dopaminergic and GABAergic systems are deeply intertwined in influencing prefrontal network processing. This mini review will focus on the dopaminergic modulation of GABAergic inhibition, which plays a significant role in shaping prefrontal cortex activity.

Resolving the prefrontal mechanisms of adaptive cognitive behaviors: A cross-species perspective

The prefrontal cortex (PFC) enables a staggering variety of complex behaviors, such as planning actions, solving problems, and adapting to new situations according to external information and internal states. These higher-order abilities, collectively defined as adaptive cognitive behavior, require cellular ensembles that coordinate the tradeoff between the stability and flexibility of neural representations. While the mechanisms underlying the function of cellular ensembles are still unclear, recent experimental and theoretical studies suggest that temporal coordination dynamically binds prefrontal neurons into functional ensembles. A so far largely separate stream of research has investigated the prefrontal efferent and afferent connectivity. These two research streams have recently converged on the hypothesis that prefrontal connectivity patterns influence ensemble formation and the function of neurons within ensembles. Here, we propose a unitary concept that, leveraging a cross-species definition of prefrontal regions, explains how prefrontal ensembles adaptively regulate and efficiently coordinate multiple processes in distinct cognitive behaviors.

Using Nonhuman Primate Models to Reverse-Engineer Prefrontal Circuit Failure Underlying Cognitive Deficits in Schizophrenia

In this chapter, I review studies in nonhuman primates that emulate the circuit failure in prefrontal cortex responsible for working memory and cognitive control deficits in schizophrenia. These studies have characterized how synaptic malfunction, typically induced by blockade of NMDAR, disrupts neural function and computation in prefrontal networks to explain errors in cognitive tasks that are seen in schizophrenia. This work is finding causal relationships between pathogenic events of relevance to schizophrenia at vastly different levels of scale, from synapses, to neurons, local, circuits, distributed networks, computation, and behavior. Pharmacological manipulation, the dominant approach in primate models, has limited construct validity for schizophrenia pathogenesis, as the disease results from a complex interplay between environmental, developmental, and genetic factors. Genetic manipulation replicating schizophrenia risk is more advanced in rodent models. Nonetheless, gene manipulation in nonhuman primates is rapidly advancing, and primate developmental models have been established. Integration of large scale neural recording, genetic manipulation, and computational modeling in nonhuman primates holds considerable potential to provide a crucial schizophrenia model moving forward. Data generated by this approach is likely to fill several crucial gaps in our understanding of the causal sequence leading to schizophrenia in humans. This causal chain presents a vexing problem largely because it requires understanding how events at very different levels of scale relate to one another, from genes to circuits to cognition to social interactions. Nonhuman primate models excel here. They optimally enable discovery of causal relationships across levels of scale in the brain that are relevant to cognitive deficits in schizophrenia. The mechanistic understanding of prefrontal circuit failure they promise to provide may point the way to more effective therapeutic interventions to restore function to prefrontal networks in the disease.

A multi-faceted role of dual-state dopamine signaling in working memory, attentional control, and intelligence

Genetic evidence strongly suggests that individual differences in intelligence will not be reducible to a single dominant cause. However, some of those variations/changes may be traced to tractable, cohesive mechanisms. One such mechanism may be the balance of dopamine D1 (D₁R) and D2 (D₂R) receptors, which regulate intrinsic currents and synaptic transmission in frontal cortical regions. Here, we review evidence from human, animal, and computational studies that suggest that this balance (in density, activity state, and/or availability) is critical to the implementation of executive functions such as attention and working memory, both of which are principal contributors to variations in intelligence. D1 receptors dominate neural responding during stable periods of short-term memory maintenance (requiring attentional focus), while D2 receptors play a more specific role during periods of instability such as changing environmental or memory states (requiring attentional disengagement). Here we bridge these observations with known properties of human intelligence. Starting from theories of intelligence that place executive functions (e.g., working memory and attentional control) at its center, we propose that dual-state dopamine signaling might be a causal contributor to at least some of the variation in intelligence across individuals and its change by experiences/training. Although it is unlikely that such a mechanism can account for more than a modest portion of the total variance in intelligence, our proposal is consistent with an array of available evidence and has a high degree of explanatory value. We suggest future directions and specific empirical tests that can further elucidate these relationships.

Assessing the coupling between local neural activity and global connectivity fluctuations: Application to human intracranial electroencephalography during a cognitive task

Cognitive-relevant information is processed by different brain areas that cooperate to eventually produce a response. The relationship between local activity and global brain states during such processes, however, remains for the most part unexplored. To address this question, we designed a simple face-recognition task performed in patients with drug-resistant epilepsy and monitored with intracranial electroencephalography (EEG). Based on our observations, we developed a novel analytical framework (named "local-global" framework) to statistically correlate the brain activity in every recorded gray-matter region with the widespread connectivity fluctuations as proxy to identify concurrent local activations and global brain phenomena that may plausibly reflect a common functional network during cognition. The application of the local-global framework to the data from three subjects showed that similar connectivity fluctuations found across patients were mainly coupled to the local activity of brain areas involved in face information processing. In particular, our findings provide preliminary evidence that the reported global measures might be a novel signature of functional brain activity reorganization when a stimulus is processed in a task context regardless of the specific recorded areas.

Multiple forms of working memory emerge from synapse-astrocyte interactions in a neuron-glia network model

Persistent activity in populations of neurons, time-varying activity across a neural population, or activity-silent mechanisms carried out by hidden internal states of the neural population have been proposed as different mechanisms of working memory (WM). Whether these mechanisms could be mutually exclusive or occur in the same neuronal circuit remains, however, elusive, and so do their biophysical underpinnings. While WM is traditionally regarded to depend purely on neuronal mechanisms, cortical networks also include astrocytes that can modulate neural activity. We propose and investigate a network model that includes both neurons and glia and show that glia-synapse interactions can lead to multiple stable states of synaptic transmission. Depending on parameters, these interactions can lead in turn to distinct patterns of network activity that can serve as substrates for WM.

Geometry of neural computation unifies working memory and planning

Real-world tasks require coordination of working memory, decision-making, and planning, yet these cognitive functions have disproportionately been studied as independent modular processes in the brain. Here, we propose that contingency representations, defined as mappings for how future behaviors depend on upcoming events, can unify working memory and planning computations. We designed a task capable of disambiguating distinct types of representations. In task-optimized recurrent neural networks, we investigated possible circuit mechanisms for contingency representations and found that these representations can explain neurophysiological observations from the prefrontal cortex during working memory tasks. Our experiments revealed that human behavior is consistent with contingency representations and not with traditional sensory models of working memory. Finally, we generated falsifiable predictions for neural data to identify contingency representations in neural data and to dissociate different models of working memory. Our findings characterize a neural representational strategy that can unify working memory, planning, and context-dependent decision-making.

Relating a Spiking Neural Network Model and the Diffusion Model of Decision-Making

Many models of decision making assume accumulation of evidence to threshold as a core mechanism to predict response probabilities and response times. A spiking neural network model (Wang, 2002) instantiates these mechanisms at the level of biophysically-plausible pools of neurons with excitatory and inhibitory connections, and has numerous model parameters tuned by physiological measures. The diffusion model (Ratcliff, 1978) is a cognitive model that can be fitted to a range of behaviors and conditions. We investigated how parameters of the cognitive-level diffusion model relate to the parameters of a neural-level spiking model. In each simulated "experiment", we generated "data" from the spiking neural network by factorially combining a manipulation of choice difficulty (via the input to the spiking model) and a manipulation of one of the core parameters of the spiking model. We then fitted the diffusion model to these simulated data to observe how manipulation of each core spiking model parameter mapped on to fitted drift rate, response threshold, and non-decision time. Manipulations of parameters in the spiking model related to input sensitivity, threshold, and stimulus processing time mapped on to their conceptual analogues in the diffusion model, namely drift rate, threshold, and non-decision time. Manipulations of parameters in the spiking model with no direct analogue to the diffusion model, non-stimulus-specific background input, strength of recurrent excitation, and receptor conductances, mapped on to threshold in the diffusion model. We discuss implications of these results for interpretations of fits of the diffusion model to behavioral data.

Dynamic branching in a neural network model for probabilistic prediction of sequences

An important function of the brain is to predict which stimulus is likely to occur based on the perceived cues. The present research studied the branching behavior of a computational network model of populations of excitatory and inhibitory neurons, both analytically and through simulations. Results show how synaptic efficacy, retroactive inhibition and short-term synaptic depression determine the dynamics of selection between different branches predicting sequences of stimuli of different probabilities. Further results show that changes in the probability of the different predictions depend on variations of neuronal gain. Such variations allow the network to optimize the probability of its predictions to changing probabilities of the sequences without changing synaptic efficacy.

Reproducing a decision-making network in a virtual visual discrimination task

We reproduced a decision-making network model using the neural simulator software neural simulation tool (NEST), and we embedded the spiking neural network in a virtual robotic agent performing a simulated behavioral task. The present work builds upon the concept of replicability in neuroscience, preserving most of the computational properties in the initial model although employing a different software tool. The proposed implementation successfully obtains equivalent results from the original study, reproducing the salient features of the neural processes underlying a binary decision. Furthermore, the resulting network is able to control a robot performing an in silico visual discrimination task, the implementation of which is openly available on the EBRAINS infrastructure through the neuro robotics platform (NRP).

Theory of the Multiregional Neocortex: Large-Scale Neural Dynamics and Distributed Cognition

The neocortex is a complex neurobiological system with many interacting regions. How these regions work together to subserve flexible behavior and cognition has become increasingly amenable to rigorous research. Here, I review recent experimental and theoretical work on the modus operandi of a multiregional cortex. These studies revealed several general principles for the neocortical interareal connectivity, low-dimensional macroscopic gradients of biological properties across cortical areas, and a hierarchy of timescales for information processing. Theoretical work suggests testable predictions regarding differential excitation and inhibition along feedforward and feedback pathways in the cortical hierarchy. Furthermore, modeling of distributed working memory and simple decision-making has given rise to a novel mathematical concept, dubbed bifurcation in space, that potentially explains how different cortical areas, with a canonical circuit organization but gradients of biological heterogeneities, are able to subserve their respective (e.g., sensory coding versus executive control) functions in a modularly organized brain.

Inhibitory control of frontal metastability sets the temporal signature of cognition

Cortical dynamics are organized over multiple anatomical and temporal scales. The mechanistic origin of the temporal organization and its contribution to cognition remain unknown. Here we demonstrate the cause of this organization by studying a specific temporal signature (time constant and latency) of neural activity. In monkey frontal areas, recorded during flexible decisions, temporal signatures display specific area-dependent ranges, as well as anatomical and cell-type distributions. Moreover, temporal signatures are functionally adapted to behaviorally relevant timescales. Fine-grained biophysical network models, constrained to account for experimentally observed temporal signatures, reveal that after-hyperpolarization potassium and inhibitory GABA-B conductances critically determine areas' specificity. They mechanistically account for temporal signatures by organizing activity into metastable states, with inhibition controlling state stability and transitions. As predicted by models, state durations non-linearly scale with temporal signatures in monkey, matching behavioral timescales. Thus, local inhibitory-controlled metastability constitutes the dynamical core specifying the temporal organization of cognitive functions in frontal areas.

The geometry of domain-general performance monitoring in the human medial frontal cortex

Controlling behavior to flexibly achieve desired goals depends on the ability to monitor one's own performance. It is unknown how performance monitoring can be both flexible, to support different tasks, and specialized, to perform each task well. We recorded single neurons in the human medial frontal cortex while subjects performed two tasks that involve three types of cognitive conflict. Neurons encoding conflict probability, conflict, and error in one or both tasks were intermixed, forming a representational geometry that simultaneously allowed task specialization and generalization. Neurons encoding conflict retrospectively served to update internal estimates of conflict probability. Population representations of conflict were compositional. These findings reveal how representations of evaluative signals can be both abstract and task-specific and suggest a neuronal mechanism for estimating control demand.

Neuromodulation and Neurophysiology on the Timescale of Learning and Decision-Making

Nervous systems evolved to effectively navigate the dynamics of the environment to achieve their goals. One framework used to study this fundamental problem arose in the study of learning and decision-making. In this framework, the demands of effective behavior require slow dynamics-on the scale of seconds to minutes-of networks of neurons. Here, we review the phenomena and mechanisms involved. Using vignettes from a few species and areas of the nervous system, we view neuromodulators as key substrates for temporal scaling of neuronal dynamics. Expected final online publication date for the Annual Review of Neuroscience, Volume 45 is July 2022. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

Effects of virtual lesions on temporal dynamics in cortical networks based on personalized dynamic models

The human brain dynamically shifts between a predominantly integrated state and a predominantly segregated state, each with different roles in supporting cognition and behavior. However, no studies to date have investigated lesions placed in different regions of the cerebral cortex to determine the effects on the temporal balance of integration and segregation. Here, we used the integrated state occurrence rate to measure the temporal balance of integration and segregation in the resting brain. Based on dynamic mean-field models coupled through the individual's structural white matter connections, neural activity was modeled. By lesioning different individual nodes from the model, our results implied that the impact of lesions reaches far beyond focal damage and could impair cognition by affecting system-level dynamics. Lesions in a portion of the nodes in the visual, somatomotor, limbic, and default networks significantly compromised the temporal balance of integration and segregation. Crucially, the effects of lesions in different regions could be predicted based on the hierarchical axis inferred from the T1w/T2w map and specific graph measures of structural brain networks. Taken together, our findings indicate the possibility of significant contributions of anatomical heterogeneity to the dynamics of functional network topology. Although still in its infancy, computational modeling may provide an entry point for understanding brain disorders at a causal mechanistic level, possibly leading to novel, more effective therapeutic interventions.

Metastable dynamics of neural circuits and networks

Cortical neurons emit seemingly erratic trains of action potentials or "spikes," and neural network dynamics emerge from the coordinated spiking activity within neural circuits. These rich dynamics manifest themselves in a variety of patterns, which emerge spontaneously or in response to incoming activity produced by sensory inputs. In this Review, we focus on neural dynamics that is best understood as a sequence of repeated activations of a number of discrete hidden states. These transiently occupied states are termed "metastable" and have been linked to important sensory and cognitive functions. In the rodent gustatory cortex, for instance, metastable dynamics have been associated with stimulus coding, with states of expectation, and with decision making. In frontal, parietal, and motor areas of macaques, metastable activity has been related to behavioral performance, choice behavior, task difficulty, and attention. In this article, we review the experimental evidence for neural metastable dynamics together with theoretical approaches to the study of metastable activity in neural circuits. These approaches include (i) a theoretical framework based on non-equilibrium statistical physics for network dynamics; (ii) statistical approaches to extract information about metastable states from a variety of neural signals; and (iii) recent neural network approaches, informed by experimental results, to model the emergence of metastable dynamics. By discussing these topics, we aim to provide a cohesive view of how transitions between different states of activity may provide the neural underpinnings for essential functions such as perception, memory, expectation, or decision making, and more generally, how the study of metastable neural activity may advance our understanding of neural circuit function in health and disease.

Mechanisms of distributed working memory in a large-scale network of macaque neocortex

Neural activity underlying working memory is not a local phenomenon but distributed across multiple brain regions. To elucidate the circuit mechanism of such distributed activity, we developed an anatomically constrained computational model of large-scale macaque cortex. We found that mnemonic internal states may emerge from inter-areal reverberation, even in a regime where none of the isolated areas is capable of generating self-sustained activity. The mnemonic activity pattern along the cortical hierarchy indicates a transition in space, separating areas engaged in working memory and those which do not. A host of spatially distinct attractor states is found, potentially subserving various internal processes. The model yields testable predictions, including the idea of counterstream inhibitory bias, the role of prefrontal areas in controlling distributed attractors, and the resilience of distributed activity to lesions or inactivation. This work provides a theoretical framework for identifying large-scale brain mechanisms and computational principles of distributed cognitive processes.

Connectivity and dynamics in the olfactory bulb

Dendrodendritic interactions between excitatory mitral cells and inhibitory granule cells in the olfactory bulb create a dense interaction network, reorganizing sensory representations of odors and, consequently, perception. Large-scale computational models are needed for revealing how the collective behavior of this network emerges from its global architecture. We propose an approach where we summarize anatomical information through dendritic geometry and density distributions which we use to calculate the connection probability between mitral and granule cells, while capturing activity patterns of each cell type in the neural dynamical systems theory of Izhikevich. In this way, we generate an efficient, anatomically and physiologically realistic large-scale model of the olfactory bulb network. Our model reproduces known connectivity between sister vs. non-sister mitral cells; measured patterns of lateral inhibition; and theta, beta, and gamma oscillations. The model in turn predicts testable relationships between network structure and several functional properties, including lateral inhibition, odor pattern decorrelation, and LFP oscillation frequency. We use the model to explore the influence of cortex on the olfactory bulb, demonstrating possible mechanisms by which cortical feedback to mitral cells or granule cells can influence bulbar activity, as well as how neurogenesis can improve bulbar decorrelation without requiring cell death. Our methodology provides a tractable tool for other researchers.

Criticality enhances the multilevel reliability of stimulus responses in cortical neural networks

Cortical neural networks exhibit high internal variability in spontaneous dynamic activities and they can robustly and reliably respond to external stimuli with multilevel features–from microscopic irregular spiking of neurons to macroscopic oscillatory local field potential. A comprehensive study integrating these multilevel features in spontaneous and stimulus–evoked dynamics with seemingly distinct mechanisms is still lacking. Here, we study the stimulus–response dynamics of biologically plausible excitation–inhibition (E–I) balanced networks. We confirm that networks around critical synchronous transition states can maintain strong internal variability but are sensitive to external stimuli. In this dynamical region, applying a stimulus to the network can reduce the trial-to-trial variability and shift the network oscillatory frequency while preserving the dynamical criticality. These multilevel features widely observed in different experiments cannot simultaneously occur in non-critical dynamical states. Furthermore, the dynamical mechanisms underlying these multilevel features are revealed using a semi-analytical mean-field theory that derives the macroscopic network field equations from the microscopic neuronal networks, enabling the analysis by nonlinear dynamics theory and linear noise approximation. The generic dynamical principle revealed here contributes to a more integrative understanding of neural systems and brain functions and incorporates multimodal and multilevel experimental observations. The E–I balanced neural network in combination with the effective mean-field theory can serve as a mechanistic modeling framework to study the multilevel neural dynamics underlying neural information and cognitive processes.

The complexity and variability of brain dynamical activity range from neuronal spiking and neural avalanches to oscillatory local field potentials of local neural circuits in both spontaneous and stimulus-evoked states. Such multilevel variable brain dynamics are functionally and behaviorally relevant and are principal components of the underlying circuit organization. To more comprehensively clarify their neural mechanisms, we use a bottom-up approach to study the stimulus–response dynamics of neural circuits. Our model assumes the following key biologically plausible components: excitation–inhibition (E–I) neuronal interaction and chemical synaptic coupling. We show that the circuits with E–I balance have a special dynamic sub-region, the critical region. Circuits around this region could account for the emergence of multilevel brain response patterns, both ongoing and stimulus-induced, observed in different experiments, including the reduction of trial-to-trial variability, effective modulation of gamma frequency, and preservation of criticality in the presence of a stimulus. We further analyze the corresponding nonlinear dynamical principles using a novel and highly generalizable semi-analytical mean-field theory. Our computational and theoretical studies explain the cross-level brain dynamical organization of spontaneous and evoked states in a more integrative manner.

Dynamic primitives of brain network interaction Special Issue "Advances in Mapping the Connectome"

What dynamic processes underly functional brain networks? Functional connectivity (FC) and functional connectivity dynamics (FCD) are used to represent the patterns and dynamics of functional brain networks. FC(D) is related to the synchrony of brain activity: when brain areas oscillate in a coordinated manner this yields a high correlation between their signal time series. To explain the processes underlying FC(D) we review how synchronized oscillations emerge from coupled neural populations in brain network models (BNMs). From detailed spiking networks to more abstract population models, there is strong support for the idea that the brain operates near critical instabilities that give rise to multistable or metastable dynamics that in turn lead to the intermittently synchronized slow oscillations underlying FC(D). We explore further consequences from these fundamental mechanisms and how they fit with reality. We conclude by highlighting the need for integrative brain models that connect separate mechanisms across levels of description and spatiotemporal scales and link them with cognitive function.

Dopamine firing plays a dual role in coding reward prediction errors and signaling motivation in a working memory task

Little is known about how dopamine (DA) neuron firing rates behave in cognitively demanding decision-making tasks. Here, we investigated midbrain DA activity in monkeys performing a discrimination task in which the animal had to use working memory (WM) to report which of two sequentially applied vibrotactile stimuli had the higher frequency. We found that perception was altered by an internal bias, likely generated by deterioration of the representation of the first frequency during the WM period. This bias greatly controlled the DA phasic response during the two stimulation periods, confirming that DA reward prediction errors reflected stimulus perception. In contrast, tonic dopamine activity during WM was not affected by the bias and did not encode the stored frequency. More interestingly, both delay-period activity and phasic responses before the second stimulus negatively correlated with reaction times of the animals after the trial start cue and thus represented motivated behavior on a trial-by-trial basis. During WM, this motivation signal underwent a ramp-like increase. At the same time, motivation positively correlated with accuracy, especially in difficult trials, probably by decreasing the effect of the bias. Overall, our results indicate that DA activity, in addition to encoding reward prediction errors, could at the same time be involved in motivation and WM. In particular, the ramping activity during the delay period suggests a possible DA role in stabilizing sustained cortical activity, hypothetically by increasing the gain communicated to prefrontal neurons in a motivation-dependent way.

The Mean Field Approach for Populations of Spiking Neurons

Mean field theory is a device to analyze the collective behavior of a dynamical system comprising many interacting particles. The theory allows to reduce the behavior of the system to the properties of a handful of parameters. In neural circuits, these parameters are typically the firing rates of distinct, homogeneous subgroups of neurons. Knowledge of the firing rates under conditions of interest can reveal essential information on both the dynamics of neural circuits and the way they can subserve brain function. The goal of this chapter is to provide an elementary introduction to the mean field approach for populations of spiking neurons. We introduce the general idea in networks of binary neurons, starting from the most basic results and then generalizing to more relevant situations. This allows to derive the mean field equations in a simplified setting. We then derive the mean field equations for populations of integrate-and-fire neurons. An effort is made to derive the main equations of the theory using only elementary methods from calculus and probability theory. The chapter ends with a discussion of the assumptions of the theory and some of the consequences of violating those assumptions. This discussion includes an introduction to balanced and metastable networks and a brief catalogue of successful applications of the mean field approach to the study of neural circuits.

Emergence of Irregular Activity in Networks of Strongly Coupled Conductance-Based Neurons

Cortical neurons are characterized by irregular firing and a broad distribution of rates. The balanced state model explains these observations with a cancellation of mean excitatory and inhibitory currents, which makes fluctuations drive firing. In networks of neurons with current-based synapses, the balanced state emerges dynamically if coupling is strong, i.e., if the mean number of synapses per neuron K is large and synaptic efficacy is of the order of 1 / K . When synapses are conductance-based, current fluctuations are suppressed when coupling is strong, questioning the applicability of the balanced state idea to biological neural networks. We analyze networks of strongly coupled conductance-based neurons and show that asynchronous irregular activity and broad distributions of rates emerge if synaptic efficacy is of the order of 1/ log(K). In such networks, unlike in the standard balanced state model, current fluctuations are small and firing is maintained by a drift-diffusion balance. This balance emerges dynamically, without fine-tuning, if inputs are smaller than a critical value, which depends on synaptic time constants and coupling strength, and is significantly more robust to connection heterogeneities than the classical balanced state model. Our analysis makes experimentally testable predictions of how the network response properties should evolve as input increases.

A data-informed mean-field approach to mapping of cortical parameter landscapes

Constraining the many biological parameters that govern cortical dynamics is computationally and conceptually difficult because of the curse of dimensionality. This paper addresses these challenges by proposing (1) a novel data-informed mean-field (MF) approach to efficiently map the parameter space of network models; and (2) an organizing principle for studying parameter space that enables the extraction biologically meaningful relations from this high-dimensional data. We illustrate these ideas using a large-scale network model of the Macaque primary visual cortex. Of the 10-20 model parameters, we identify 7 that are especially poorly constrained, and use the MF algorithm in (1) to discover the firing rate contours in this 7D parameter cube. Defining a "biologically plausible" region to consist of parameters that exhibit spontaneous Excitatory and Inhibitory firing rates compatible with experimental values, we find that this region is a slightly thickened codimension-1 submanifold. An implication of this finding is that while plausible regimes depend sensitively on parameters, they are also robust and flexible provided one compensates appropriately when parameters are varied. Our organizing principle for conceptualizing parameter dependence is to focus on certain 2D parameter planes that govern lateral inhibition: Intersecting these planes with the biologically plausible region leads to very simple geometric structures which, when suitably scaled, have a universal character independent of where the intersections are taken. In addition to elucidating the geometry of the plausible region, this invariance suggests useful approximate scaling relations. Our study offers, for the first time, a complete characterization of the set of all biologically plausible parameters for a detailed cortical model, which has been out of reach due to the high dimensionality of parameter space.

Attention improves information flow between neuronal populations without changing the communication subspace

Visual attention allows observers to change the influence of different parts of a visual scene on their behavior, suggesting that information can be flexibly shared between visual cortex and neurons involved in decision making. We investigated the neural substrate of flexible information routing by analyzing the activity of populations of visual neurons in the medial temporal area (MT) and oculo-motor neurons in the superior colliculus (SC) while rhesus monkeys switched spatial attention. We demonstrated that attention increases the efficacy of visuomotor communication: trial-to-trial variability in SC population activity could be better predicted by the activity of the MT population (and vice versa) when attention was directed toward their joint receptive fields. Surprisingly, this improvement in prediction was not explained by changes in the dimensionality of the shared subspace or in the magnitude of local or shared pairwise noise correlations. These results lay a foundation for future theoretical and experimental studies into how visual attention can flexibly change information flow between sensory and decision neurons.

Towards biologically constrained attractor models of schizophrenia

Alterations in neuromodulation or synaptic transmission in biophysical attractor network models, as proposed by the dominant dopaminergic and glutamatergic theories of schizophrenia, successfully mimic working memory (WM) deficits in people with schizophrenia (PSZ). Yet, multiple, often opposing alterations in memory circuits can lead to the same behavioral patterns in these network models. Here, we critically revise the computational and experimental literature that links NMDAR hypofunction to WM precision loss in PSZ. We show in network simulations that currently available experimental evidence cannot set apart competing biophysical accounts. Critical points to resolve are the effects of increases vs. decreases in E/I ratio (e.g. through NMDAR blockade) on firing rate tuning and shared noise modulations and possible concomitant deficits in short-term plasticity. We argue that these concerted experimental and computational efforts will lead to a better understanding of the neurobiology underlying cognitive deficits in PSZ.

Memory retrieval dynamics and storage capacity of a modular network model of association cortex with featural decomposition

The primate heteromodal cortex presents an evident functional modularity at a mesoscopic level, with physiological and anatomical evidence pointing to it as likely substrate of long-term memory. In order to investigate some of its properties, a model of multimodular autoassociator is studied. Each of the many modules represents a neocortical functional ensemble of recurrently connected neurons and operates as a Hebbian autoassociator, storing a number of local features which it can recall upon cue. The global memory patterns are made of combinations of features sparsely distributed across the modules. Intermodular connections are modelled as a finite-connectivity random graph. Any pair of features in any respective pair of modules is allowed to be involved in several memory patterns; the coarse-grained modular network dynamics is defined in such a way as to overcome the consequent ambiguity of associations. Effects of long-range homeostatic synaptic scaling on network performance are also assessed. The dynamical process of cued retrieval almost saturates a natural upper bound while producing negligible spurious activation. The extent of finite-size effects on storage capacity is quantitatively evaluated. In the limit of infinite size, the functional relationship between storage capacity and number of features per module reduces to that which other authors found by methods from equilibrium statistical mechanics, which suggests that the origin of the functional form is of a combinatorial nature. In contrast with its apparent inevitability at intramodular level, long-range synaptic scaling results to be of minor relevance to both retrieval and storage capacity, casting doubt on its existence in the neocortex. A conjecture is also posited about how statistical fluctuation of connectivity across the network may underpin spontaneous emergence of semantic hierarchies through learning.