Cortical information flow in Parkinson's disease: a composite network/field model

Sensitivity-analysis-guided Bayesian parameter estimation for neural mass models: Applications in epilepsy

It is well established that neural mass models (NMMs) can effectively simulate the mesoscopic and macroscopic dynamics of electroencephalography (EEG), including epileptic EEG. However, the use of NMMs to gain insight on the neuronal system by parameter estimation is hampered by their high dimensionality and the lack of knowledge on what NMM parameters can be reliably estimated. In this article, we analyze the parameter sensitivity of the Jansen and Rit NMM (JR NMM) in order to identify the most sensitive JR-NMM parameters for reliable parameter estimation from EEG data. We then propose a Bayesian approach for estimating the JR-NMM states and parameters based on an expectation-maximization algorithm combined with the unscented Kalman smoother (UKS EM). Global sensitivity analysis methods including the Morris method and the Sobol method are used to perform sensitivity analysis. Results from both the Morris and the Sobol method show that the average inhibitory synaptic gain, B, and the reciprocal of the time constant of the average inhibitory postsynaptic potentials, b, have a significant impact on the JR-NMM output along with having the least interaction with other model parameters. The UKS-EM method for estimating the parameters B and b is validated using simulations under varying levels of measurement noise. Finally we apply the UKS-EM algorithm to intracranial EEG data from 16 epileptic patients. Our results, both at individual and group level show that the parameters B and b change significantly between the preseizure and seizure period, and between the seizure and postseizure period, with the transition to seizure characterized by a decrease in the average B, and the high frequency activity in seizure characterized by an increase in b. These results establish a sensitivity analysis guided Bayesian parameter estimation as a powerful tool for reducing the parameter space of high-dimensional NMMs enabling reliable and efficient estimation of the most sensitive NMM parameters, with the potential for online and fast tracking of NMM parameters in applications such as seizure tracking and control.

A Model with Dopamine Depletion in Basal Ganglia and Cerebellum Predicts Changes in Thalamocortical Beta Oscillations

Parkinsonism is presented as a motor syndrome characterized by rigidity, tremors, and bradykinesia, with Parkinson's disease (PD) being the predominant cause. The discovery that those motor symptoms result from the death of dopaminergic cells in the substantia nigra led to focus most of parkinsonism research on the basal ganglia (BG). However, recent findings point to an active involvement of the cerebellum in this motor syndrome. Here, we have developed a multiscale computational model of the rodent brain's BG-cerebellar network. Simulations showed that a direct effect of dopamine depletion on the cerebellum must be taken into account to reproduce the alterations of neural activity in parkinsonism, particularly the increased beta oscillations widely reported in PD patients. Moreover, dopamine depletion indirectly impacted spike-time-dependent plasticity at the parallel fiber-Purkinje cell synapses, degrading associative motor learning as observed in parkinsonism. Overall, these results suggest a relevant involvement of cerebellum in parkinsonism associative motor symptoms.

Dysfunction of motor cortices in Parkinson's disease

The cerebral cortex has long been thought to be involved in the pathophysiology of motor symptoms of Parkinson's disease. The impaired cortical function is believed to be a direct and immediate effect of pathologically patterned basal ganglia output, mediated to the cerebral cortex by way of the ventral motor thalamus. However, recent studies in humans with Parkinson's disease and in animal models of the disease have provided strong evidence suggesting that the involvement of the cerebral cortex is much broader than merely serving as a passive conduit for subcortical disturbances. In the present review, we discuss Parkinson's disease-related changes in frontal cortical motor regions, focusing on neuropathology, plasticity, changes in neurotransmission, and altered network interactions. We will also examine recent studies exploring the cortical circuits as potential targets for neuromodulation to treat Parkinson's disease.

Ultrafast simulation of large-scale neocortical microcircuitry with biophysically realistic neurons

Understanding the activity of the mammalian brain requires an integrative knowledge of circuits at distinct scales, ranging from ion channel gating to circuit connectomics. Computational models are regularly employed to understand how multiple parameters contribute synergistically to circuit behavior. However, traditional models of anatomically and biophysically realistic neurons are computationally demanding, especially when scaled to model local circuits. To overcome this limitation, we trained several artificial neural network (ANN) architectures to model the activity of realistic multicompartmental cortical neurons. We identified an ANN architecture that accurately predicted subthreshold activity and action potential firing. The ANN could correctly generalize to previously unobserved synaptic input, including in models containing nonlinear dendritic properties. When scaled, processing times were orders of magnitude faster compared with traditional approaches, allowing for rapid parameter-space mapping in a circuit model of Rett syndrome. Thus, we present a novel ANN approach allowing for rapid, detailed network experiments using inexpensive and commonly available computational resources.

Neuro4PD: An Initial Neurorobotics Model of Parkinson's Disease

In this work, we present the first steps toward the creation of a new neurorobotics model of Parkinson's Disease (PD) that embeds, for the first time in a real robot, a well-established computational model of PD. PD mostly affects the modulation of movement in humans. The number of people suffering from this neurodegenerative disease is set to double in the next 15 years and there is still no cure. With the new model we were capable to further explore the dynamics of the disease using a humanoid robot. Results show that the embedded model under both conditions, healthy and parkinsonian, was capable of performing a simple behavioural task with different levels of motor disturbance. We believe that this neurorobotics model is a stepping stone to the development of more sophisticated models that could eventually test and inform new PD therapies and help to reduce and replace animals in research.

Dysregulation of excitatory neural firing replicates physiological and functional changes in aging visual cortex

The mammalian visual system has been the focus of countless experimental and theoretical studies designed to elucidate principles of neural computation and sensory coding. Most theoretical work has focused on networks intended to reflect developing or mature neural circuitry, in both health and disease. Few computational studies have attempted to model changes that occur in neural circuitry as an organism ages non-pathologically. In this work we contribute to closing this gap, studying how physiological changes correlated with advanced age impact the computational performance of a spiking network model of primary visual cortex (V1). Our results demonstrate that deterioration of homeostatic regulation of excitatory firing, coupled with long-term synaptic plasticity, is a sufficient mechanism to reproduce features of observed physiological and functional changes in neural activity data, specifically declines in inhibition and in selectivity to oriented stimuli. This suggests a potential causality between dysregulation of neuron firing and age-induced changes in brain physiology and functional performance. While this does not rule out deeper underlying causes or other mechanisms that could give rise to these changes, our approach opens new avenues for exploring these underlying mechanisms in greater depth and making predictions for future experiments.

Evoked response activity eigenmode analysis in a convoluted cortex via neural field theory

Neural field theory of the corticothalamic system is used to explore evoked response potentials (ERPs) caused by spatially localized impulse stimuli on the convoluted cortex and on a spherical cortex. Eigenfunctions are calculated analytically on the spherical cortex and numerically on the convoluted cortex via eigenfunction expansions. Eigenmodes on a convoluted cortex are similar to those of the spherical cortex, and a few such modes are found to be sufficient to reproduce the main ERP features. It is found that the ERP peak is stronger in spherical cortex than convoluted cortex, but in both cases the peak decreases monotonically with increasing distance from the stimulus point. In the convoluted case, cortical folding causes ERPs to differ between locations at the same distance from the stimulus point and spherical symmetries are only approximately preserved.

A review of computational modeling and deep brain stimulation: applications to Parkinson's disease

Biophysical computational models are complementary to experiments and theories, providing powerful tools for the study of neurological diseases. The focus of this review is the dynamic modeling and control strategies of Parkinson's disease (PD). In previous studies, the development of parkinsonian network dynamics modeling has made great progress. Modeling mainly focuses on the cortex-thalamus-basal ganglia (CTBG) circuit and its sub-circuits, which helps to explore the dynamic behavior of the parkinsonian network, such as synchronization. Deep brain stimulation (DBS) is an effective strategy for the treatment of PD. At present, many studies are based on the side effects of the DBS. However, the translation from modeling results to clinical disease mitigation therapy still faces huge challenges. Here, we introduce the progress of DBS improvement. Its specific purpose is to develop novel DBS treatment methods, optimize the treatment effect of DBS for each patient, and focus on the study in closed-loop DBS. Our goal is to review the inspiration and insights gained by combining the system theory with these computational models to analyze neurodynamics and optimize DBS treatment.

Neural field theory of evoked response potentials in a spherical brain geometry

Evoked response potentials (ERPs) are calculated in spherical and planar geometries using neural field theory of the corticothalamic system. The ERP is modeled as an impulse response and the resulting modal effects of spherical corticothalamic dynamics are explored, showing that results for spherical and planar geometries converge in the limit of large brain size. Cortical modal effects can lead to a double-peak structure in the ERP time series. It is found that the main difference between infinite planar geometry and spherical geometry is that the ERP peak is sharper and stronger in the spherical geometry. It is also found that the magnitude of the response decreases with increasing spatial width of the stimulus at the cortex. The peak is slightly delayed at large angles from the stimulus point, corresponding to group velocities of 6-10 m s^{-1}. Strong modal effects are found in the spherical geometry, with the lowest few modes sufficing to describe the main features of ERPs, except very near to spatially narrow stimuli.

Physiology-based ERPs in normal and abnormal states

Evoked response potentials (ERPs) and other transients are modeled as impulse responses using physiology-based neural field theory (NFT) of the corticothalamic system of neural activity in the human brain that incorporates synaptic and dendritic dynamics, firing response, axonal propagation, and corticocortical and corticothalamic pathways. The properties of model-predicted ERPs are explored throughout the stability zone of the corticothalamic system, and predicted time series and wavelet spectra are also analyzed. This provides a unified treatment of predicted ERPs for both normal and abnormal states within the brain's stability zone, including likely parameters to represent abnormal states of reduced arousal.

NFTsim: Theory and Simulation of Multiscale Neural Field Dynamics

A user ready, portable, documented software package, NFTsim, is presented to facilitate numerical simulations of a wide range of brain systems using continuum neural field modeling. NFTsim enables users to simulate key aspects of brain activity at multiple scales. At the microscopic scale, it incorporates characteristics of local interactions between cells, neurotransmitter effects, synaptodendritic delays and feedbacks. At the mesoscopic scale, it incorporates information about medium to large scale axonal ranges of fibers, which are essential to model dissipative wave transmission and to produce synchronous oscillations and associated cross-correlation patterns as observed in local field potential recordings of active tissue. At the scale of the whole brain, NFTsim allows for the inclusion of long range pathways, such as thalamocortical projections, when generating macroscopic activity fields. The multiscale nature of the neural activity produced by NFTsim has the potential to enable the modeling of resulting quantities measurable via various neuroimaging techniques. In this work, we give a comprehensive description of the design and implementation of the software. Due to its modularity and flexibility, NFTsim enables the systematic study of an unlimited number of neural systems with multiple neural populations under a unified framework and allows for direct comparison with analytic and experimental predictions. The code is written in C++ and bundled with Matlab routines for a rapid quantitative analysis and visualization of the outputs. The output of NFTsim is stored in plain text file enabling users to select from a broad range of tools for offline analysis. This software enables a wide and convenient use of powerful physiologically-based neural field approaches to brain modeling. NFTsim is distributed under the Apache 2.0 license.

Credibility, Replicability, and Reproducibility in Simulation for Biomedicine and Clinical Applications in Neuroscience

Modeling and simulation in computational neuroscience is currently a research enterprise to better understand neural systems. It is not yet directly applicable to the problems of patients with brain disease. To be used for clinical applications, there must not only be considerable progress in the field but also a concerted effort to use best practices in order to demonstrate model credibility to regulatory bodies, to clinics and hospitals, to doctors, and to patients. In doing this for neuroscience, we can learn lessons from long-standing practices in other areas of simulation (aircraft, computer chips), from software engineering, and from other biomedical disciplines. In this manuscript, we introduce some basic concepts that will be important in the development of credible clinical neuroscience models: reproducibility and replicability; verification and validation; model configuration; and procedures and processes for credible mechanistic multiscale modeling. We also discuss how garnering strong community involvement can promote model credibility. Finally, in addition to direct usage with patients, we note the potential for simulation usage in the area of Simulation-Based Medical Education, an area which to date has been primarily reliant on physical models (mannequins) and scenario-based simulations rather than on numerical simulations.

Optimization by Adaptive Stochastic Descent

When standard optimization methods fail to find a satisfactory solution for a parameter fitting problem, a tempting recourse is to adjust parameters manually. While tedious, this approach can be surprisingly powerful in terms of achieving optimal or near-optimal solutions. This paper outlines an optimization algorithm, Adaptive Stochastic Descent (ASD), that has been designed to replicate the essential aspects of manual parameter fitting in an automated way. Specifically, ASD uses simple principles to form probabilistic assumptions about (a) which parameters have the greatest effect on the objective function, and (b) optimal step sizes for each parameter. We show that for a certain class of optimization problems (namely, those with a moderate to large number of scalar parameter dimensions, especially if some dimensions are more important than others), ASD is capable of minimizing the objective function with far fewer function evaluations than classic optimization methods, such as the Nelder-Mead nonlinear simplex, Levenberg-Marquardt gradient descent, simulated annealing, and genetic algorithms. As a case study, we show that ASD outperforms standard algorithms when used to determine how resources should be allocated in order to minimize new HIV infections in Swaziland.

Computational modelling of pathogenic protein spread in neurodegenerative diseases

Pathogenic protein accumulation and spread are fundamental principles of neurodegenerative diseases and ultimately account for the atrophy patterns that distinguish these diseases clinically. However, the biological mechanisms that link pathogenic proteins to specific neural network damage patterns have not been defined. We developed computational models for mechanisms of pathogenic protein accumulation, spread and toxic effects in an artificial neural network of cortical columns. By varying simulation parameters we assessed the effects of modelled mechanisms on network breakdown patterns. Our findings suggest that patterns of network breakdown and the convergence of patterns follow rules determined by particular protein parameters. These rules can account for empirical data on pathogenic protein spread in neural networks. This work provides a basis for understanding the effects of pathogenic proteins on neural circuits and predicting progression of neurodegeneration.

The role of cortical oscillations in a spiking neural network model of the basal ganglia

Although brain oscillations involving the basal ganglia (BG) have been the target of extensive research, the main focus lies disproportionally on oscillations generated within the BG circuit rather than other sources, such as cortical areas. We remedy this here by investigating the influence of various cortical frequency bands on the intrinsic effective connectivity of the BG, as well as the role of the latter in regulating cortical behaviour. To do this, we construct a detailed neural model of the complete BG circuit based on fine-tuned spiking neurons, with both electrical and chemical synapses as well as short-term plasticity between structures. As a measure of effective connectivity, we estimate information transfer between nuclei by means of transfer entropy. Our model successfully reproduces firing and oscillatory behaviour found in both the healthy and Parkinsonian BG. We found that, indeed, effective connectivity changes dramatically for different cortical frequency bands and phase offsets, which are able to modulate (or even block) information flow in the three major BG pathways. In particular, alpha (8–12Hz) and beta (13–30Hz) oscillations activate the direct BG pathway, and favour the modulation of the indirect and hyper-direct pathways via the subthalamic nucleus—globus pallidus loop. In contrast, gamma (30–90Hz) frequencies block the information flow from the cortex completely through activation of the indirect pathway. Finally, below alpha, all pathways decay gradually and the system gives rise to spontaneous activity generated in the globus pallidus. Our results indicate the existence of a multimodal gating mechanism at the level of the BG that can be entirely controlled by cortical oscillations, and provide evidence for the hypothesis of cortically-entrained but locally-generated subthalamic beta activity. These two findings suggest new insights into the pathophysiology of specific BG disorders.

Computational models of basal ganglia dysfunction: the dynamics is in the details

The development, simulation, and analysis of mathematical models offer helpful tools for integrating experimental findings and exploring or suggesting possible explanatory mechanisms. As models relating to basal ganglia dysfunction have proliferated, however, there has not always been consistency among their findings. This work points out several ways in which biological details, relating to ionic currents and synaptic pathways, can influence the dynamics of models of the basal ganglia under parkinsonian conditions and hence may be important for inclusion in models. It also suggests some additional useful directions for future modeling studies relating to basal ganglia dysfunction.

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.

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.

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.

Functional correlates of exaggerated oscillatory activity in basal ganglia output in hemiparkinsonian rats

Exaggerated beta range (13-30Hz) synchronized activity is observed in the basal ganglia of Parkinson's disease (PD) patients during implantation of deep brain stimulation electrodes and is thought to contribute to the motor symptoms of this disorder. To explore the translational potential of similar activity observed in a rat model of PD, local field potentials (LFPs) and spiking activity in basal ganglia output were characterized in rats with unilateral dopamine cell lesion during a range of behaviors. A circular treadmill was used to assess activity during walking; hemiparkinsonian rats could maintain a steady gait when oriented ipsiversive to the lesioned hemisphere, but were less effective at walking when oriented contraversive to lesion. Dramatic increases in substantia nigra pars reticulata (SNpr) LFP oscillatory activity and spike-LFP synchronization were observed within the beta/low gamma range (12-40Hz) in the lesioned hemisphere, relative to the non-lesioned hemisphere, with the dominant frequency of spike-LFP entrainment and LFP power varying with behavioral state. At 3weeks postlesion, the mean dominant entrainment frequency during ipsiversive treadmill walking and grooming was 34Hz. Other behaviors were associated with lower mean entrainment frequencies: 27-28Hz during alert non-walking and REM, 17Hz during rest and 21Hz during urethane anesthesia with sensory stimulation. SNpr spike-LFP entrainment frequency was stable during individual treadmill walking epochs, but increased gradually over weeks postlesion. In contrast, SNpr LFP power in the 25-40Hz range was greatest at the initiation of each walking epoch, and decreased during walking to stabilize by 6min at 49% of initial values. Power was further modulated in conjunction with the 1.5s stepping rhythm. Administration of l-dopa improved contraversive treadmill walking in correlation with a reduction in SNpr 25-40Hz LFP power and spike synchronization in the dopamine cell lesioned hemisphere. These effects were reversed by the serotonergic 1A agonist, 8-OH-DPAT. While the prominent spike-LFP phase locking observed during ongoing motor activity in the hemiparkinsonian rats occurs at frequencies intriguingly higher than in PD patients, the synchronized activity in the SNpr of this animal model has much in common with oscillatory activity recorded from the basal ganglia of the PD patients. Results support the potential of this model for providing insight into relationships between synchronization of basal ganglia output induced by loss of dopamine and motor symptoms in PD.

Multiscale modeling for clinical translation in neuropsychiatric disease

Multiscale modeling of neuropsychiatric illness bridges scales of clinical importance: from the highest scales (presentation of behavioral signs and symptoms), through intermediate scales (clinical testing and surgical intervention), down to the molecular scale of pharmacotherapy. Modeling of brain disease is difficult compared to modeling of other organs, because dysfunction manifests at scales where measurements are rudimentary due both to inadequate access (memory and cognition) and to complexity (behavior). Nonetheless, we can begin to explore these aspects through the use of information-theoretic measures as stand-ins for meaning at the top scales. We here describe efforts across five disorders: Parkinson’s, Alzheimer’s, stroke, schizophrenia, and epilepsy. We look at the use of therapeutic brain stimulation to replace lost neural signals, a loss that produces diaschisis, defined as activity changes in other brain areas due to missing inputs. These changes may in some cases be compensatory, hence beneficial, but in many cases a primary pathology, whether itself static or dynamic, sets in motion a series of dynamic consequences that produce further pathology. The simulations presented here suggest how diaschisis can be reversed by using a neuroprosthetic signal. Despite having none of the information content of the lost physiological signal, the simplified neuroprosthetic signal can restore a diaschitic area to near-normal patterns of activity. Computer simulation thus begins to explain the remarkable success of stimulation technologies - deep brain stimulation, transcranial magnetic stimulation, ultrasound stimulation, transcranial direct current stimulation - across an extremely broad range of pathologies. Multiscale modeling can help us to optimize and integrate these neuroprosthetic therapies by taking into consideration effects of different stimulation protocols, combinations of stimulation with neuropharmacological therapy, and interplay of these therapeutic modalities with particular patterns of disease focality, dynamics, and prior therapies.

Closed-loop cortical neuromodulation in Parkinson's disease: An alternative to deep brain stimulation?

Deep brain stimulation (DBS) is usually performed to treat advanced Parkinson's disease (PD) patients with electrodes permanently implanted in basal ganglia while the stimulator delivers electrical impulses continuously and independently of any feedback (open-loop stimulation). Conversely, in closed-loop stimulation, electrical stimulation is delivered as a function of neuronal activities recorded and analyzed online. There is an emerging development of closed-loop DBS in the treatment of PD and a growing discussion about proposing cortical stimulation rather than DBS for this purpose. Why does it make sense to "close the loop" to treat parkinsonian symptoms? Could closed-loop stimulation applied to the cortex become a valuable therapeutic strategy for PD? Can mathematical modeling contribute to the development of this technique? We review the various evidences in favor of the use of closed-loop cortical stimulation for the treatment of advanced PD, as an emerging technique which might offer substantial clinical benefits for PD patients.

Reaction-diffusion in the NEURON simulator

In order to support research on the role of cell biological principles (genomics, proteomics, signaling cascades and reaction dynamics) on the dynamics of neuronal response in health and disease, NEURON's Reaction-Diffusion (rxd) module in Python provides specification and simulation for these dynamics, coupled with the electrophysiological dynamics of the cell membrane. Arithmetic operations on species and parameters are overloaded, allowing arbitrary reaction formulas to be specified using Python syntax. These expressions are then transparently compiled into bytecode that uses NumPy for fast vectorized calculations. At each time step, rxd combines NEURON's integrators with SciPy's sparse linear algebra library.

Reinforcement learning of two-joint virtual arm reaching in a computer model of sensorimotor cortex

Neocortical mechanisms of learning sensorimotor control involve a complex series of interactions at multiple levels, from synaptic mechanisms to cellular dynamics to network connectomics. We developed a model of sensory and motor neocortex consisting of 704 spiking model neurons. Sensory and motor populations included excitatory cells and two types of interneurons. Neurons were interconnected with AMPA/NMDA and GABAA synapses. We trained our model using spike-timing-dependent reinforcement learning to control a two-joint virtual arm to reach to a fixed target. For each of 125 trained networks, we used 200 training sessions, each involving 15 s reaches to the target from 16 starting positions. Learning altered network dynamics, with enhancements to neuronal synchrony and behaviorally relevant information flow between neurons. After learning, networks demonstrated retention of behaviorally relevant memories by using proprioceptive information to perform reach-to-target from multiple starting positions. Networks dynamically controlled which joint rotations to use to reach a target, depending on current arm position. Learning-dependent network reorganization was evident in both sensory and motor populations: learned synaptic weights showed target-specific patterning optimized for particular reach movements. Our model embodies an integrative hypothesis of sensorimotor cortical learning that could be used to interpret future electrophysiological data recorded in vivo from sensorimotor learning experiments. We used our model to make the following predictions: learning enhances synchrony in neuronal populations and behaviorally relevant information flow across neuronal populations, enhanced sensory processing aids task-relevant motor performance and the relative ease of a particular movement in vivo depends on the amount of sensory information required to complete the movement.