Molecular mechanisms underlying striatal synaptic plasticity: relevance to chronic alcohol consumption and seeking

A stochastic model of hippocampal synaptic plasticity with geometrical readout of enzyme dynamics

Discovering the rules of synaptic plasticity is an important step for understanding brain learning. Existing plasticity models are either (1) top-down and interpretable, but not flexible enough to account for experimental data, or (2) bottom-up and biologically realistic, but too intricate to interpret and hard to fit to data. To avoid the shortcomings of these approaches, we present a new plasticity rule based on a geometrical readout mechanism that flexibly maps synaptic enzyme dynamics to predict plasticity outcomes. We apply this readout to a multi-timescale model of hippocampal synaptic plasticity induction that includes electrical dynamics, calcium, CaMKII and calcineurin, and accurate representation of intrinsic noise sources. Using a single set of model parameters, we demonstrate the robustness of this plasticity rule by reproducing nine published ex vivo experiments covering various spike-timing and frequency-dependent plasticity induction protocols, animal ages, and experimental conditions. Our model also predicts that in vivo-like spike timing irregularity strongly shapes plasticity outcome. This geometrical readout modelling approach can be readily applied to other excitatory or inhibitory synapses to discover their synaptic plasticity rules.

Genetic mechanisms for impaired synaptic plasticity in schizophrenia revealed by computational modelling

Schizophrenia phenotypes are suggestive of impaired cortical plasticity in the disease, but the mechanisms of these deficits are unknown. Genomic association studies have implicated a large number of genes that regulate neuromodulation and plasticity, indicating that the plasticity deficits have a genetic origin. Here, we used biochemically detailed computational modelling of post-synaptic plasticity to investigate how schizophrenia-associated genes regulate long-term potentiation (LTP) and depression (LTD). We combined our model with data from post-mortem mRNA expression studies (CommonMind gene-expression datasets) to assess the consequences of altered expression of plasticity-regulating genes for the amplitude of LTP and LTD. Our results show that the expression alterations observed post mortem, especially those in anterior cingulate cortex, lead to impaired PKA-pathway-mediated LTP in synapses containing GluR1 receptors. We validated these findings using a genotyped EEG dataset where polygenic risk scores for synaptic and ion channel-encoding genes as well as modulation of visual evoked potentials (VEP) were determined for 286 healthy controls. Our results provide a possible genetic mechanism for plasticity impairments in schizophrenia, which can lead to improved understanding and, ultimately, treatment of the disorder.

Computational analysis of cortical neuronal excitotoxicity in a large animal model of neonatal brain injury

BACKGROUND

Neonatal hypoxic brain injury is a major cause of intellectual and developmental disability. Hypoxia causes neuronal dysfunction and death in the developing cerebral cortex due to excitotoxic Ca²⁺-influx. In the translational piglet model of hypoxic encephalopathy, we have previously shown that hypoxia overactivates Ca²⁺/Calmodulin (CaM) signaling via Sarcoma (Src) kinase in cortical neurons, resulting in overexpression of proapoptotic genes. However, identifying the exact relationship between alterations in neuronal Ca²⁺-influx, molecular determinants of cell death, and the degree of hypoxia in a dynamic system represents a significant challenge.

METHODS

We used experimental and computational methods to identify molecular events critical to the onset of excitotoxicity-induced apoptosis in the cerebral cortex of newborn piglets. We used 2-3-day-old piglets (normoxic [Nx], hypoxic [Hx], and hypoxic + Src-inhibitor-treatment [Hx+PP2] groups) for biochemical analysis of ATP production, Ca²⁺-influx, and Ca²⁺/CaM-dependent protein kinase kinase 2 (CaMKK2) expression. We then used SimBiology to build a computational model of the Ca²⁺/CaM-Src-kinase signaling cascade, simulating Nx, Hx, and Hx+PP2 conditions. To evaluate our model, we used Sobol variance decomposition, multiparametric global sensitivity analysis, and parameter scanning.

RESULTS

Our model captures important molecular trends caused by hypoxia in the piglet brain. Incorporating the action of Src kinase inhibitor PP2 further validated our model and enabled predictive analysis of the effect of hypoxia on CaMKK2. We determined the impact of a feedback loop related to Src phosphorylation of NMDA receptors and activation kinetics of CaMKII. We also identified distinct modes of signaling wherein Ca²⁺ level alterations following Src kinase inhibition may not be a linear predictor of changes in Bax expression. Importantly, our model indicates that while pharmacological pre-treatment significantly reduces the onset of abnormal Ca²⁺-influx, there exists a window of intervention after hypoxia during which targeted modulation of Src-NMDAR interaction kinetics in combination with PP2 administration can reduce Ca²⁺-influx and Bax expression to similar levels as pre-treatment.

CONCLUSIONS

Our model identifies new dynamics of critical components in the Ca²⁺/CaM-Src signaling pathway leading to neuronal injury and provides a feasible framework for drug efficacy studies in translational models of neonatal brain injury for the prevention of intellectual and developmental disabilities.

A Multisubcellular Compartment Model of AMPA Receptor Trafficking for Neuromodulation of Hebbian Synaptic Plasticity

Neuromodulation can profoundly impact the gain and polarity of postsynaptic changes in Hebbian synaptic plasticity. An emerging pattern observed in multiple central synapses is a pull–push type of control in which activation of receptors coupled to the G-protein Gs promote long-term potentiation (LTP) at the expense of long-term depression (LTD), whereas receptors coupled to Gq promote LTD at the expense of LTP. Notably, coactivation of both Gs- and Gq-coupled receptors enhances the gain of both LTP and LTD. To account for these observations, we propose a simple kinetic model in which AMPA receptors (AMPARs) are trafficked between multiple subcompartments in and around the postsynaptic spine. In the model AMPARs in the postsynaptic density compartment (PSD) are the primary contributors to synaptic conductance. During LTP induction, AMPARs are trafficked to the PSD primarily from a relatively small perisynaptic (peri-PSD) compartment. Gs-coupled receptors promote LTP by replenishing peri-PSD through increased AMPAR exocytosis from a pool of endocytic AMPAR. During LTD induction AMPARs are trafficked in the reverse direction, from the PSD to the peri-PSD compartment, and Gq-coupled receptors promote LTD by clearing the peri-PSD compartment through increased AMPAR endocytosis. We claim that the model not only captures essential features of the pull–push neuromodulation of synaptic plasticity, but it is also consistent with other actions of neuromodulators observed in slice experiments and is compatible with the current understanding of AMPAR trafficking.

A computational model of dopaminergic modulation of hippocampal Schaffer collateral-CA1 long-term plasticity

Dopamine plays a critical role in modulating the long-term synaptic plasticity of the hippocampal Schaffer collateral-CA1 pyramidal neuron synapses (SC-CA1), a widely accepted cellular model of learning and memory. Limited results from hippocampal slice experiments over the last four decades have shown that the timing of the activation of dopamine D1/D5 receptors relative to a high/low-frequency stimulation (HFS/LFS) in SC-CA1 synapses regulates the modulation of HFS/LFS-induced long-term potentiation/depression (LTP/LTD) in these synapses. However, the existing literature lacks a complete picture of how various concentrations of D1/D5 agonists and the relative timing between the activation of D1/D5 receptors and LTP/LTD induction by HFS/LFS, affect the spatiotemporal modulation of SC-CA1 synaptic dynamics. In this paper, we have developed a computational model, a first of its kind, to make quantitative predictions of the temporal dose-dependent modulation of the HFS/LFS induced LTP/LTD in SC-CA1 synapses by various D1/D5 agonists. Our model combines the biochemical effects with the electrical effects at the electrophysiological level. We have estimated the model parameters from the published electrophysiological data, available from diverse hippocampal CA1 slice experiments, in a Bayesian framework. Our modeling results demonstrate the capability of our model in making quantitative predictions of the available experimental results under diverse HFS/LFS protocols. The predictions from our model show a strong nonlinear dependency of the modulated LTP/LTD by D1/D5 agonists on the relative timing between the activated D1/D5 receptors and the HFS/LFS protocol and the applied concentration of D1/D5 agonists.

Chronic and Acute Manipulation of Cortical Glutamate Transmission Induces Structural and Synaptic Changes in Co-cultured Striatal Neurons

In contrast to the prenatal topographic development of sensory cortices, striatal circuit organization is slow and requires the functional maturation of cortical and thalamic excitatory inputs throughout the first postnatal month. While mechanisms regulating synapse development and plasticity are quite well described at excitatory synapses of glutamatergic neurons in the neocortex, comparatively little is known of how this translates to glutamate synapses onto GABAergic neurons in the striatum. Here we investigate excitatory striatal synapse plasticity in an in vitro system, where glutamate can be studied in isolation from dopamine and other neuromodulators. We examined pre-and post-synaptic structural and functional plasticity in GABAergic striatal spiny projection neurons (SPNs), co-cultured with glutamatergic cortical neurons. After synapse formation, medium-term (24 h) TTX silencing increased the density of filopodia, and modestly decreased dendritic spine density, when assayed at 21 days in vitro (DIV). Spine reductions appeared to require residual spontaneous activation of ionotropic glutamate receptors. Conversely, chronic (14 days) TTX silencing markedly reduced spine density without any observed increase in filopodia density. Time-dependent, biphasic changes to the presynaptic marker Synapsin-1 were also observed, independent of residual spontaneous activity. Acute silencing (3 h) did not affect presynaptic markers or postsynaptic structures. To induce rapid, activity-dependent plasticity in striatal neurons, a chemical NMDA receptor-dependent “long-term potentiation (LTP)” paradigm was employed. Within 30 min, this increased spine and GluA1 cluster densities, and the percentage of spines containing GluA1 clusters, without altering the presynaptic signal. The results demonstrate that the growth and pruning of dendritic protrusions is an active process, requiring glutamate receptor activity in striatal projection neurons. Furthermore, NMDA receptor activation is sufficient to drive glutamatergic structural plasticity in SPNs, in the absence of dopamine or other neuromodulators.

A unified computational model for cortical post-synaptic plasticity

Signalling pathways leading to post-synaptic plasticity have been examined in many types of experimental studies, but a unified picture on how multiple biochemical pathways collectively shape neocortical plasticity is missing. We built a biochemically detailed model of post-synaptic plasticity describing CaMKII, PKA, and PKC pathways and their contribution to synaptic potentiation or depression. We developed a statistical AMPA-receptor-tetramer model, which permits the estimation of the AMPA-receptor-mediated maximal synaptic conductance based on numbers of GluR1s and GluR2s predicted by the biochemical signalling model. We show that our model reproduces neuromodulator-gated spike-timing-dependent plasticity as observed in the visual cortex and can be fit to data from many cortical areas, uncovering the biochemical contributions of the pathways pinpointed by the underlying experimental studies. Our model explains the dependence of different forms of plasticity on the availability of different proteins and can be used for the study of mental disorder-associated impairments of cortical plasticity.

Signaling models for dopamine-dependent temporal contiguity in striatal synaptic plasticity

Animals remember temporal links between their actions and subsequent rewards. We previously discovered a synaptic mechanism underlying such reward learning in D1 receptor (D1R)-expressing spiny projection neurons (D1 SPN) of the striatum. Dopamine (DA) bursts promote dendritic spine enlargement in a time window of only a few seconds after paired pre- and post-synaptic spiking (pre-post pairing), which is termed as reinforcement plasticity (RP). The previous study has also identified underlying signaling pathways; however, it still remains unclear how the signaling dynamics results in RP. In the present study, we first developed a computational model of signaling dynamics of D1 SPNs. The D1 RP model successfully reproduced experimentally observed protein kinase A (PKA) activity, including its critical time window. In this model, adenylate cyclase type 1 (AC1) in the spines/thin dendrites played a pivotal role as a coincidence detector against pre-post pairing and DA burst. In particular, pre-post pairing (Ca2+ signal) stimulated AC1 with a delay, and the Ca2+-stimulated AC1 was activated by the DA burst for the asymmetric time window. Moreover, the smallness of the spines/thin dendrites is crucial to the short time window for the PKA activity. We then developed a RP model for D2 SPNs, which also predicted the critical time window for RP that depended on the timing of pre-post pairing and phasic DA dip. AC1 worked for the coincidence detector in the D2 RP model as well. We further simulated the signaling pathway leading to Ca2+/calmodulin-dependent protein kinase II (CaMKII) activation and clarified the role of the downstream molecules of AC1 as the integrators that turn transient input signals into persistent spine enlargement. Finally, we discuss how such timing windows guide animals' reward learning.

The microcircuits of striatum in silico

The basal ganglia play an important role in decision making and selection of action primarily based on input from cortex, thalamus, and the dopamine system. Their main input structure, striatum, is central to this process. It consists of two types of projection neurons, together representing 95% of the neurons, and 5% of interneurons, among which are the cholinergic, fast-spiking, and low threshold-spiking subtypes. The membrane properties, soma-dendritic shape, and intrastriatal and extrastriatal synaptic interactions of these neurons are quite well described in the mouse, and therefore they can be simulated in sufficient detail to capture their intrinsic properties, as well as the connectivity. We focus on simulation at the striatal cellular/microcircuit level, in which the molecular/subcellular and systems levels meet. We present a nearly full-scale model of the mouse striatum using available data on synaptic connectivity, cellular morphology, and electrophysiological properties to create a microcircuit mimicking the real network. A striatal volume is populated with reconstructed neuronal morphologies with appropriate cell densities, and then we connect neurons together based on appositions between neurites as possible synapses and constrain them further with available connectivity data. Moreover, we simulate a subset of the striatum involving 10,000 neurons, with input from cortex, thalamus, and the dopamine system, as a proof of principle. Simulation at this biological scale should serve as an invaluable tool to understand the mode of operation of this complex structure. This platform will be updated with new data and expanded to simulate the entire striatum.

A Dual Role Hypothesis of the Cortico-Basal-Ganglia Pathways: Opponency and Temporal Difference Through Dopamine and Adenosine

The hypothesis that the basal-ganglia direct and indirect pathways represent goodness (or benefit) and badness (or cost) of options, respectively, explains a wide range of phenomena. However, this hypothesis, named the Opponent Actor Learning (OpAL), still has limitations. Structurally, the OpAL model does not incorporate differentiation of the two types of cortical inputs to the basal-ganglia pathways received from intratelencephalic (IT) and pyramidal-tract (PT) neurons. Functionally, the OpAL model does not describe the temporal-difference (TD)-type reward-prediction-error (RPE), nor explains how RPE is calculated in the circuitry connecting to the DA neurons. In fact, there is a different hypothesis on the basal-ganglia pathways and DA, named the Cortico-Striatal-Temporal-Difference (CS-TD) model. The CS-TD model differentiates the IT and PT inputs, describes the TD-type RPE, and explains how TD-RPE is calculated. However, a critical difficulty in this model lies in its assumption that DA induces the same direction of plasticity in both direct and indirect pathways, which apparently contradicts the experimentally observed opposite effects of DA on these pathways. Here, we propose a new hypothesis that integrates the OpAL and CS-TD models. Specifically, we propose that the IT-basal-ganglia pathways represent goodness/badness of current options while the PT-indirect pathway represents the overall value of the previously chosen option, and both of these have influence on the DA neurons, through the basal-ganglia output, so that a variant of TD-RPE is calculated. A key assumption is that opposite directions of plasticity are induced upon phasic activation of DA neurons in the IT-indirect pathway and PT-indirect pathway because of different profiles of IT and PT inputs. Specifically, at PT→indirect-pathway-medium-spiny-neuron (iMSN) synapses, sustained glutamatergic inputs generate rich adenosine, which allosterically prevents DA-D2 receptor signaling and instead favors adenosine-A2A receptor signaling. Then, phasic DA-induced phasic adenosine, which reflects TD-RPE, causes long-term synaptic potentiation. In contrast, at IT→iMSN synapses where adenosine is scarce, phasic DA causes long-term synaptic depression via D2 receptor signaling. This new Opponency and Temporal-Difference (OTD) model provides unique predictions, part of which is potentially in line with recently reported activity patterns of neurons in the globus pallidus externus on the indirect pathway.

Using NEURON for Reaction-Diffusion Modeling of Extracellular Dynamics

Development of credible clinically-relevant brain simulations has been slowed due to a focus on electrophysiology in computational neuroscience, neglecting the multiscale whole-tissue modeling approach used for simulation in most other organ systems. We have now begun to extend the NEURON simulation platform in this direction by adding extracellular modeling. The extracellular medium of neural tissue is an active medium of neuromodulators, ions, inflammatory cells, oxygen, NO and other gases, with additional physiological, pharmacological and pathological agents. These extracellular agents influence, and are influenced by, cellular electrophysiology, and cellular chemophysiology-the complex internal cellular milieu of second-messenger signaling and cascades. NEURON's extracellular reaction-diffusion is supported by an intuitive Python-based where/who/what command sequence, derived from that used for intracellular reaction diffusion, to support coarse-grained macroscopic extracellular models. This simulation specification separates the expression of the conceptual model and parameters from the underlying numerical methods. In the volume-averaging approach used, the macroscopic model of tissue is characterized by free volume fraction-the proportion of space in which species are able to diffuse, and tortuosity-the average increase in path length due to obstacles. These tissue characteristics can be defined within particular spatial regions, enabling the modeler to account for regional differences, due either to intrinsic organization, particularly gray vs. white matter, or to pathology such as edema. We illustrate simulation development using spreading depression, a pathological phenomenon thought to play roles in migraine, epilepsy and stroke. Simulation results were verified against analytic results and against the extracellular portion of the simulation run under FiPy. The creation of this NEURON interface provides a pathway for interoperability that can be used to automatically export this class of models into complex intracellular/extracellular simulations and future cross-simulator standardization.