A mechanism for deviance detection and contextual routing in the thalamus: a review and theoretical proposal

Predictive processing theories conceptualize neocortical feedback as conveying expectations and contextual attention signals derived from internal cortical models, playing an essential role in the perception and interpretation of sensory information. However, few predictive processing frameworks outline concrete mechanistic roles for the corticothalamic (CT) feedback from layer 6 (L6), despite the fact that the number of CT axons is an order of magnitude greater than that of feedforward thalamocortical (TC) axons. Here we review the functional architecture of CT circuits and propose a mechanism through which L6 could regulate thalamic firing modes (burst, tonic) to detect unexpected inputs. Using simulations in a model of a TC cell, we show how the CT feedback could support prediction-based input discrimination in TC cells by promoting burst firing. This type of CT control can enable the thalamic circuit to implement spatial and context selective attention mechanisms. The proposed mechanism generates specific experimentally testable hypotheses. We suggest that the L6 CT feedback allows the thalamus to detect deviance from predictions of internal cortical models, thereby supporting contextual attention and routing operations, a far more powerful role than traditionally assumed.

Data-driven multiscale model of macaque auditory thalamocortical circuits reproduces in vivo dynamics

We developed a detailed model of macaque auditory thalamocortical circuits, including primary auditory cortex (A1), medial geniculate body (MGB), and thalamic reticular nucleus, utilizing the NEURON simulator and NetPyNE tool. The A1 model simulates a cortical column with over 12,000 neurons and 25 million synapses, incorporating data on cell-type-specific neuron densities, morphology, and connectivity across six cortical layers. It is reciprocally connected to the MGB thalamus, which includes interneurons and core and matrix-layer-specific projections to A1. The model simulates multiscale measures, including physiological firing rates, local field potentials (LFPs), current source densities (CSDs), and electroencephalography (EEG) signals. Laminar CSD patterns, during spontaneous activity and in response to broadband noise stimulus trains, mirror experimental findings. Physiological oscillations emerge spontaneously across frequency bands comparable to those recorded in vivo. We elucidate population-specific contributions to observed oscillation events and relate them to firing and presynaptic input patterns. The model offers a quantitative theoretical framework to integrate and interpret experimental data and predict its underlying cellular and circuit mechanisms.

The Roles of Potassium and Calcium Currents in the Bistable Firing Transition

Healthy brains display a wide range of firing patterns, from synchronized oscillations during slow-wave sleep to desynchronized firing during movement. These physiological activities coexist with periods of pathological hyperactivity in the epileptic brain, where neurons can fire in synchronized bursts. Most cortical neurons are pyramidal regular spiking (RS) cells with frequency adaptation and do not exhibit bursts in current-clamp experiments (in vitro). In this work, we investigate the transition mechanism of spike-to-burst patterns due to slow potassium and calcium currents, considering a conductance-based model of a cortical RS cell. The joint influence of potassium and calcium ion channels on high synchronous patterns is investigated for different synaptic couplings (gsyn) and external current inputs (I). Our results suggest that slow potassium currents play an important role in the emergence of high-synchronous activities, as well as in the spike-to-burst firing pattern transitions. This transition is related to the bistable dynamics of the neuronal network, where physiological asynchronous states coexist with pathological burst synchronization. The hysteresis curve of the coefficient of variation of the inter-spike interval demonstrates that a burst can be initiated by firing states with neuronal synchronization. Furthermore, we notice that high-threshold (IL) and low-threshold (IT) ion channels play a role in increasing and decreasing the parameter conditions (gsyn and I) in which bistable dynamics occur, respectively. For high values of IL conductance, a synchronous burst appears when neurons are weakly coupled and receive more external input. On the other hand, when the conductance IT increases, higher coupling and lower I are necessary to produce burst synchronization. In light of our results, we suggest that channel subtype-specific pharmacological interactions can be useful to induce transitions from pathological high bursting states to healthy states.

The Role of Potassium and Calcium Currents in the Bistable Firing Transition

Healthy brains display a wide range of firing patterns, from synchronized oscillations during slowwave sleep to desynchronized firing during movement. These physiological activities coexist with periods of pathological hyperactivity in the epileptic brain, where neurons can fire in synchronized bursts. Most cortical neurons are pyramidal regular spiking cells (RS) with frequency adaptation and do not exhibit bursts in current-clamp experiments ( in vitro ). In this work, we investigate the transition mechanism of spike-to-burst patterns due to slow potassium and calcium currents, considering a conductance-based model of a cortical RS cell. The joint influence of potassium and calcium ion channels on high synchronous patterns is investigated for different synaptic couplings ( g syn ) and external current inputs ( I ). Our results suggest that slow potassium currents play an important role in the emergence of high-synchronous activities, as well as in the spike-to-burst firing pattern transitions. This transition is related to bistable dynamics of the neuronal network, where physiological asynchronous states coexist with pathological burst synchronization. The hysteresis curve of the coefficient of variation of the inter-spike interval demonstrates that a burst can be initiated by firing states with neuronal synchronization. Furthermore, we notice that high-threshold ( I L ) and low-threshold ( I T ) ion channels play a role in increasing and decreasing the parameter conditions ( g syn and I ) in which bistable dynamics occur, respectively. For high values of I L conductance, a synchronous burst appears when neurons are weakly coupled and receive more external input. On the other hand, when the conductance I T increases, higher coupling and lower I are necessary to produce burst synchronization. In light of our results, we suggest that channel subtype-specific pharmacological interactions can be useful to induce transitions from pathological high bursting states to healthy states.

Multiscale model of primary motor cortex circuits predicts in vivo cell-type-specific, behavioral state-dependent dynamics

Understanding cortical function requires studying multiple scales: molecular, cellular, circuit, and behavioral. We develop a multiscale, biophysically detailed model of mouse primary motor cortex (M1) with over 10,000 neurons and 30 million synapses. Neuron types, densities, spatial distributions, morphologies, biophysics, connectivity, and dendritic synapse locations are constrained by experimental data. The model includes long-range inputs from seven thalamic and cortical regions and noradrenergic inputs. Connectivity depends on cell class and cortical depth at sublaminar resolution. The model accurately predicts in vivo layer- and cell-type-specific responses (firing rates and LFP) associated with behavioral states (quiet wakefulness and movement) and experimental manipulations (noradrenaline receptor blockade and thalamus inactivation). We generate mechanistic hypotheses underlying the observed activity and analyzed low-dimensional population latent dynamics. This quantitative theoretical framework can be used to integrate and interpret M1 experimental data and sheds light on the cell-type-specific multiscale dynamics associated with several experimental conditions and behaviors.

Theta-gamma phase amplitude coupling in a hippocampal CA1 microcircuit

Phase amplitude coupling (PAC) between slow and fast oscillations is found throughout the brain and plays important functional roles. Its neural origin remains unclear. Experimental findings are often puzzling and sometimes contradictory. Most computational models rely on pairs of pacemaker neurons or neural populations tuned at different frequencies to produce PAC. Here, using a data-driven model of a hippocampal microcircuit, we demonstrate that PAC can naturally emerge from a single feedback mechanism involving an inhibitory and excitatory neuron population, which interplay to generate theta frequency periodic bursts of higher frequency gamma. The model suggests the conditions under which a CA1 microcircuit can operate to elicit theta-gamma PAC, and highlights the modulatory role of OLM and PVBC cells, recurrent connectivity, and short term synaptic plasticity. Surprisingly, the results suggest the experimentally testable prediction that the generation of the slow population oscillation requires the fast one and cannot occur without it.

Large-scale biophysically detailed model of somatosensory thalamocortical circuits in NetPyNE

The primary somatosensory cortex (S1) of mammals is critically important in the perception of touch and related sensorimotor behaviors. In 2015, the Blue Brain Project (BBP) developed a groundbreaking rat S1 microcircuit simulation with over 31,000 neurons with 207 morpho-electrical neuron types, and 37 million synapses, incorporating anatomical and physiological information from a wide range of experimental studies. We have implemented this highly detailed and complex S1 model in NetPyNE, using the data available in the Neocortical Microcircuit Collaboration Portal. NetPyNE provides a Python high-level interface to NEURON and allows defining complicated multiscale models using an intuitive declarative standardized language. It also facilitates running parallel simulations, automates the optimization and exploration of parameters using supercomputers, and provides a wide range of built-in analysis functions. This will make the S1 model more accessible and simpler to scale, modify and extend in order to explore research questions or interconnect to other existing models. Despite some implementation differences, the NetPyNE model preserved the original cell morphologies, electrophysiological responses and spatial distribution for all 207 cell types; and the connectivity properties of all 1941 pathways, including synaptic dynamics and short-term plasticity (STP). The NetPyNE S1 simulations produced reasonable physiological firing rates and activity patterns across all populations. When STP was included, the network generated a 1 Hz oscillation comparable to the original model in vitro-like state. By then reducing the extracellular calcium concentration, the model reproduced the original S1 in vivo-like states with asynchronous activity. These results validate the original study using a new modeling tool. Simulated local field potentials (LFPs) exhibited realistic oscillatory patterns and features, including distance- and frequency-dependent attenuation. The model was extended by adding thalamic circuits, including 6 distinct thalamic populations with intrathalamic, thalamocortical (TC) and corticothalamic connectivity derived from experimental data. The thalamic model reproduced single known cell and circuit-level dynamics, including burst and tonic firing modes and oscillatory patterns, providing a more realistic input to cortex and enabling study of TC interactions. Overall, our work provides a widely accessible, data-driven and biophysically-detailed model of the somatosensory TC circuits that can be employed as a community tool for researchers to study neural dynamics, function and disease.

Modernizing the NEURON Simulator for Sustainability, Portability, and Performance

The need for reproducible, credible, multiscale biological modeling has led to the development of standardized simulation platforms, such as the widely-used NEURON environment for computational neuroscience. Developing and maintaining NEURON over several decades has required attention to the competing needs of backwards compatibility, evolving computer architectures, the addition of new scales and physical processes, accessibility to new users, and efficiency and flexibility for specialists. In order to meet these challenges, we have now substantially modernized NEURON, providing continuous integration, an improved build system and release workflow, and better documentation. With the help of a new source-to-source compiler of the NMODL domain-specific language we have enhanced NEURON's ability to run efficiently, via the CoreNEURON simulation engine, on a variety of hardware platforms, including GPUs. Through the implementation of an optimized in-memory transfer mechanism this performance optimized backend is made easily accessible to users, providing training and model-development paths from laptop to workstation to supercomputer and cloud platform. Similarly, we have been able to accelerate NEURON's reaction-diffusion simulation performance through the use of just-in-time compilation. We show that these efforts have led to a growing developer base, a simpler and more robust software distribution, a wider range of supported computer architectures, a better integration of NEURON with other scientific workflows, and substantially improved performance for the simulation of biophysical and biochemical models.

Training a spiking neuronal network model of visual-motor cortex to play a virtual racket-ball game using reinforcement learning

Recent models of spiking neuronal networks have been trained to perform behaviors in static environments using a variety of learning rules, with varying degrees of biological realism. Most of these models have not been tested in dynamic visual environments where models must make predictions on future states and adjust their behavior accordingly. The models using these learning rules are often treated as black boxes, with little analysis on circuit architectures and learning mechanisms supporting optimal performance. Here we developed visual/motor spiking neuronal network models and trained them to play a virtual racket-ball game using several reinforcement learning algorithms inspired by the dopaminergic reward system. We systematically investigated how different architectures and circuit-motifs (feed-forward, recurrent, feedback) contributed to learning and performance. We also developed a new biologically-inspired learning rule that significantly enhanced performance, while reducing training time. Our models included visual areas encoding game inputs and relaying the information to motor areas, which used this information to learn to move the racket to hit the ball. Neurons in the early visual area relayed information encoding object location and motion direction across the network. Neuronal association areas encoded spatial relationships between objects in the visual scene. Motor populations received inputs from visual and association areas representing the dorsal pathway. Two populations of motor neurons generated commands to move the racket up or down. Model-generated actions updated the environment and triggered reward or punishment signals that adjusted synaptic weights so that the models could learn which actions led to reward. Here we demonstrate that our biologically-plausible learning rules were effective in training spiking neuronal network models to solve problems in dynamic environments. We used our models to dissect the circuit architectures and learning rules most effective for learning. Our model shows that learning mechanisms involving different neural circuits produce similar performance in sensory-motor tasks. In biological networks, all learning mechanisms may complement one another, accelerating the learning capabilities of animals. Furthermore, this also highlights the resilience and redundancy in biological systems.

Multiscale Computer Model of the Spinal Dorsal Horn Reveals Changes in Network Processing Associated with Chronic Pain

Pain-related sensory input is processed in the spinal dorsal horn (SDH) before being relayed to the brain. That processing profoundly influences whether stimuli are correctly or incorrectly perceived as painful. Significant advances have been made in identifying the types of excitatory and inhibitory neurons that comprise the SDH, and there is some information about how neuron types are connected, but it remains unclear how the overall circuit processes sensory input or how that processing is disrupted under chronic pain conditions. To explore SDH function, we developed a computational model of the circuit that is tightly constrained by experimental data. Our model comprises conductance-based neuron models that reproduce the characteristic firing patterns of spinal neurons. Excitatory and inhibitory neuron populations, defined by their expression of genetic markers, spiking pattern, or morphology, were synaptically connected according to available qualitative data. Using a genetic algorithm, synaptic weights were tuned to reproduce projection neuron firing rates (model output) based on primary afferent firing rates (model input) across a range of mechanical stimulus intensities. Disparate synaptic weight combinations could produce equivalent circuit function, revealing degeneracy that may underlie heterogeneous responses of different circuits to perturbations or pathologic insults. To validate our model, we verified that it responded to the reduction of inhibition (i.e., disinhibition) and ablation of specific neuron types in a manner consistent with experiments. Thus validated, our model offers a valuable resource for interpreting experimental results and testing hypotheses in silico to plan experiments for examining normal and pathologic SDH circuit function.SIGNIFICANCE STATEMENT We developed a multiscale computer model of the posterior part of spinal cord gray matter (spinal dorsal horn), which is involved in perceiving touch and pain. The model reproduces several experimental observations and makes predictions about how specific types of spinal neurons and synapses influence projection neurons that send information to the brain. Misfiring of these projection neurons can produce anomalous sensations associated with chronic pain. Our computer model will not only assist in planning future experiments, but will also be useful for developing new pharmacotherapy for chronic pain disorders, connecting the effect of drugs acting at the molecular scale with emergent properties of neurons and circuits that shape the pain experience.

NetPyNE Implementation and Scaling of the Potjans-Diesmann Cortical Microcircuit Model

The Potjans-Diesmann cortical microcircuit model is a widely used model originally implemented in NEST. Here, we reimplemented the model using NetPyNE, a high-level Python interface to the NEURON simulator, and reproduced the findings of the original publication. We also implemented a method for scaling the network size that preserves first- and second-order statistics, building on existing work on network theory. Our new implementation enabled the use of more detailed neuron models with multicompartmental morphologies and multiple biophysically realistic ion channels. This opens the model to new research, including the study of dendritic processing, the influence of individual channel parameters, the relation to local field potentials, and other multiscale interactions. The scaling method we used provides flexibility to increase or decrease the network size as needed when running these CPU-intensive detailed simulations. Finally, NetPyNE facilitates modifying or extending the model using its declarative language; optimizing model parameters; running efficient, large-scale parallelized simulations; and analyzing the model through built-in methods, including local field potential calculation and information flow measures.

Multiscale modeling meets machine learning: What can we learn?

Machine learning is increasingly recognized as a promising technology in the biological, biomedical, and behavioral sciences. There can be no argument that this technique is incredibly successful in image recognition with immediate applications in diagnostics including electrophysiology, radiology, or pathology, where we have access to massive amounts of annotated data. However, machine learning often performs poorly in prognosis, especially when dealing with sparse data. This is a field where classical physics-based simulation seems to remain irreplaceable. In this review, we identify areas in the biomedical sciences where machine learning and multiscale modeling can mutually benefit from one another: Machine learning can integrate physics-based knowledge in the form of governing equations, boundary conditions, or constraints to manage ill-posted problems and robustly handle sparse and noisy data; multiscale modeling can integrate machine learning to create surrogate models, identify system dynamics and parameters, analyze sensitivities, and quantify uncertainty to bridge the scales and understand the emergence of function. With a view towards applications in the life sciences, we discuss the state of the art of combining machine learning and multiscale modeling, identify applications and opportunities, raise open questions, and address potential challenges and limitations. We anticipate that it will stimulate discussion within the community of computational mechanics and reach out to other disciplines including mathematics, statistics, computer science, artificial intelligence, biomedicine, systems biology, and precision medicine to join forces towards creating robust and efficient models for biological systems.

Effects of Ih and TASK-like shunting current on dendritic impedance in layer 5 pyramidal-tract neurons

Pyramidal neurons in neocortex have complex input-output relationships that depend on their morphologies, ion channel distributions, and the nature of their inputs, but which cannot be replicated by simple integrate-and-fire models. The impedance properties of their dendritic arbors, such as resonance and phase shift, shape neuronal responses to synaptic inputs and provide intraneuronal functional maps reflecting their intrinsic dynamics and excitability. Experimental studies of dendritic impedance have shown that neocortical pyramidal tract neurons exhibit distance-dependent changes in resonance and impedance phase with respect to the soma. We, therefore, investigated how well several biophysically detailed multicompartment models of neocortical layer 5 pyramidal tract neurons reproduce the location-dependent impedance profiles observed experimentally. Each model tested here exhibited location-dependent impedance profiles, but most captured either the observed impedance amplitude or phase, not both. The only model that captured features from both incorporates hyperpolarization-activated cyclic nucleotide-gated (HCN) channels and a shunting current, such as that produced by Twik-related acid-sensitive K⁺ (TASK) channels. TASK-like channel density in this model was proportional to local HCN channel density. We found that although this shunting current alone is insufficient to produce resonance or realistic phase response, it modulates all features of dendritic impedance, including resonance frequencies, resonance strength, synchronous frequencies, and total inductive phase. We also explored how the interaction of HCN channel current (I_h) and a TASK-like shunting current shape synaptic potentials and produce degeneracy in dendritic impedance profiles, wherein different combinations of I_h and shunting current can produce the same impedance profile.NEW & NOTEWORTHY We simulated chirp current stimulation in the apical dendrites of 5 biophysically detailed multicompartment models of neocortical pyramidal tract neurons and found that a combination of HCN channels and TASK-like channels produced the best fit to experimental measurements of dendritic impedance. We then explored how HCN and TASK-like channels can shape the dendritic impedance as well as the voltage response to synaptic currents.

Local glutamate-mediated dendritic plateau potentials change the state of the cortical pyramidal neuron

Dendritic spikes in thin dendritic branches (basal and oblique dendrites) are traditionally inferred from spikelets measured in the cell body. Here, we used laser-spot voltage-sensitive dye imaging in cortical pyramidal neurons (rat brain slices) to investigate the voltage waveforms of dendritic potentials occurring in response to spatially restricted glutamatergic inputs. Local dendritic potentials lasted 200-500 ms and propagated to the cell body, where they caused sustained 10- to 20-mV depolarizations. Plateau potentials propagating from dendrite to soma and action potentials propagating from soma to dendrite created complex voltage waveforms in the middle of the thin basal dendrite, comprised of local sodium spikelets, local plateau potentials, and backpropagating action potentials, superimposed on each other. Our model replicated these voltage waveforms across a gradient of glutamatergic stimulation intensities. The model then predicted that somatic input resistance (Rin) and membrane time constant (tau) may be reduced during dendritic plateau potential. We then tested these model predictions in real neurons and found that the model correctly predicted the direction of Rin and tau change but not the magnitude. In summary, dendritic plateau potentials occurring in basal and oblique branches put pyramidal neurons into an activated neuronal state ("prepared state"), characterized by depolarized membrane potential and smaller but faster membrane responses. The prepared state provides a time window of 200-500 ms, during which cortical neurons are particularly excitable and capable of following afferent inputs. At the network level, this predicts that sets of cells with simultaneous plateaus would provide cellular substrate for the formation of functional neuronal ensembles.NEW & NOTEWORTHY In cortical pyramidal neurons, we recorded glutamate-mediated dendritic plateau potentials with voltage imaging and created a computer model that recreated experimental measures from dendrite and cell body. Our model made new predictions, which were then tested in experiments. Plateau potentials profoundly change neuronal state: a plateau potential triggered in one basal dendrite depolarizes the soma and shortens membrane time constant, making the cell more susceptible to firing triggered by other afferent inputs.

Simulating Large-scale Models of Brain Neuronal Circuits using Google Cloud Platform

Biophysically detailed modeling provides an unmatched method to integrate data from many disparate experimental studies, and manipulate and explore with high precision the resultin brain circuit simulation. We developed a detailed model of the brain motor cortex circuits, simulating over 10,000 biophysically detailed neurons and 30 million synaptic connections. Optimization and evaluation of the cortical model parameters and responses was achieved via parameter exploration using grid search parameter sweeps and evolutionary algorithms. This involves running tens of thousands of simulations requiring significant computational resources. This paper describes our experience in setting up and using Google Compute Platform (GCP) with Slurm to run these large-scale simulations. We describe the best practices and solutions to the issues that arose during the process, and present preliminary results from running simulations on GCP.

The SONATA data format for efficient description of large-scale network models

Increasing availability of comprehensive experimental datasets and of high-performance computing resources are driving rapid growth in scale, complexity, and biological realism of computational models in neuroscience. To support construction and simulation, as well as sharing of such large-scale models, a broadly applicable, flexible, and high-performance data format is necessary. To address this need, we have developed the Scalable Open Network Architecture TemplAte (SONATA) data format. It is designed for memory and computational efficiency and works across multiple platforms. The format represents neuronal circuits and simulation inputs and outputs via standardized files and provides much flexibility for adding new conventions or extensions. SONATA is used in multiple modeling and visualization tools, and we also provide reference Application Programming Interfaces and model examples to catalyze further adoption. SONATA format is free and open for the community to use and build upon with the goal of enabling efficient model building, sharing, and reproducibility.

Neuroscience is experiencing a rapid growth of data streams characterizing composition, connectivity, and activity of brain networks in ever increasing details. Data-driven modeling will be essential to integrate these multimodal and complex data into predictive simulations to advance our understanding of brain function and mechanisms. To enable efficient development and sharing of such large-scale models utilizing diverse data types, we have developed the Scalable Open Network Architecture TemplAte (SONATA) data format. The format represents neuronal circuits and simulation inputs and outputs via standardized files and provides much flexibility for adding new conventions or extensions. SONATA is already supported by several popular tools for model building, simulations, and visualization. It is free and open for everyone to use and build upon and will enable increased efficiency, reproducibility, and scientific exchange in the community.

Integrating machine learning and multiscale modeling-perspectives, challenges, and opportunities in the biological, biomedical, and behavioral sciences

Fueled by breakthrough technology developments, the biological, biomedical, and behavioral sciences are now collecting more data than ever before. There is a critical need for time- and cost-efficient strategies to analyze and interpret these data to advance human health. The recent rise of machine learning as a powerful technique to integrate multimodality, multifidelity data, and reveal correlations between intertwined phenomena presents a special opportunity in this regard. However, machine learning alone ignores the fundamental laws of physics and can result in ill-posed problems or non-physical solutions. Multiscale modeling is a successful strategy to integrate multiscale, multiphysics data and uncover mechanisms that explain the emergence of function. However, multiscale modeling alone often fails to efficiently combine large datasets from different sources and different levels of resolution. Here we demonstrate that machine learning and multiscale modeling can naturally complement each other to create robust predictive models that integrate the underlying physics to manage ill-posed problems and explore massive design spaces. We review the current literature, highlight applications and opportunities, address open questions, and discuss potential challenges and limitations in four overarching topical areas: ordinary differential equations, partial differential equations, data-driven approaches, and theory-driven approaches. Towards these goals, we leverage expertise in applied mathematics, computer science, computational biology, biophysics, biomechanics, engineering mechanics, experimentation, and medicine. Our multidisciplinary perspective suggests that integrating machine learning and multiscale modeling can provide new insights into disease mechanisms, help identify new targets and treatment strategies, and inform decision making for the benefit of human health.

Open Source Brain: A Collaborative Resource for Visualizing, Analyzing, Simulating, and Developing Standardized Models of Neurons and Circuits

Computational models are powerful tools for exploring the properties of complex biological systems. In neuroscience, data-driven models of neural circuits that span multiple scales are increasingly being used to understand brain function in health and disease. But their adoption and reuse has been limited by the specialist knowledge required to evaluate and use them. To address this, we have developed Open Source Brain, a platform for sharing, viewing, analyzing, and simulating standardized models from different brain regions and species. Model structure and parameters can be automatically visualized and their dynamical properties explored through browser-based simulations. Infrastructure and tools for collaborative interaction, development, and testing are also provided. We demonstrate how existing components can be reused by constructing new models of inhibition-stabilized cortical networks that match recent experimental results. These features of Open Source Brain improve the accessibility, transparency, and reproducibility of models and facilitate their reuse by the wider community.

NetPyNE, a tool for data-driven multiscale modeling of brain circuits

Biophysical modeling of neuronal networks helps to integrate and interpret rapidly growing and disparate experimental datasets at multiple scales. The NetPyNE tool (www.netpyne.org) provides both programmatic and graphical interfaces to develop data-driven multiscale network models in NEURON. NetPyNE clearly separates model parameters from implementation code. Users provide specifications at a high level via a standardized declarative language, for example connectivity rules, to create millions of cell-to-cell connections. NetPyNE then enables users to generate the NEURON network, run efficiently parallelized simulations, optimize and explore network parameters through automated batch runs, and use built-in functions for visualization and analysis - connectivity matrices, voltage traces, spike raster plots, local field potentials, and information theoretic measures. NetPyNE also facilitates model sharing by exporting and importing standardized formats (NeuroML and SONATA). NetPyNE is already being used to teach computational neuroscience students and by modelers to investigate brain regions and phenomena.

Geppetto: a reusable modular open platform for exploring neuroscience data and models

Geppetto is an open-source platform that provides generic middleware infrastructure for building both online and desktop tools for visualizing neuroscience models and data and managing simulations. Geppetto underpins a number of neuroscience applications, including Open Source Brain (OSB), Virtual Fly Brain (VFB), NEURON-UI and NetPyNE-UI. OSB is used by researchers to create and visualize computational neuroscience models described in NeuroML and simulate them through the browser. VFB is the reference hub for Drosophila melanogaster neural anatomy and imaging data including neuropil, segmented neurons, microscopy stacks and gene expression pattern data. Geppetto is also being used to build a new user interface for NEURON, a widely used neuronal simulation environment, and for NetPyNE, a Python package for network modelling using NEURON. Geppetto defines domain agnostic abstractions used by all these applications to represent their models and data and offers a set of modules and components to integrate, visualize and control simulations in a highly accessible way. The platform comprises a backend which can connect to external data sources, model repositories and simulators together with a highly customizable frontend.

This article is part of a discussion meeting issue ‘Connectome to behaviour: modelling C. elegans at cellular resolution’.

Computer modeling for pharmacological treatments for dystonia

Dystonia is a movement disorder that produces involuntary muscle contractions. Current pharmacological treatments are of limited efficacy. Dystonia, like epilepsy is a disorder involving excessive activty of motor areas including motor cortex and several causal gene mutations have been identified. In order to evaluate potential novel agents for multitarget therapy for dystonia, we have developed a computer model of cortex that includes some of the complex array of molecular interactions that, along with membrane ion channels, control cell excitability.

Evolutionary algorithm optimization of biological learning parameters in a biomimetic neuroprosthesis

Biomimetic simulation permits neuroscientists to better understand the complex neuronal dynamics of the brain. Embedding a biomimetic simulation in a closed-loop neuroprosthesis, which can read and write signals from the brain, will permit applications for amelioration of motor, psychiatric, and memory-related brain disorders. Biomimetic neuroprostheses require real-time adaptation to changes in the external environment, thus constituting an example of a dynamic data-driven application system. As model fidelity increases, so does the number of parameters and the complexity of finding appropriate parameter configurations. Instead of adapting synaptic weights via machine learning, we employed major biological learning methods: spike-timing dependent plasticity and reinforcement learning. We optimized the learning metaparameters using evolutionary algorithms, which were implemented in parallel and which used an island model approach to obtain sufficient speed. We employed these methods to train a cortical spiking model to utilize macaque brain activity, indicating a selected target, to drive a virtual musculoskeletal arm with realistic anatomical and biomechanical properties to reach to that target. The optimized system was able to reproduce macaque data from a comparable experimental motor task. These techniques can be used to efficiently tune the parameters of multiscale systems, linking realistic neuronal dynamics to behavior, and thus providing a useful tool for neuroscience and neuroprosthetics.

Optimizing computer models of corticospinal neurons to replicate in vitro dynamics

Corticospinal neurons (SPI), thick-tufted pyramidal neurons in motor cortex layer 5B that project caudally via the medullary pyramids, display distinct class-specific electrophysiological properties in vitro: strong sag with hyperpolarization, lack of adaptation, and a nearly linear frequency-current (F-I) relationship. We used our electrophysiological data to produce a pair of large archives of SPI neuron computer models in two model classes: 1) detailed models with full reconstruction; and 2) simplified models with six compartments. We used a PRAXIS and an evolutionary multiobjective optimization (EMO) in sequence to determine ion channel conductances. EMO selected good models from each of the two model classes to form the two model archives. Archived models showed tradeoffs across fitness functions. For example, parameters that produced excellent F-I fit produced a less-optimal fit for interspike voltage trajectory. Because of these tradeoffs, there was no single best model but rather models that would be best for particular usages for either single neuron or network explorations. Further exploration of exemplar models with strong F-I fit demonstrated that both the detailed and simple models produced excellent matches to the experimental data. Although dendritic ion identities and densities cannot yet be fully determined experimentally, we explored the consequences of a demonstrated proximal to distal density gradient of I_h, demonstrating that this would lead to a gradient of resonance properties with increased resonant frequencies more distally. We suggest that this dynamical feature could serve to make the cell particularly responsive to major frequency bands that differ by cortical layer.

NEW & NOTEWORTHY

We developed models of motor cortex corticospinal neurons that replicate in vitro dynamics, including hyperpolarization-induced sag and realistic firing patterns. Models demonstrated resonance in response to synaptic stimulation, with resonance frequency increasing in apical dendrites with increasing distance from soma, matching the increasing oscillation frequencies spanning deep to superficial cortical layers. This gradient may enable specific corticospinal neuron dendrites to entrain to relevant oscillations in different cortical layers, contributing to appropriate motor output commands.

Simulation Neurotechnologies for Advancing Brain Research: Parallelizing Large Networks in NEURON

Large multiscale neuronal network simulations are of increasing value as more big data are gathered about brain wiring and organization under the auspices of a current major research initiative, such as Brain Research through Advancing Innovative Neurotechnologies. The development of these models requires new simulation technologies. We describe here the current use of the NEURON simulator with message passing interface (MPI) for simulation in the domain of moderately large networks on commonly available high-performance computers (HPCs). We discuss the basic layout of such simulations, including the methods of simulation setup, the run-time spike-passing paradigm, and postsimulation data storage and data management approaches. Using the Neuroscience Gateway, a portal for computational neuroscience that provides access to large HPCs, we benchmark simulations of neuronal networks of different sizes (500-100,000 cells), and using different numbers of nodes (1-256). We compare three types of networks, composed of either Izhikevich integrate-and-fire neurons (I&F), single-compartment Hodgkin-Huxley (HH) cells, or a hybrid network with half of each. Results show simulation run time increased approximately linearly with network size and decreased almost linearly with the number of nodes. Networks with I&F neurons were faster than HH networks, although differences were small since all tested cells were point neurons with a single compartment.

25th Annual Computational Neuroscience Meeting: CNS-2016

A1 Functional advantages of cell-type heterogeneity in neural circuits

Tatyana O. Sharpee

A2 Mesoscopic modeling of propagating waves in visual cortex

Alain Destexhe

A3 Dynamics and biomarkers of mental disorders

Mitsuo Kawato

F1 Precise recruitment of spiking output at theta frequencies requires dendritic h-channels in multi-compartment models of oriens-lacunosum/moleculare hippocampal interneurons

Vladislav Sekulić, Frances K. Skinner

F2 Kernel methods in reconstruction of current sources from extracellular potentials for single cells and the whole brains

Daniel K. Wójcik, Chaitanya Chintaluri, Dorottya Cserpán, Zoltán Somogyvári

F3 The synchronized periods depend on intracellular transcriptional repression mechanisms in circadian clocks.

Jae Kyoung Kim, Zachary P. Kilpatrick, Matthew R. Bennett, Kresimir Josić

O1 Assessing irregularity and coordination of spiking-bursting rhythms in central pattern generators

Irene Elices, David Arroyo, Rafael Levi, Francisco B. Rodriguez, Pablo Varona

O2 Regulation of top-down processing by cortically-projecting parvalbumin positive neurons in basal forebrain

Eunjin Hwang, Bowon Kim, Hio-Been Han, Tae Kim, James T. McKenna, Ritchie E. Brown, Robert W. McCarley, Jee Hyun Choi

O3 Modeling auditory stream segregation, build-up and bistability

James Rankin, Pamela Osborn Popp, John Rinzel

O4 Strong competition between tonotopic neural ensembles explains pitch-related dynamics of auditory cortex evoked fields

Alejandro Tabas, André Rupp, Emili Balaguer-Ballester

O5 A simple model of retinal response to multi-electrode stimulation

Matias I. Maturana, David B. Grayden, Shaun L. Cloherty, Tatiana Kameneva, Michael R. Ibbotson, Hamish Meffin

O6 Noise correlations in V4 area correlate with behavioral performance in visual discrimination task

Veronika Koren, Timm Lochmann, Valentin Dragoi, Klaus Obermayer

O7 Input-location dependent gain modulation in cerebellar nucleus neurons

Maria Psarrou, Maria Schilstra, Neil Davey, Benjamin Torben-Nielsen, Volker Steuber

O8 Analytic solution of cable energy function for cortical axons and dendrites

Huiwen Ju, Jiao Yu, Michael L. Hines, Liang Chen, Yuguo Yu

O9 C. elegans interactome: interactive visualization of Caenorhabditis elegans worm neuronal network

Jimin Kim, Will Leahy, Eli Shlizerman

O10 Is the model any good? Objective criteria for computational neuroscience model selection

Justas Birgiolas, Richard C. Gerkin, Sharon M. Crook

O11 Cooperation and competition of gamma oscillation mechanisms

Atthaphon Viriyopase, Raoul-Martin Memmesheimer, Stan Gielen

O12 A discrete structure of the brain waves

Yuri Dabaghian, Justin DeVito, Luca Perotti

O13 Direction-specific silencing of the Drosophila gaze stabilization system

Anmo J. Kim, Lisa M. Fenk, Cheng Lyu, Gaby Maimon

O14 What does the fruit fly think about values? A model of olfactory associative learning

Chang Zhao, Yves Widmer, Simon Sprecher,Walter Senn

O15 Effects of ionic diffusion on power spectra of local field potentials (LFP)

Geir Halnes, Tuomo Mäki-Marttunen, Daniel Keller, Klas H. Pettersen,Ole A. Andreassen, Gaute T. Einevoll

O16 Large-scale cortical models towards understanding relationship between brain structure abnormalities and cognitive deficits

Yasunori Yamada

O17 Spatial coarse-graining the brain: origin of minicolumns

Moira L. Steyn-Ross, D. Alistair Steyn-Ross

O18 Modeling large-scale cortical networks with laminar structure

Jorge F. Mejias, John D. Murray, Henry Kennedy, Xiao-Jing Wang

O19 Information filtering by partial synchronous spikes in a neural population

Alexandra Kruscha, Jan Grewe, Jan Benda, Benjamin Lindner

O20 Decoding context-dependent olfactory valence in Drosophila

Laurent Badel, Kazumi Ohta, Yoshiko Tsuchimoto, Hokto Kazama

P1 Neural network as a scale-free network: the role of a hub

B. Kahng

P2 Hemodynamic responses to emotions and decisions using near-infrared spectroscopy optical imaging

Nicoladie D. Tam

P3 Phase space analysis of hemodynamic responses to intentional movement directions using functional near-infrared spectroscopy (fNIRS) optical imaging technique

Nicoladie D.Tam, Luca Pollonini, George Zouridakis

P4 Modeling jamming avoidance of weakly electric fish

Jaehyun Soh, DaeEun Kim

P5 Synergy and redundancy of retinal ganglion cells in prediction

Minsu Yoo, S. E. Palmer

P6 A neural field model with a third dimension representing cortical depth

Viviana Culmone, Ingo Bojak

P7 Network analysis of a probabilistic connectivity model of the Xenopus tadpole spinal cord

Andrea Ferrario, Robert Merrison-Hort, Roman Borisyuk

P8 The recognition dynamics in the brain

Chang Sub Kim

P9 Multivariate spike train analysis using a positive definite kernel

Taro Tezuka

P10 Synchronization of burst periods may govern slow brain dynamics during general anesthesia

Pangyu Joo

P11 The ionic basis of heterogeneity affects stochastic synchrony

Young-Ah Rho, Shawn D. Burton, G. Bard Ermentrout, Jaeseung Jeong, Nathaniel N. Urban

P12 Circular statistics of noise in spike trains with a periodic component

Petr Marsalek

P14 Representations of directions in EEG-BCI using Gaussian readouts

Hoon-Hee Kim, Seok-hyun Moon, Do-won Lee, Sung-beom Lee, Ji-yong Lee, Jaeseung Jeong

P15 Action selection and reinforcement learning in basal ganglia during reaching movements

Yaroslav I. Molkov, Khaldoun Hamade, Wondimu Teka, William H. Barnett, Taegyo Kim, Sergey Markin, Ilya A. Rybak

P17 Axon guidance: modeling axonal growth in T-Junction assay

Csaba Forro, Harald Dermutz, László Demkó, János Vörös

P19 Transient cell assembly networks encode persistent spatial memories

Yuri Dabaghian, Andrey Babichev

P20 Theory of population coupling and applications to describe high order correlations in large populations of interacting neurons

Haiping Huang

P21 Design of biologically-realistic simulations for motor control

Sergio Verduzco-Flores

P22 Towards understanding the functional impact of the behavioural variability of neurons

Filipa Dos Santos, Peter Andras

P23 Different oscillatory dynamics underlying gamma entrainment deficits in schizophrenia

Christoph Metzner, Achim Schweikard, Bartosz Zurowski

P24 Memory recall and spike frequency adaptation

James P. Roach, Leonard M. Sander, Michal R. Zochowski

P25 Stability of neural networks and memory consolidation preferentially occur near criticality

Quinton M. Skilling, Nicolette Ognjanovski, Sara J. Aton, Michal Zochowski

P26 Stochastic Oscillation in Self-Organized Critical States of Small Systems: Sensitive Resting State in Neural Systems

Sheng-Jun Wang, Guang Ouyang, Jing Guang, Mingsha Zhang, K. Y. Michael Wong, Changsong Zhou

P27 Neurofield: a C++ library for fast simulation of 2D neural field models

Peter A. Robinson, Paula Sanz-Leon, Peter M. Drysdale, Felix Fung, Romesh G. Abeysuriya, Chris J. Rennie, Xuelong Zhao

P28 Action-based grounding: Beyond encoding/decoding in neural code

Yoonsuck Choe, Huei-Fang Yang

P29 Neural computation in a dynamical system with multiple time scales

Yuanyuan Mi, Xiaohan Lin, Si Wu

P30 Maximum entropy models for 3D layouts of orientation selectivity

Joscha Liedtke, Manuel Schottdorf, Fred Wolf

P31 A behavioral assay for probing computations underlying curiosity in rodents

Yoriko Yamamura, Jeffery R. Wickens

P32 Using statistical sampling to balance error function contributions to optimization of conductance-based models

Timothy Rumbell, Julia Ramsey, Amy Reyes, Danel Draguljić, Patrick R. Hof, Jennifer Luebke, Christina M. Weaver

P33 Exploration and implementation of a self-growing and self-organizing neuron network building algorithm

Hu He, Xu Yang, Hailin Ma, Zhiheng Xu, Yuzhe Wang

P34 Disrupted resting state brain network in obese subjects: a data-driven graph theory analysis

Kwangyeol Baek, Laurel S. Morris, Prantik Kundu, Valerie Voon

P35 Dynamics of cooperative excitatory and inhibitory plasticity

Everton J. Agnes, Tim P. Vogels

P36 Frequency-dependent oscillatory signal gating in feed-forward networks of integrate-and-fire neurons

William F. Podlaski, Tim P. Vogels

P37 Phenomenological neural model for adaptation of neurons in area IT

Martin Giese, Pradeep Kuravi, Rufin Vogels

P38 ICGenealogy: towards a common topology of neuronal ion channel function and genealogy in model and experiment

Alexander Seeholzer, William Podlaski, Rajnish Ranjan, Tim Vogels

P39 Temporal input discrimination from the interaction between dynamic synapses and neural subthreshold oscillations

Joaquin J. Torres, Fabiano Baroni, Roberto Latorre, Pablo Varona

P40 Different roles for transient and sustained activity during active visual processing

Bart Gips, Eric Lowet, Mark J. Roberts, Peter de Weerd, Ole Jensen, Jan van der Eerden

P41 Scale-free functional networks of 2D Ising model are highly robust against structural defects: neuroscience implications

Abdorreza Goodarzinick, Mohammad D. Niry, Alireza Valizadeh

P42 High frequency neuron can facilitate propagation of signal in neural networks

Aref Pariz, Shervin S. Parsi, Alireza Valizadeh

P43 Investigating the effect of Alzheimer’s disease related amyloidopathy on gamma oscillations in the CA1 region of the hippocampus

Julia M. Warburton, Lucia Marucci, Francesco Tamagnini, Jon Brown, Krasimira Tsaneva-Atanasova

P44 Long-tailed distributions of inhibitory and excitatory weights in a balanced network with eSTDP and iSTDP

Florence I. Kleberg, Jochen Triesch

P45 Simulation of EMG recording from hand muscle due to TMS of motor cortex

Bahar Moezzi, Nicolangelo Iannella, Natalie Schaworonkow, Lukas Plogmacher, Mitchell R. Goldsworthy, Brenton Hordacre, Mark D. McDonnell, Michael C. Ridding, Jochen Triesch

P46 Structure and dynamics of axon network formed in primary cell culture

Martin Zapotocky, Daniel Smit, Coralie Fouquet, Alain Trembleau

P47 Efficient signal processing and sampling in random networks that generate variability

Sakyasingha Dasgupta, Isao Nishikawa, Kazuyuki Aihara, Taro Toyoizumi

P48 Modeling the effect of riluzole on bursting in respiratory neural networks

Daniel T. Robb, Nick Mellen, Natalia Toporikova

P49 Mapping relaxation training using effective connectivity analysis

Rongxiang Tang, Yi-Yuan Tang

P50 Modeling neuron oscillation of implicit sequence learning

Guangsheng Liang, Seth A. Kiser, James H. Howard, Jr., Yi-Yuan Tang

P51 The role of cerebellar short-term synaptic plasticity in the pathology and medication of downbeat nystagmus

Julia Goncharenko, Neil Davey, Maria Schilstra, Volker Steuber

P52 Nonlinear response of noisy neurons

Sergej O. Voronenko, Benjamin Lindner

P53 Behavioral embedding suggests multiple chaotic dimensions underlie C. elegans locomotion

Tosif Ahamed, Greg Stephens

P54 Fast and scalable spike sorting for large and dense multi-electrodes recordings

Pierre Yger, Baptiste Lefebvre, Giulia Lia Beatrice Spampinato, Elric Esposito, Marcel Stimberg et Olivier Marre

P55 Sufficient sampling rates for fast hand motion tracking

Hansol Choi, Min-Ho Song

P56 Linear readout of object manifolds

SueYeon Chung, Dan D. Lee, Haim Sompolinsky

P57 Differentiating models of intrinsic bursting and rhythm generation of the respiratory pre-Bötzinger complex using phase response curves

Ryan S. Phillips, Jeffrey Smith

P58 The effect of inhibitory cell network interactions during theta rhythms on extracellular field potentials in CA1 hippocampus

Alexandra Pierri Chatzikalymniou, Katie Ferguson, Frances K. Skinner

P59 Expansion recoding through sparse sampling in the cerebellar input layer speeds learning

N. Alex Cayco Gajic, Claudia Clopath, R. Angus Silver

P60 A set of curated cortical models at multiple scales on Open Source Brain

Padraig Gleeson, Boris Marin, Sadra Sadeh, Adrian Quintana, Matteo Cantarelli, Salvador Dura-Bernal, William W. Lytton, Andrew Davison, R. Angus Silver

P61 A synaptic story of dynamical information encoding in neural adaptation

Luozheng Li, Wenhao Zhang, Yuanyuan Mi, Dahui Wang, Si Wu

P62 Physical modeling of rule-observant rodent behavior

Youngjo Song, Sol Park, Ilhwan Choi, Jaeseung Jeong, Hee-sup Shin

P64 Predictive coding in area V4 and prefrontal cortex explains dynamic discrimination of partially occluded shapes

Hannah Choi, Anitha Pasupathy, Eric Shea-Brown

P65 Stability of FORCE learning on spiking and rate-based networks

Dongsung Huh, Terrence J. Sejnowski

P66 Stabilising STDP in striatal neurons for reliable fast state recognition in noisy environments

Simon M. Vogt, Arvind Kumar, Robert Schmidt

P67 Electrodiffusion in one- and two-compartment neuron models for characterizing cellular effects of electrical stimulation

Stephen Van Wert, Steven J. Schiff

P68 STDP improves speech recognition capabilities in spiking recurrent circuits parameterized via differential evolution Markov Chain Monte Carlo

Richard Veale, Matthias Scheutz

P69 Bidirectional transformation between dominant cortical neural activities and phase difference distributions

Sang Wan Lee

P70 Maturation of sensory networks through homeostatic structural plasticity

Júlia Gallinaro, Stefan Rotter

P71 Corticothalamic dynamics: structure, number of solutions and stability of steady-state solutions in the space of synaptic couplings

Paula Sanz-Leon, Peter A. Robinson

P72 Optogenetic versus electrical stimulation of the parkinsonian basal ganglia. Computational study

Leonid L. Rubchinsky, Chung Ching Cheung, Shivakeshavan Ratnadurai-Giridharan

P73 Exact spike-timing distribution reveals higher-order interactions of neurons

Safura Rashid Shomali, Majid Nili Ahmadabadi, Hideaki Shimazaki, S. Nader Rasuli

P74 Neural mechanism of visual perceptual learning using a multi-layered neural network

Xiaochen Zhao, Malte J. Rasch

P75 Inferring collective spiking dynamics from mostly unobserved systems

Jens Wilting, Viola Priesemann

P76 How to infer distributions in the brain from subsampled observations

Anna Levina, Viola Priesemann

P77 Influences of embedding and estimation strategies on the inferred memory of single spiking neurons

Lucas Rudelt, Joseph T. Lizier, Viola Priesemann

P78 A nearest-neighbours based estimator for transfer entropy between spike trains

Joseph T. Lizier, Richard E. Spinney, Mikail Rubinov, Michael Wibral, Viola Priesemann

P79 Active learning of psychometric functions with multinomial logistic models

Ji Hyun Bak, Jonathan Pillow

P81 Inferring low-dimensional network dynamics with variational latent Gaussian process

Yuan Zaho, Il Memming Park

P82 Computational investigation of energy landscapes in the resting state subcortical brain network

Jiyoung Kang, Hae-Jeong Park

P83 Local repulsive interaction between retinal ganglion cells can generate a consistent spatial periodicity of orientation map

Jaeson Jang, Se-Bum Paik

P84 Phase duration of bistable perception reveals intrinsic time scale of perceptual decision under noisy condition

Woochul Choi, Se-Bum Paik

P85 Feedforward convergence between retina and primary visual cortex can determine the structure of orientation map

Changju Lee, Jaeson Jang, Se-Bum Paik

P86 Computational method classifying neural network activity patterns for imaging data

Min Song, Hyeonsu Lee, Se-Bum Paik

P87 Symmetry of spike-timing-dependent-plasticity kernels regulates volatility of memory

Youngjin Park, Woochul Choi^, Se-Bum Paik

P88 Effects of time-periodic coupling strength on the first-spike latency dynamics of a scale-free network of stochastic Hodgkin-Huxley neurons

Ergin Yilmaz, Veli Baysal, Mahmut Ozer

P89 Spectral properties of spiking responses in V1 and V4 change within the trial and are highly relevant for behavioral performance

Veronika Koren, Klaus Obermayer

P90 Methods for building accurate models of individual neurons

Daniel Saska, Thomas Nowotny

P91 A full size mathematical model of the early olfactory system of honeybees

Ho Ka Chan, Alan Diamond, Thomas Nowotny

P92 Stimulation-induced tuning of ongoing oscillations in spiking neural networks

Christoph S. Herrmann, Micah M. Murray, Silvio Ionta, Axel Hutt, Jérémie Lefebvre

P93 Decision-specific sequences of neural activity in balanced random networks driven by structured sensory input

Philipp Weidel, Renato Duarte, Abigail Morrison

P94 Modulation of tuning induced by abrupt reduction of SST cell activity

Jung H. Lee, Ramakrishnan Iyer, Stefan Mihalas

P95 The functional role of VIP cell activation during locomotion

Jung H. Lee, Ramakrishnan Iyer, Christof Koch, Stefan Mihalas

P96 Stochastic inference with spiking neural networks

Mihai A. Petrovici, Luziwei Leng, Oliver Breitwieser, David Stöckel, Ilja Bytschok, Roman Martel, Johannes Bill, Johannes Schemmel, Karlheinz Meier

P97 Modeling orientation-selective electrical stimulation with retinal prostheses

Timothy B. Esler, Anthony N. Burkitt, David B. Grayden, Robert R. Kerr, Bahman Tahayori, Hamish Meffin

P98 Ion channel noise can explain firing correlation in auditory nerves

Bahar Moezzi, Nicolangelo Iannella, Mark D. McDonnell

P99 Limits of temporal encoding of thalamocortical inputs in a neocortical microcircuit

Max Nolte, Michael W. Reimann, Eilif Muller, Henry Markram

P100 On the representation of arm reaching movements: a computational model

Antonio Parziale, Rosa Senatore, Angelo Marcelli

P101 A computational model for investigating the role of cerebellum in acquisition and retention of motor behavior

Rosa Senatore, Antonio Parziale, Angelo Marcelli

P102 The emergence of semantic categories from a large-scale brain network of semantic knowledge

K. Skiker, M. Maouene

P103 Multiscale modeling of M1 multitarget pharmacotherapy for dystonia

Samuel A. Neymotin, Salvador Dura-Bernal, Alexandra Seidenstein, Peter Lakatos, Terence D. Sanger, William W. Lytton

P104 Effect of network size on computational capacity

Salvador Dura-Bernal, Rosemary J. Menzies, Campbell McLauchlan, Sacha J. van Albada, David J. Kedziora, Samuel Neymotin, William W. Lytton, Cliff C. Kerr

P105 NetPyNE: a Python package for NEURON to facilitate development and parallel simulation of biological neuronal networks

Salvador Dura-Bernal, Benjamin A. Suter, Samuel A. Neymotin, Cliff C. Kerr, Adrian Quintana, Padraig Gleeson, Gordon M. G. Shepherd, William W. Lytton

P107 Inter-areal and inter-regional inhomogeneity in co-axial anisotropy of Cortical Point Spread in human visual areas

Juhyoung Ryu, Sang-Hun Lee

P108 Two bayesian quanta of uncertainty explain the temporal dynamics of cortical activity in the non-sensory areas during bistable perception

Joonwon Lee, Sang-Hun Lee

P109 Optimal and suboptimal integration of sensory and value information in perceptual decision making

Hyang Jung Lee, Sang-Hun Lee

P110 A Bayesian algorithm for phoneme Perception and its neural implementation

Daeseob Lim, Sang-Hun Lee

P111 Complexity of EEG signals is reduced during unconsciousness induced by ketamine and propofol

Jisung Wang, Heonsoo Lee

P112 Self-organized criticality of neural avalanche in a neural model on complex networks

Nam Jung, Le Anh Quang, Seung Eun Maeng, Tae Ho Lee, Jae Woo Lee

P113 Dynamic alterations in connection topology of the hippocampal network during ictal-like epileptiform activity in an in vitro rat model

Chang-hyun Park, Sora Ahn, Jangsup Moon, Yun Seo Choi, Juhee Kim, Sang Beom Jun, Seungjun Lee, Hyang Woon Lee

P114 Computational model to replicate seizure suppression effect by electrical stimulation

Sora Ahn, Sumin Jo, Eunji Jun, Suin Yu, Hyang Woon Lee, Sang Beom Jun, Seungjun Lee

P115 Identifying excitatory and inhibitory synapses in neuronal networks from spike trains using sorted local transfer entropy

Felix Goetze, Pik-Yin Lai

P116 Neural network model for obstacle avoidance based on neuromorphic computational model of boundary vector cell and head direction cell

Seonghyun Kim, Jeehyun Kwag

P117 Dynamic gating of spike pattern propagation by Hebbian and anti-Hebbian spike timing-dependent plasticity in excitatory feedforward network model

Hyun Jae Jang, Jeehyun Kwag

P118 Inferring characteristics of input correlations of cells exhibiting up-down state transitions in the rat striatum

Marko Filipović, Ramon Reig, Ad Aertsen, Gilad Silberberg, Arvind Kumar

P119 Graph properties of the functional connected brain under the influence of Alzheimer’s disease

Claudia Bachmann, Simone Buttler, Heidi Jacobs, Kim Dillen, Gereon R. Fink, Juraj Kukolja, Abigail Morrison

P120 Learning sparse representations in the olfactory bulb

Daniel Kepple, Hamza Giaffar, Dima Rinberg, Steven Shea, Alex Koulakov

P121 Functional classification of homologous basal-ganglia networks

Jyotika Bahuguna,Tom Tetzlaff, Abigail Morrison, Arvind Kumar, Jeanette Hellgren Kotaleski

P122 Short term memory based on multistability

Tim Kunze, Andre Peterson, Thomas Knösche

P123 A physiologically plausible, computationally efficient model and simulation software for mammalian motor units

Minjung Kim, Hojeong Kim

P125 Decoding laser-induced somatosensory information from EEG

Ji Sung Park, Ji Won Yeon, Sung-Phil Kim

P126 Phase synchronization of alpha activity for EEG-based personal authentication

Jae-Hwan Kang, Chungho Lee, Sung-Phil Kim

P129 Investigating phase-lags in sEEG data using spatially distributed time delays in a large-scale brain network model

Andreas Spiegler, Spase Petkoski, Matias J. Palva, Viktor K. Jirsa

P130 Epileptic seizures in the unfolding of a codimension-3 singularity

Maria L. Saggio, Silvan F. Siep, Andreas Spiegler, William C. Stacey, Christophe Bernard, Viktor K. Jirsa

P131 Incremental dimensional exploratory reasoning under multi-dimensional environment

Oh-hyeon Choung, Yong Jeong

P132 A low-cost model of eye movements and memory in personal visual cognition

Yong-il Lee, Jaeseung Jeong

P133 Complex network analysis of structural connectome of autism spectrum disorder patients

Su Hyun Kim, Mir Jeong, Jaeseung Jeong

P134 Cognitive motives and the neural correlates underlying human social information transmission, gossip

Jeungmin Lee, Jaehyung Kwon, Jerald D. Kralik, Jaeseung Jeong

P135 EEG hyperscanning detects neural oscillation for the social interaction during the economic decision-making

Jaehwan Jahng, Dong-Uk Hwang, Jaeseung Jeong

P136 Detecting purchase decision based on hyperfrontality of the EEG

Jae-Hyung Kwon, Sang-Min Park, Jaeseung Jeong

P137 Vulnerability-based critical neurons, synapses, and pathways in the Caenorhabditis elegans connectome

Seongkyun Kim, Hyoungkyu Kim, Jerald D. Kralik, Jaeseung Jeong

P138 Motif analysis reveals functionally asymmetrical neurons in C. elegans

Pyeong Soo Kim, Seongkyun Kim, Hyoungkyu Kim, Jaeseung Jeong

P139 Computational approach to preference-based serial decision dynamics: do temporal discounting and working memory affect it?

Sangsup Yoon, Jaehyung Kwon, Sewoong Lim, Jaeseung Jeong

P141 Social stress induced neural network reconfiguration affects decision making and learning in zebrafish

Choongseok Park, Thomas Miller, Katie Clements, Sungwoo Ahn, Eoon Hye Ji, Fadi A. Issa

P142 Descriptive, generative, and hybrid approaches for neural connectivity inference from neural activity data

JeongHun Baek, Shigeyuki Oba, Junichiro Yoshimoto, Kenji Doya, Shin Ishii

P145 Divergent-convergent synaptic connectivities accelerate coding in multilayered sensory systems

Thiago S. Mosqueiro, Martin F. Strube-Bloss, Brian Smith, Ramon Huerta

P146 Swinging networks

Michal Hadrava, Jaroslav Hlinka

P147 Inferring dynamically relevant motifs from oscillatory stimuli: challenges, pitfalls, and solutions

Hannah Bos, Moritz Helias

P148 Spatiotemporal mapping of brain network dynamics during cognitive tasks using magnetoencephalography and deep learning

Charles M. Welzig, Zachary J. Harper

P149 Multiscale complexity analysis for the segmentation of MRI images

Won Sup Kim, In-Seob Shin, Hyeon-Man Baek, Seung Kee Han

P150 A neuro-computational model of emotional attention

René Richter, Julien Vitay, Frederick Beuth, Fred H. Hamker

P151 Multi-site delayed feedback stimulation in parkinsonian networks

Kelly Toppin, Yixin Guo

P152 Bistability in Hodgkin–Huxley-type equations

Tatiana Kameneva, Hamish Meffin, Anthony N. Burkitt, David B. Grayden

P153 Phase changes in postsynaptic spiking due to synaptic connectivity and short term plasticity: mathematical analysis of frequency dependency

Mark D. McDonnell, Bruce P. Graham

P154 Quantifying resilience patterns in brain networks: the importance of directionality

Penelope J. Kale, Leonardo L. Gollo

P155 Dynamics of rate-model networks with separate excitatory and inhibitory populations

Merav Stern, L. F. Abbott

P156 A model for multi-stable dynamics in action recognition modulated by integration of silhouette and shading cues

Leonid A. Fedorov, Martin A. Giese

P157 Spiking model for the interaction between action recognition and action execution

Mohammad Hovaidi Ardestani, Martin Giese

P158 Surprise-modulated belief update: how to learn within changing environments?

Mohammad Javad Faraji, Kerstin Preuschoff, Wulfram Gerstner

P159 A fast, stochastic and adaptive model of auditory nerve responses to cochlear implant stimulation

Margriet J. van Gendt, Jeroen J. Briaire, Randy K. Kalkman, Johan H. M. Frijns

P160 Quantitative comparison of graph theoretical measures of simulated and empirical functional brain networks

Won Hee Lee, Sophia Frangou

P161 Determining discriminative properties of fMRI signals in schizophrenia using highly comparative time-series analysis

Ben D. Fulcher, Patricia H. P. Tran, Alex Fornito

P162 Emergence of narrowband LFP oscillations from completely asynchronous activity during seizures and high-frequency oscillations

Stephen V. Gliske, William C. Stacey, Eugene Lim, Katherine A. Holman, Christian G. Fink

P163 Neuronal diversity in structure and function: cross-validation of anatomical and physiological classification of retinal ganglion cells in the mouse

Jinseop S. Kim, Shang Mu, Kevin L. Briggman, H. Sebastian Seung, the EyeWirers

P164 Analysis and modelling of transient firing rate changes in area MT in response to rapid stimulus feature changes

Detlef Wegener, Lisa Bohnenkamp, Udo A. Ernst

P165 Step-wise model fitting accounting for high-resolution spatial measurements: construction of a layer V pyramidal cell model with reduced morphology

Tuomo Mäki-Marttunen, Geir Halnes, Anna Devor, Christoph Metzner, Anders M. Dale, Ole A. Andreassen, Gaute T. Einevoll

P166 Contributions of schizophrenia-associated genes to neuron firing and cardiac pacemaking: a polygenic modeling approach

Tuomo Mäki-Marttunen, Glenn T. Lines, Andy Edwards, Aslak Tveito, Anders M. Dale, Gaute T. Einevoll, Ole A. Andreassen

P167 Local field potentials in a 4 × 4 mm² multi-layered network model

Espen Hagen, Johanna Senk, Sacha J. van Albada, Markus Diesmann

P168 A spiking network model explains multi-scale properties of cortical dynamics

Maximilian Schmidt, Rembrandt Bakker, Kelly Shen, Gleb Bezgin, Claus-Christian Hilgetag, Markus Diesmann, Sacha Jennifer van Albada

P169 Using joint weight-delay spike-timing dependent plasticity to find polychronous neuronal groups

Haoqi Sun, Olga Sourina, Guang-Bin Huang, Felix Klanner, Cornelia Denk

P170 Tensor decomposition reveals RSNs in simulated resting state fMRI

Katharina Glomb, Adrián Ponce-Alvarez, Matthieu Gilson, Petra Ritter, Gustavo Deco

P171 Getting in the groove: testing a new model-based method for comparing task-evoked vs resting-state activity in fMRI data on music listening

Matthieu Gilson, Maria AG Witek, Eric F. Clarke, Mads Hansen, Mikkel Wallentin, Gustavo Deco, Morten L. Kringelbach^, Peter Vuust

P172 STochastic engine for pathway simulation (STEPS) on massively parallel processors

Guido Klingbeil, Erik De Schutter

P173 Toolkit support for complex parallel spatial stochastic reaction–diffusion simulation in STEPS

Weiliang Chen, Erik De Schutter

P174 Modeling the generation and propagation of Purkinje cell dendritic spikes caused by parallel fiber synaptic input

Yunliang Zang, Erik De Schutter

P175 Dendritic morphology determines how dendrites are organized into functional subunits

Sungho Hong, Akira Takashima, Erik De Schutter

P176 A model of Ca²⁺/calmodulin-dependent protein kinase II activity in long term depression at Purkinje cells

Criseida Zamora, Andrew R. Gallimore, Erik De Schutter

P177 Reward-modulated learning of population-encoded vectors for insect-like navigation in embodied agents

Dennis Goldschmidt, Poramate Manoonpong, Sakyasingha Dasgupta

P178 Data-driven neural models part II: connectivity patterns of human seizures

Philippa J. Karoly, Dean R. Freestone, Daniel Soundry, Levin Kuhlmann, Liam Paninski, Mark Cook

P179 Data-driven neural models part I: state and parameter estimation

Dean R. Freestone, Philippa J. Karoly, Daniel Soundry, Levin Kuhlmann, Mark Cook

P180 Spectral and spatial information processing in human auditory streaming

Jaejin Lee, Yonatan I. Fishman, Yale E. Cohen

P181 A tuning curve for the global effects of local perturbations in neural activity: Mapping the systems-level susceptibility of the brain

Leonardo L. Gollo, James A. Roberts, Luca Cocchi

P182 Diverse homeostatic responses to visual deprivation mediated by neural ensembles

Yann Sweeney, Claudia Clopath

P183 Opto-EEG: a novel method for investigating functional connectome in mouse brain based on optogenetics and high density electroencephalography

Soohyun Lee, Woo-Sung Jung, Jee Hyun Choi

P184 Biphasic responses of frontal gamma network to repetitive sleep deprivation during REM sleep

Bowon Kim, Youngsoo Kim, Eunjin Hwang, Jee Hyun Choi

P185 Brain-state correlate and cortical connectivity for frontal gamma oscillations in top-down fashion assessed by auditory steady-state response

Younginha Jung, Eunjin Hwang, Yoon-Kyu Song, Jee Hyun Choi

P186 Neural field model of localized orientation selective activation in V1

James Rankin, Frédéric Chavane

P187 An oscillatory network model of Head direction and Grid cells using locomotor inputs

Karthik Soman, Vignesh Muralidharan, V. Srinivasa Chakravarthy

P188 A computational model of hippocampus inspired by the functional architecture of basal ganglia

Karthik Soman, Vignesh Muralidharan, V. Srinivasa Chakravarthy

P189 A computational architecture to model the microanatomy of the striatum and its functional properties

Sabyasachi Shivkumar, Vignesh Muralidharan, V. Srinivasa Chakravarthy

P190 A scalable cortico-basal ganglia model to understand the neural dynamics of targeted reaching

Vignesh Muralidharan, Alekhya Mandali, B. Pragathi Priyadharsini, Hima Mehta, V. Srinivasa Chakravarthy

P191 Emergence of radial orientation selectivity from synaptic plasticity

Catherine E. Davey, David B. Grayden, Anthony N. Burkitt

P192 How do hidden units shape effective connections between neurons?

Braden A. W. Brinkman, Tyler Kekona, Fred Rieke, Eric Shea-Brown, Michael Buice

P193 Characterization of neural firing in the presence of astrocyte-synapse signaling

Maurizio De Pittà, Hugues Berry, Nicolas Brunel

P194 Metastability of spatiotemporal patterns in a large-scale network model of brain dynamics

James A. Roberts, Leonardo L. Gollo, Michael Breakspear

P195 Comparison of three methods to quantify detection and discrimination capacity estimated from neural population recordings

Gary Marsat, Jordan Drew, Phillip D. Chapman, Kevin C. Daly, Samual P. Bradley

P196 Quantifying the constraints for independent evoked and spontaneous NMDA receptor mediated synaptic transmission at individual synapses

Sat Byul Seo, Jianzhong Su, Ege T. Kavalali, Justin Blackwell

P199 Gamma oscillation via adaptive exponential integrate-and-fire neurons

LieJune Shiau, Laure Buhry, Kanishka Basnayake

P200 Visual face representations during memory retrieval compared to perception

Sue-Hyun Lee, Brandon A. Levy, Chris I. Baker

P201 Top-down modulation of sequential activity within packets modeled using avalanche dynamics

Timothée Leleu, Kazuyuki Aihara

Q28 An auto-encoder network realizes sparse features under the influence of desynchronized vascular dynamics

Ryan T. Philips, Karishma Chhabria, V. Srinivasa Chakravarthy

Multitarget Multiscale Simulation for Pharmacological Treatment of Dystonia in Motor Cortex

A large number of physiomic pathologies can produce hyperexcitability in cortex. Depending on severity, cortical hyperexcitability may manifest clinically as a hyperkinetic movement disorder or as epilpesy. We focus here on dystonia, a movement disorder that produces involuntary muscle contractions and involves pathology in multiple brain areas including basal ganglia, thalamus, cerebellum, and sensory and motor cortices. Most research in dystonia has focused on basal ganglia, while much pharmacological treatment is provided directly at muscles to prevent contraction. Motor cortex is another potential target for therapy that exhibits pathological dynamics in dystonia, including heightened activity and altered beta oscillations. We developed a multiscale model of primary motor cortex, ranging from molecular, up to cellular, and network levels, containing 1715 compartmental model neurons with multiple ion channels and intracellular molecular dynamics. We wired the model based on electrophysiological data obtained from mouse motor cortex circuit mapping experiments. We used the model to reproduce patterns of heightened activity seen in dystonia by applying independent random variations in parameters to identify pathological parameter sets. These models demonstrated degeneracy, meaning that there were many ways of obtaining the pathological syndrome. There was no single parameter alteration which would consistently distinguish pathological from physiological dynamics. At higher dimensions in parameter space, we were able to use support vector machines to distinguish the two patterns in different regions of space and thereby trace multitarget routes from dystonic to physiological dynamics. These results suggest the use of in silico models for discovery of multitarget drug cocktails.

Restoring Behavior via Inverse Neurocontroller in a Lesioned Cortical Spiking Model Driving a Virtual Arm

Neural stimulation can be used as a tool to elicit natural sensations or behaviors by modulating neural activity. This can be potentially used to mitigate the damage of brain lesions or neural disorders. However, in order to obtain the optimal stimulation sequences, it is necessary to develop neural control methods, for example by constructing an inverse model of the target system. For real brains, this can be very challenging, and often unfeasible, as it requires repeatedly stimulating the neural system to obtain enough probing data, and depends on an unwarranted assumption of stationarity. By contrast, detailed brain simulations may provide an alternative testbed for understanding the interactions between ongoing neural activity and external stimulation. Unlike real brains, the artificial system can be probed extensively and precisely, and detailed output information is readily available. Here we employed a spiking network model of sensorimotor cortex trained to drive a realistic virtual musculoskeletal arm to reach a target. The network was then perturbed, in order to simulate a lesion, by either silencing neurons or removing synaptic connections. All lesions led to significant behvaioral impairments during the reaching task. The remaining cells were then systematically probed with a set of single and multiple-cell stimulations, and results were used to build an inverse model of the neural system. The inverse model was constructed using a kernel adaptive filtering method, and was used to predict the neural stimulation pattern required to recover the pre-lesion neural activity. Applying the derived neurostimulation to the lesioned network improved the reaching behavior performance. This work proposes a novel neurocontrol method, and provides theoretical groundwork on the use biomimetic brain models to develop and evaluate neurocontrollers that restore the function of damaged brain regions and the corresponding motor behaviors.

Towards real-time communication between in vivo neurophysiological data sources and simulator-based brain biomimetic models

Development of more sophisticated implantable brain-machine interface (BMI) will require both interpretation of the neurophysiological data being measured and subsequent determination of signals to be delivered back to the brain. Computational models are the heart of the machine of BMI and therefore an essential tool in both of these processes. One approach is to utilize brain biomimetic models (BMMs) to develop and instantiate these algorithms. These then must be connected as hybrid systems in order to interface the BMM with in vivo data acquisition devices and prosthetic devices. The combined system then provides a test bed for neuroprosthetic rehabilitative solutions and medical devices for the repair and enhancement of damaged brain. We propose here a computer network-based design for this purpose, detailing its internal modules and data flows. We describe a prototype implementation of the design, enabling interaction between the Plexon Multichannel Acquisition Processor (MAP) server, a commercial tool to collect signals from microelectrodes implanted in a live subject and a BMM, a NEURON-based model of sensorimotor cortex capable of controlling a virtual arm. The prototype implementation supports an online mode for real-time simulations, as well as an offline mode for data analysis and simulations without real-time constraints, and provides binning operations to discretize continuous input to the BMM and filtering operations for dealing with noise. Evaluation demonstrated that the implementation successfully delivered monkey spiking activity to the BMM through LAN environments, respecting real-time constraints.

Cortical Spiking Network Interfaced with Virtual Musculoskeletal Arm and Robotic Arm

Embedding computational models in the physical world is a critical step towards constraining their behavior and building practical applications. Here we aim to drive a realistic musculoskeletal arm model using a biomimetic cortical spiking model, and make a robot arm reproduce the same trajectories in real time. Our cortical model consisted of a 3-layered cortex, composed of several hundred spiking model-neurons, which display physiologically realistic dynamics. We interconnected the cortical model to a two-joint musculoskeletal model of a human arm, with realistic anatomical and biomechanical properties. The virtual arm received muscle excitations from the neuronal model, and fed back proprioceptive information, forming a closed-loop system. The cortical model was trained using spike timing-dependent reinforcement learning to drive the virtual arm in a 2D reaching task. Limb position was used to simultaneously control a robot arm using an improved network interface. Virtual arm muscle activations responded to motoneuron firing rates, with virtual arm muscles lengths encoded via population coding in the proprioceptive population. After training, the virtual arm performed reaching movements which were smoother and more realistic than those obtained using a simplistic arm model. This system provided access to both spiking network properties and to arm biophysical properties, including muscle forces. The use of a musculoskeletal virtual arm and the improved control system allowed the robot arm to perform movements which were smoother than those reported in our previous paper using a simplistic arm. This work provides a novel approach consisting of bidirectionally connecting a cortical model to a realistic virtual arm, and using the system output to drive a robotic arm in real time. Our techniques are applicable to the future development of brain neuroprosthetic control systems, and may enable enhanced brain-machine interfaces with the possibility for finer control of limb prosthetics.

Towards a real-time interface between a biomimetic model of sensorimotor cortex and a robotic arm

Brain-machine interfaces can greatly improve the performance of prosthetics. Utilizing biomimetic neuronal modeling in brain machine interfaces (BMI) offers the possibility of providing naturalistic motor-control algorithms for control of a robotic limb. This will allow finer control of a robot, while also giving us new tools to better understand the brain's use of electrical signals. However, the biomimetic approach presents challenges in integrating technologies across multiple hardware and software platforms, so that the different components can communicate in real-time. We present the first steps in an ongoing effort to integrate a biomimetic spiking neuronal model of motor learning with a robotic arm. The biomimetic model (BMM) was used to drive a simple kinematic two-joint virtual arm in a motor task requiring trial-and-error convergence on a single target. We utilized the output of this model in real time to drive mirroring motion of a Barrett Technology WAM robotic arm through a user datagram protocol (UDP) interface. The robotic arm sent back information on its joint positions, which was then used by a visualization tool on the remote computer to display a realistic 3D virtual model of the moving robotic arm in real time. This work paves the way towards a full closed-loop biomimetic brain-effector system that can be incorporated in a neural decoder for prosthetic control, to be used as a platform for developing biomimetic learning algorithms for controlling real-time devices.

Multimodal integration of micro-Doppler sonar and auditory signals for behavior classification with convolutional networks

The ability to recognize the behavior of individuals is of great interest in the general field of safety (e.g. building security, crowd control, transport analysis, independent living for the elderly). Here we report a new real-time acoustic system for human action and behavior recognition that integrates passive audio and active micro-Doppler sonar signatures over multiple time scales. The system architecture is based on a six-layer convolutional neural network, trained and evaluated using a dataset of 10 subjects performing seven different behaviors. Probabilistic combination of system output through time for each modality separately yields 94% (passive audio) and 91% (micro-Doppler sonar) correct behavior classification; probabilistic multimodal integration increases classification performance to 98%. This study supports the efficacy of micro-Doppler sonar systems in characterizing human actions, which can then be efficiently classified using ConvNets. It also demonstrates that the integration of multiple sources of acoustic information can significantly improve the system's performance.

Top-down feedback in an HMAX-like cortical model of object perception based on hierarchical Bayesian networks and belief propagation

Hierarchical generative models, such as Bayesian networks, and belief propagation have been shown to provide a theoretical framework that can account for perceptual processes, including feedforward recognition and feedback modulation. The framework explains both psychophysical and physiological experimental data and maps well onto the hierarchical distributed cortical anatomy. However, the complexity required to model cortical processes makes inference, even using approximate methods, very computationally expensive. Thus, existing object perception models based on this approach are typically limited to tree-structured networks with no loops, use small toy examples or fail to account for certain perceptual aspects such as invariance to transformations or feedback reconstruction. In this study we develop a Bayesian network with an architecture similar to that of HMAX, a biologically-inspired hierarchical model of object recognition, and use loopy belief propagation to approximate the model operations (selectivity and invariance). Crucially, the resulting Bayesian network extends the functionality of HMAX by including top-down recursive feedback. Thus, the proposed model not only achieves successful feedforward recognition invariant to noise, occlusions, and changes in position and size, but is also able to reproduce modulatory effects such as illusory contour completion and attention. Our novel and rigorous methodology covers key aspects such as learning using a layerwise greedy algorithm, combining feedback information from multiple parents and reducing the number of operations required. Overall, this work extends an established model of object recognition to include high-level feedback modulation, based on state-of-the-art probabilistic approaches. The methodology employed, consistent with evidence from the visual cortex, can be potentially generalized to build models of hierarchical perceptual organization that include top-down and bottom-up interactions, for example, in other sensory modalities.