Postsynaptic signal transduction models for long-term potentiation and depression

A calcium-based plasticity model for predicting long-term potentiation and depression in the neocortex

Pyramidal cells (PCs) form the backbone of the layered structure of the neocortex, and plasticity of their synapses is thought to underlie learning in the brain. However, such long-term synaptic changes have been experimentally characterized between only a few types of PCs, posing a significant barrier for studying neocortical learning mechanisms. Here we introduce a model of synaptic plasticity based on data-constrained postsynaptic calcium dynamics, and show in a neocortical microcircuit model that a single parameter set is sufficient to unify the available experimental findings on long-term potentiation (LTP) and long-term depression (LTD) of PC connections. In particular, we find that the diverse plasticity outcomes across the different PC types can be explained by cell-type-specific synaptic physiology, cell morphology and innervation patterns, without requiring type-specific plasticity. Generalizing the model to in vivo extracellular calcium concentrations, we predict qualitatively different plasticity dynamics from those observed in vitro. This work provides a first comprehensive null model for LTP/LTD between neocortical PC types in vivo, and an open framework for further developing models of cortical synaptic plasticity.

Stable continual learning through structured multiscale plasticity manifolds

Biological plasticity is ubiquitous. How does the brain navigate this complex plasticity space, where any component can seemingly change, in adapting to an ever-changing environment? We build a systematic case that stable continuous learning is achieved by structured rules that enforce multiple, but not all, components to change together in specific directions. This rule-based low-dimensional plasticity manifold of permitted plasticity combinations emerges from cell type-specific molecular signaling and triggers cascading impacts that span multiple scales. These multiscale plasticity manifolds form the basis for behavioral learning and are dynamic entities that are altered by neuromodulation, metaplasticity, and pathology. We explore the strong links between heterogeneities, degeneracy, and plasticity manifolds and emphasize the need to incorporate plasticity manifolds into learning-theoretical frameworks and experimental designs.

Online Learning and Memory of Neural Trajectory Replays for Prefrontal Persistent and Dynamic Representations in the Irregular Asynchronous State

In the prefrontal cortex (PFC), higher-order cognitive functions and adaptive flexible behaviors rely on continuous dynamical sequences of spiking activity that constitute neural trajectories in the state space of activity. Neural trajectories subserve diverse representations, from explicit mappings in physical spaces to generalized mappings in the task space, and up to complex abstract transformations such as working memory, decision-making and behavioral planning. Computational models have separately assessed learning and replay of neural trajectories, often using unrealistic learning rules or decoupling simulations for learning from replay. Hence, the question remains open of how neural trajectories are learned, memorized and replayed online, with permanently acting biological plasticity rules. The asynchronous irregular regime characterizing cortical dynamics in awake conditions exerts a major source of disorder that may jeopardize plasticity and replay of locally ordered activity. Here, we show that a recurrent model of local PFC circuitry endowed with realistic synaptic spike timing-dependent plasticity and scaling processes can learn, memorize and replay large-size neural trajectories online under asynchronous irregular dynamics, at regular or fast (sped-up) timescale. Presented trajectories are quickly learned (within seconds) as synaptic engrams in the network, and the model is able to chunk overlapping trajectories presented separately. These trajectory engrams last long-term (dozen hours) and trajectory replays can be triggered over an hour. In turn, we show the conditions under which trajectory engrams and replays preserve asynchronous irregular dynamics in the network. Functionally, spiking activity during trajectory replays at regular timescale accounts for both dynamical coding with temporal tuning in individual neurons, persistent activity at the population level, and large levels of variability consistent with observed cognitive-related PFC dynamics. Together, these results offer a consistent theoretical framework accounting for how neural trajectories can be learned, memorized and replayed in PFC networks circuits to subserve flexible dynamic representations and adaptive behaviors.

Astrocyte-mediated spike-timing-dependent long-term depression modulates synaptic properties in the developing cortex

Astrocytes have been shown to modulate synaptic transmission and plasticity in specific cortical synapses, but our understanding of the underlying molecular and cellular mechanisms remains limited. Here we present a new biophysicochemical model of a somatosensory cortical layer 4 to layer 2/3 synapse to study the role of astrocytes in spike-timing-dependent long-term depression (t-LTD) in vivo. By applying the synapse model and electrophysiological data recorded from rodent somatosensory cortex, we show that a signal from a postsynaptic neuron, orchestrated by endocannabinoids, astrocytic calcium signaling, and presynaptic N-methyl-D-aspartate receptors coupled with calcineurin signaling, induces t-LTD which is sensitive to the temporal difference between post- and presynaptic firing. We predict for the first time the dynamics of astrocyte-mediated molecular mechanisms underlying t-LTD and link complex biochemical networks at presynaptic, postsynaptic, and astrocytic sites to the time window of t-LTD induction. During t-LTD a single astrocyte acts as a delay factor for fast neuronal activity and integrates fast neuronal sensory processing with slow non-neuronal processing to modulate synaptic properties in the brain. Our results suggest that astrocytes play a critical role in synaptic computation during postnatal development and are of paramount importance in guiding the development of brain circuit functions, learning and memory.

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.

Effect of NMDAR-NMNAT1/2 pathway on neuronal cell damage and cognitive impairment of sevoflurane-induced aged rats

Objective: The possible effect of NMDAR (N-methyl-D-aspartate receptor)-NMNAT1/2 (nicotinamide/nicotinic acid mono-nucleotide adenylyltransferase) signaling pathway on the neuronal cell damage and cognitive impairment of aged rats anesthetized by sevoflurane was explored.Methods: Adult male Wistar rats were selected and divided into Control, Sevo (Sevoflurane), Sevo+DCS (NMDAR agonist D-cycloserine) 30 mg/kg, Sevo+DCS 100 mg/kg, and Sevo+DCS 200 mg/kg groups. Morris water maze and fear conditioning text were used to observe cognitive function changes of rats. The inflammatory cytokines were determined by quantitative real-time polymerase chain reaction (qRT-PCR) and enzyme-linked immunosorbent assay (ELISA) assay, neuronal apoptosis by terminal deoxynucleotidyl transferase (TdT)-mediated dUTP nick-end labelling (TUNEL) staining and MDAR-NMNAT1/2 pathway-related proteins by Western blotting.Results: The longer escape latency, decreased platform crossing times and reduced staying time spent in platform quadrant were found in rats from Sevo group, with decreased percentage of freezing time in contextual test and tone cued test; and meanwhile, these rats had increased inflammatory cytokines (interleukin (IL)-1β, tumor necrosis factor (TNF-α), IL-6, and IL-8) and neuronal apoptosis, but declined expressions of MDAR-NMNAT1/2 pathway-related proteins. However, the above changes were exhibited an opposite tendency in those Sevo rats treated with different concentrations of DCS (including 30, 100, and 200 mg/kg, respectively). Particularly, the improving effect of low-dose DCS on each aspect in aged rats was better than high-dose ones.Conclusion: Activation of NMDAR-NMNAT1/2 signaling pathway could not only reduce neuronal apoptosis, but also alleviate sevoflurane-induced neuronal inflammation and cognitive impairment in aged rats.

Degeneracy in hippocampal physiology and plasticity

Degeneracy, defined as the ability of structurally disparate elements to perform analogous function, has largely been assessed from the perspective of maintaining robustness of physiology or plasticity. How does the framework of degeneracy assimilate into an encoding system where the ability to change is an essential ingredient for storing new incoming information? Could degeneracy maintain the balance between the apparently contradictory goals of the need to change for encoding and the need to resist change towards maintaining homeostasis? In this review, we explore these fundamental questions with the mammalian hippocampus as an example encoding system. We systematically catalog lines of evidence, spanning multiple scales of analysis that point to the expression of degeneracy in hippocampal physiology and plasticity. We assess the potential of degeneracy as a framework to achieve the conjoint goals of encoding and homeostasis without cross-interferences. We postulate that biological complexity, involving interactions among the numerous parameters spanning different scales of analysis, could establish disparate routes towards accomplishing these conjoint goals. These disparate routes then provide several degrees of freedom to the encoding-homeostasis system in accomplishing its tasks in an input- and state-dependent manner. Finally, the expression of degeneracy spanning multiple scales offers an ideal reconciliation to several outstanding controversies, through the recognition that the seemingly contradictory disparate observations are merely alternate routes that the system might recruit towards accomplishment of its goals.

Multi-scale modeling of the circadian modulation of learning and memory

We propose a multi-scale model to explain the time-of-day effects on learning and memory. We specifically model the circadian variation of hippocampus (HC) dependent long-term potentiation (LTP), depression (LTD), and the fear conditioning paradigm in amygdala. The model we built has both Goodwin type circadian gene regulatory network (GRN) and the conductance model of Morris-Lecar (ML) type to explain the spontaneous firing patterns (SFR) in suprachiasmatic nucleus (SCN). In the conductance model, we also include N-Methyl-D-aspartic acid receptor (NMDAR) to study the circadian dependent changes in LTP/LTD in hippocampus and include both NMDAR and α-amino-3-hydroxy-5-methyl-4-isoxazolepropionic acid receptor (AMPAR) dynamics to explain the circadian modulation of fear conditioning paradigm in memory acquisition, recall, and extinction as seen in amygdala. Our multi-scale model captures the essential dynamics seen in the experiments and strongly supports the circadian time-of-the-day effects on learning and memory.

Postsynaptic Stability and Variability Described by a Stochastic Model of Endosomal Trafficking

Neurons undergo dynamic processes of constitutive AMPA-type glutamate receptor (AMPAR) trafficking, such as the insertion and internalization of AMPARs by exocytosis and endocytosis, while stably maintaining synaptic efficacy. Studies using advanced imaging techniques have suggested that the frequency of these constitutive trafficking processes, as well as the number of AMPARs that are involved in a particular event highly fluctuate. In addition, mechanisms that trigger some forms of synaptic plasticity have been shown to include not only these processes but also additional fluctuating processes, such as the sorting of AMPARs to late endosomes (LEs). Thus, the regulation of postsynaptic AMPARs by the endosomal trafficking system appears to have superficially conflicting properties between the stability or organized control of plasticity and highly fluctuating or stochastic processes. However, it is not clear how the endosomal trafficking system reconciles and utilizes such conflicting properties. Although deterministic models have been effective to describe the stable maintenance of synaptic AMPAR numbers by constitutive recycling, as well as the involvement of endosomal trafficking in synaptic plasticity, they do not take stochasticity into account. In this study, we introduced the stochasticity into the model of each crucial machinery of the endosomal trafficking system. The specific questions we solved by our improved model are whether stability is accomplished even with a combination of fluctuating processes, and how overall variability occurs while controlling long-term synaptic depression (LTD). Our new stochastic model indeed demonstrated the stable regulation of postsynaptic AMPAR numbers at the basal state and during LTD maintenance, despite fast fluctuations in AMPAR numbers as well as high variability in the time course and amounts of LTD. In addition, our analysis suggested that the high variability arising from this stochasticity is beneficial for reproducing the relatively constant timing of LE sorting for LTD. We therefore propose that the coexistence of stability and stochasticity in the endosomal trafficking system is suitable for stable synaptic transmission and the reliable induction of synaptic plasticity, with variable properties that have been observed experimentally.

FindSim: A Framework for Integrating Neuronal Data and Signaling Models

Current experiments touch only small but overlapping parts of very complex subcellular signaling networks in neurons. Even with modern optical reporters and pharmacological manipulations, a given experiment can only monitor and control a very small subset of the diverse, multiscale processes of neuronal signaling. We have developed FindSim (Framework for Integrating Neuronal Data and SIgnaling Models) to anchor models to structured experimental datasets. FindSim is a framework for integrating many individual electrophysiological and biochemical experiments with large, multiscale models so as to systematically refine and validate the model. We use a structured format for encoding the conditions of many standard physiological and pharmacological experiments, specifying which parts of the model are involved, and comparing experiment outcomes with model output. A database of such experiments is run against successive generations of composite cellular models to iteratively improve the model against each experiment, while retaining global model validity. We suggest that this toolchain provides a principled and scalable way to tackle model complexity and diversity of data sources.

Challenges in Reproducibility, Replicability, and Comparability of Computational Models and Tools for Neuronal and Glial Networks, Cells, and Subcellular Structures

The possibility to replicate and reproduce published research results is one of the biggest challenges in all areas of science. In computational neuroscience, there are thousands of models available. However, it is rarely possible to reimplement the models based on the information in the original publication, let alone rerun the models just because the model implementations have not been made publicly available. We evaluate and discuss the comparability of a versatile choice of simulation tools: tools for biochemical reactions and spiking neuronal networks, and relatively new tools for growth in cell cultures. The replicability and reproducibility issues are considered for computational models that are equally diverse, including the models for intracellular signal transduction of neurons and glial cells, in addition to single glial cells, neuron-glia interactions, and selected examples of spiking neuronal networks. We also address the comparability of the simulation results with one another to comprehend if the studied models can be used to answer similar research questions. In addition to presenting the challenges in reproducibility and replicability of published results in computational neuroscience, we highlight the need for developing recommendations and good practices for publishing simulation tools and computational models. Model validation and flexible model description must be an integral part of the tool used to simulate and develop computational models. Constant improvement on experimental techniques and recording protocols leads to increasing knowledge about the biophysical mechanisms in neural systems. This poses new challenges for computational neuroscience: extended or completely new computational methods and models may be required. Careful evaluation and categorization of the existing models and tools provide a foundation for these future needs, for constructing multiscale models or extending the models to incorporate additional or more detailed biophysical mechanisms. Improving the quality of publications in computational neuroscience, enabling progressive building of advanced computational models and tools, can be achieved only through adopting publishing standards which underline replicability and reproducibility of research results.

Computational Models for Calcium-Mediated Astrocyte Functions

The computational neuroscience field has heavily concentrated on the modeling of neuronal functions, largely ignoring other brain cells, including one type of glial cell, the astrocytes. Despite the short history of modeling astrocytic functions, we were delighted about the hundreds of models developed so far to study the role of astrocytes, most often in calcium dynamics, synchronization, information transfer, and plasticity in vitro, but also in vascular events, hyperexcitability, and homeostasis. Our goal here is to present the state-of-the-art in computational modeling of astrocytes in order to facilitate better understanding of the functions and dynamics of astrocytes in the brain. Due to the large number of models, we concentrated on a hundred models that include biophysical descriptions for calcium signaling and dynamics in astrocytes. We categorized the models into four groups: single astrocyte models, astrocyte network models, neuron-astrocyte synapse models, and neuron-astrocyte network models to ease their use in future modeling projects. We characterized the models based on which earlier models were used for building the models and which type of biological entities were described in the astrocyte models. Features of the models were compared and contrasted so that similarities and differences were more readily apparent. We discovered that most of the models were basically generated from a small set of previously published models with small variations. However, neither citations to all the previous models with similar core structure nor explanations of what was built on top of the previous models were provided, which made it possible, in some cases, to have the same models published several times without an explicit intention to make new predictions about the roles of astrocytes in brain functions. Furthermore, only a few of the models are available online which makes it difficult to reproduce the simulation results and further develop the models. Thus, we would like to emphasize that only via reproducible research are we able to build better computational models for astrocytes, which truly advance science. Our study is the first to characterize in detail the biophysical and biochemical mechanisms that have been modeled for astrocytes.

Multi-scale modeling of altered synaptic plasticity related to Amyloid β effects

As suggested by Palop and Mucke (2010) pathologically elevated β-amyloid (Aβ) impairs long term potentiation (LTP) and enhances long term depression (LTD) possible underlying mechanisms in Alzheimer's Disease (AD). In the present paper we adopt and further elaborate a phenomenological computational model of bidirectional plasticity based on the calcium control hypothesis of Shouval et al. (2002). First, to account for Aβ effects the activation function Ω was modified assuming competition between LTP and LTD, and parameter sets were identified that well describe both normal and pathological synaptic plasticity processes. Second, a biophysically plausible kinetic model of bidirectional synaptic plasticity by D'Alcantara et al. (2003) was used to support findings of the phenomenological model and to further explain underlying kinetic processes. Model fitting pointed out molecular contributors, particularly calcineurin and type 1 protein phosphatase that might contribute to observed physiological disturbances in AD.

Impact of synaptic localization and subunit composition of ionotropic glutamate receptors on synaptic function: modeling and simulation studies

Ionotropic NMDA and AMPA glutamate receptors (iGluRs) play important roles in synaptic function under physiological and pathological conditions. iGluRs sub-synaptic localization and subunit composition are dynamically regulated by activity-dependent insertion and internalization. However, understanding the impact on synaptic transmission of changes in composition and localization of iGluRs is difficult to address experimentally. To address this question, we developed a detailed computational model of glutamatergic synapses, including spine and dendritic compartments, elementary models of subtypes of NMDA and AMPA receptors, glial glutamate transporters, intracellular calcium and a calcium-dependent signaling cascade underlying the development of long-term potentiation (LTP). These synapses were distributed on a neuron model and numerical simulations were performed to assess the impact of changes in composition and localization (synaptic vs extrasynaptic) of iGluRs on synaptic transmission and plasticity following various patterns of presynaptic stimulation. In addition, the effects of various pharmacological compounds targeting NMDARs or AMPARs were determined. Our results showed that changes in NMDAR localization have a greater impact on synaptic plasticity than changes in AMPARs. Moreover, the results suggest that modulators of AMPA and NMDA receptors have differential effects on restoring synaptic plasticity under different experimental situations mimicking various human diseases.

Reproducibility and Comparability of Computational Models for Astrocyte Calcium Excitability

The scientific community across all disciplines faces the same challenges of ensuring accessibility, reproducibility, and efficient comparability of scientific results. Computational neuroscience is a rapidly developing field, where reproducibility and comparability of research results have gained increasing interest over the past years. As the number of computational models of brain functions is increasing, we chose to address reproducibility using four previously published computational models of astrocyte excitability as an example. Although not conventionally taken into account when modeling neuronal systems, astrocytes have been shown to take part in a variety of in vitro and in vivo phenomena including synaptic transmission. Two of the selected astrocyte models describe spontaneous calcium excitability, and the other two neurotransmitter-evoked calcium excitability. We specifically addressed how well the original simulation results can be reproduced with a reimplementation of the models. Additionally, we studied how well the selected models can be reused and whether they are comparable in other stimulation conditions and research settings. Unexpectedly, we found out that three of the model publications did not give all the necessary information required to reimplement the models. In addition, we were able to reproduce the original results of only one of the models completely based on the information given in the original publications and in the errata. We actually found errors in the equations provided by two of the model publications; after modifying the equations accordingly, the original results were reproduced more accurately. Even though the selected models were developed to describe the same biological event, namely astrocyte calcium excitability, the models behaved quite differently compared to one another. Our findings on a specific set of published astrocyte models stress the importance of proper validation of the models against experimental wet-lab data from astrocytes as well as the careful review process of models. A variety of aspects of model development could be improved, including the presentation of models in publications and databases. Specifically, all necessary mathematical equations, as well as parameter values, initial values of variables, and stimuli used should be given precisely for successful reproduction of scientific results.

A Cognitive Model Based on Neuromodulated Plasticity

Associative learning, including classical conditioning and operant conditioning, is regarded as the most fundamental type of learning for animals and human beings. Many models have been proposed surrounding classical conditioning or operant conditioning. However, a unified and integrated model to explain the two types of conditioning is much less studied. Here, a model based on neuromodulated synaptic plasticity is presented. The model is bioinspired including multistored memory module and simulated VTA dopaminergic neurons to produce reward signal. The synaptic weights are modified according to the reward signal, which simulates the change of associative strengths in associative learning. The experiment results in real robots prove the suitability and validity of the proposed model.

Understanding the role of synaptopodin and the spine apparatus in Hebbian synaptic plasticity - New perspectives and the need for computational modeling

Synaptopodin (SP) is a proline-rich actin-associated protein essential for the formation of a spine apparatus (SA) in dendritic spines. The SA consists of stacks of smooth endoplasmic reticulum (sER) contiguous with the meshwork of somatodendritic ER. Spines of SP-deficient mice contain sER but no SA, demonstrating that SP is necessary for the assembly of ER cisterns into the more complex SA organelle. Although the SA was described decades ago, its function was difficult to investigate and remained elusive, in part because reliable markers for the SA were missing. After SP was identified as an essential component and a reliable marker of the SA, a role of SP/SA in hippocampal synaptic plasticity could be firmly established using loss-of-function approaches. Further studies revealed that SP/SA participate in the regulation of Ca²⁺-dependent spine-specific Hebbian plasticity and in activity-dependent changes in the spine actin cytoskeleton. In this review we are summarizing recent progress made on SP/SA in Hebbian plasticity and discuss open questions such as causality, spatiotemporal dynamics and complementarity of SP/SA-dependent mechanisms. We are proposing that computational modeling of spine Ca²⁺-signaling and actin remodeling pathways could address some of these issues and could indicate future research directions. Moreover, reaction-diffusion simulations could help to identify key feedforward and feedback regulatory motifs regulating the switch between an LTP and an LTD signaling module in SP/SA-containing spines, thus helping to find a unified view of SP/SA action in Hebbian plasticity.

Modelling bidirectional modulations in synaptic plasticity: A biochemical pathway model to understand the emergence of long term potentiation (LTP) and long term depression (LTD)

Synaptic plasticity induces bidirectional modulations of the postsynaptic response following a synaptic transmission. The long term forms of synaptic plasticity, named long term potentiation (LTP) and long term depression (LTD), are critical for the antithetic functions of the memory system, memory formation and removal, respectively. A common Ca(2+) signalling upstream triggers both LTP and LTD, and the critical proteins and factors coordinating the LTP/LTD inductions are not well understood. We develop an integrated model based on the sub-models of the indispensable synaptic proteins in the emergence of synaptic plasticity to validate and understand their potential roles in the expression of synaptic plasticity. The model explains Ca(2+)/calmodulin (CaM) complex dependent coordination of LTP/LTD expressions by the interactions among the indispensable proteins using the experimentally estimated kinetic parameters. Analysis of the integrated model provides us with insights into the effective timescales of the key proteins and we conclude that the CaM pool size is critical for the coordination between LTP/LTD expressions.

Pancreatic Beta Cell G-Protein Coupled Receptors and Second Messenger Interactions: A Systems Biology Computational Analysis

Insulin secretory in pancreatic beta-cells responses to nutrient stimuli and hormonal modulators include multiple messengers and signaling pathways with complex interdependencies. Here we present a computational model that incorporates recent data on glucose metabolism, plasma membrane potential, G-protein-coupled-receptors (GPCR), cytoplasmic and endoplasmic reticulum calcium dynamics, cAMP and phospholipase C pathways that regulate interactions between second messengers in pancreatic beta-cells. The values of key model parameters were inferred from published experimental data. The model gives a reasonable fit to important aspects of experimentally measured metabolic and second messenger concentrations and provides a framework for analyzing the role of metabolic, hormones and neurotransmitters changes on insulin secretion. Our analysis of the dynamic data provides support for the hypothesis that activation of Ca²⁺-dependent adenylyl cyclases play a critical role in modulating the effects of glucagon-like peptide 1 (GLP-1), glucose-dependent insulinotropic polypeptide (GIP) and catecholamines. The regulatory properties of adenylyl cyclase isoforms determine fluctuations in cytoplasmic cAMP concentration and reveal a synergistic action of glucose, GLP-1 and GIP on insulin secretion. On the other hand, the regulatory properties of phospholipase C isoforms determine the interaction of glucose, acetylcholine and free fatty acids (FFA) (that act through the FFA receptors) on insulin secretion. We found that a combination of GPCR agonists activating different messenger pathways can stimulate insulin secretion more effectively than a combination of GPCR agonists for a single pathway. This analysis also suggests that the activators of GLP-1, GIP and FFA receptors may have a relatively low risk of hypoglycemia in fasting conditions whereas an activator of muscarinic receptors can increase this risk. This computational analysis demonstrates that study of second messenger pathway interactions will improve understanding of critical regulatory sites, how different GPCRs interact and pharmacological targets for modulating insulin secretion in type 2 diabetes.

A Computational Model for the AMPA Receptor Phosphorylation Master Switch Regulating Cerebellar Long-Term Depression

The expression of long-term depression (LTD) in cerebellar Purkinje cells results from the internalisation of α-amino-3-hydroxy-5-methylisoxazole-4-propionic acid receptors (AMPARs) from the postsynaptic membrane. This process is regulated by a complex signalling pathway involving sustained protein kinase C (PKC) activation, inhibition of serine/threonine phosphatase, and an active protein tyrosine phosphatase, PTPMEG. In addition, two AMPAR-interacting proteins-glutamate receptor-interacting protein (GRIP) and protein interacting with C kinase 1 (PICK1)-regulate the availability of AMPARs for trafficking between the postsynaptic membrane and the endosome. Here we present a new computational model of these overlapping signalling pathways. The model reveals how PTPMEG cooperates with PKC to drive LTD expression by facilitating the effect of PKC on the dissociation of AMPARs from GRIP and thus their availability for trafficking. Model simulations show that LTD expression is increased by serine/threonine phosphatase inhibition, and negatively regulated by Src-family tyrosine kinase activity, which restricts the dissociation of AMPARs from GRIP under basal conditions. We use the model to expose the dynamic balance between AMPAR internalisation and reinsertion, and the phosphorylation switch responsible for the perturbation of this balance and for the rapid plasticity initiation and regulation. Our model advances the understanding of PF-PC LTD regulation and induction, and provides a validated extensible platform for more detailed studies of this fundamental synaptic process.

Systems biology of synaptic plasticity: a review on N-methyl-D-aspartate receptor mediated biochemical pathways and related mathematical models

Synaptic plasticity, an emergent property of synaptic networks, has shown strong correlation to one of the essential functions of the brain, memory formation. Through understanding synaptic plasticity, we hope to discover the modulators and mechanisms that trigger memory formation. In this paper, we first review the well understood modulators and mechanisms underlying N-methyl-D-aspartate receptor dependent synaptic plasticity, a major form of synaptic plasticity in hippocampus, and then comment on the key mathematical modelling approaches available in the literature to understand synaptic plasticity as the integration of the established functionalities of synaptic components.

The delicate bistability of CaMKII

Calcium/calmodulin-dependent protein kinase II (CaMKII) is a synaptic, autophosphorylating kinase that is essential for learning and memory. Previous models have suggested that CaMKII functions as a bistable switch that could be the molecular correlate of long-term memory, but experiments have failed to validate these predictions. These models involved significant approximations to overcome the combinatorial complexity inherent in a multisubunit, multistate system. Here, we develop a stochastic particle-based model of CaMKII activation and dynamics that overcomes combinatorial complexity without significant approximations. We report four major findings. First, the CaMKII model system is never bistable at resting calcium concentrations, which suggests that CaMKII activity does not function as the biochemical switch underlying long-term memory. Second, the steady-state activation curves are either laserlike or steplike. Both are characterized by a well-defined threshold for activation, which suggests that thresholding is a robust feature of this system. Third, transiently activated CaMKII can maintain its activity over the time course of many experiments, and such slow deactivation may account for the few reports of bistability in the literature. And fourth, under in vivo conditions, increases in phosphatase activity can increase CaMKII activity. This is a surprising and counterintuitive effect, as dephosphorylation is generally associated with CaMKII deactivation.

Shaping synaptic learning by the duration of postsynaptic action potential in a new STDP model

Single spikes and their timing matter in changing synaptic efficacy, which is known as spike-timing-dependent plasticity (STDP). Most previous studies treated spikes as all-or-none events, and considered their duration and magnitude as negligible. Here we explore the effects of action potential (AP) duration on synaptic plasticity in a simplified model neuron using computer simulations. We propose a novel STDP model that depresses synapses using an AP duration dependent LTD window and induces potentiation of synaptic strength when presynaptic spikes arrive before and during a postsynaptic AP (dSTDP). We demonstrate that AP duration is another key factor for insensitizing the postsynaptic neural firing and for controlling the shape of synaptic weight distribution. Extended AP durations produce a wide unimodal weight distribution that resembles the ones reported experimentally and make the postsynaptic neuron tranquil when disturbed by poisson noise spike trains, while equivalently sensitive to the synchronized. Our results suggest that the impact of AP duration, modeled here as an AP-dependent STDP window, on synaptic plasticity can be dramatic and should motivate future STDP studies.

Astrocyte-neuron interactions: from experimental research-based models to translational medicine

In this chapter, we review the principal astrocyte functions and the interactions between neurons and astrocytes. We then address how the experimentally observed functions have been verified in computational models and review recent experimental literature on astrocyte-neuron interactions. Benefits of computational neuroscience work are highlighted through selected studies with neurons and astrocytes by analyzing the existing models qualitatively and assessing the relevance of these models to experimental data. Common strategies to mathematical modeling and computer simulation in neuroscience are summarized for the nontechnical reader. The astrocyte-neuron interactions are then further illustrated by examples of some neurological and neurodegenerative diseases, where the miscommunication between glia and neurons is found to be increasingly important.

Heterosynaptic plasticity prevents runaway synaptic dynamics

Spike timing-dependent plasticity (STDP) and other conventional Hebbian-type plasticity rules are prone to produce runaway dynamics of synaptic weights. Once potentiated, a synapse would have higher probability to lead to spikes and thus to be further potentiated, but once depressed, a synapse would tend to be further depressed. The runaway synaptic dynamics can be prevented by precisely balancing STDP rules for potentiation and depression; however, experimental evidence shows a great variety of potentiation and depression windows and magnitudes. Here we show that modifications of synapses to layer 2/3 pyramidal neurons from rat visual and auditory cortices in slices can be induced by intracellular tetanization: bursts of postsynaptic spikes without presynaptic stimulation. Induction of these heterosynaptic changes depended on the rise of intracellular calcium, and their direction and magnitude correlated with initial state of release mechanisms. We suggest that this type of plasticity serves as a mechanism that stabilizes the distribution of synaptic weights and prevents their runaway dynamics. To test this hypothesis, we develop a cortical neuron model implementing both homosynaptic (STDP) and heterosynaptic plasticity with properties matching the experimental data. We find that heterosynaptic plasticity effectively prevented runaway dynamics for the tested range of STDP and input parameters. Synaptic weights, although shifted from the original, remained normally distributed and nonsaturated. Our study presents a biophysically constrained model of how the interaction of different forms of plasticity--Hebbian and heterosynaptic--may prevent runaway synaptic dynamics and keep synaptic weights unsaturated and thus capable of further plastic changes and formation of new memories.

System Identification of mGluR-Dependent Long-Term Depression

Recent advances have started to uncover the underlying mechanisms of metabotropic glutamate receptor (mGluR)–dependent long-term depression (LTD). However, it is not completely clear how these mechanisms are linked, and it is believed that several crucial mechanisms remain to be revealed. In this study, we investigated whether system identification (SI) methods can be used to gain insight into the mechanisms of synaptic plasticity. SI methods have been shown to be an objective and powerful approach for describing how sensory neurons encode information about stimuli. However, to our knowledge, it is the first time that SI methods have been applied to electrophysiological brain slice recordings of synaptic plasticity responses. The results indicate that the SI approach is a valuable tool for reverse-engineering of mGluR-LTD responses. We suggest that such SI methods can aid in unraveling the complexities of synaptic function.

Transfer functions for protein signal transduction: application to a model of striatal neural plasticity

We present a novel formulation for biochemical reaction networks in the context of protein signal transduction. The model consists of input-output transfer functions, which are derived from differential equations, using stable equilibria. We select a set of “source” species, which are interpreted as input signals. Signals are transmitted to all other species in the system (the “target” species) with a specific delay and with a specific transmission strength. The delay is computed as the maximal reaction time until a stable equilibrium for the target species is reached, in the context of all other reactions in the system. The transmission strength is the concentration change of the target species. The computed input-output transfer functions can be stored in a matrix, fitted with parameters, and even recalled to build dynamical models on the basis of state changes. By separating the temporal and the magnitudinal domain we can greatly simplify the computational model, circumventing typical problems of complex dynamical systems. The transfer function transformation of biochemical reaction systems can be applied to mass-action kinetic models of signal transduction. The paper shows that this approach yields significant novel insights while remaining a fully testable and executable dynamical model for signal transduction. In particular we can deconstruct the complex system into local transfer functions between individual species. As an example, we examine modularity and signal integration using a published model of striatal neural plasticity. The modularizations that emerge correspond to a known biological distinction between calcium-dependent and cAMP-dependent pathways. Remarkably, we found that overall interconnectedness depends on the magnitude of inputs, with higher connectivity at low input concentrations and significant modularization at moderate to high input concentrations. This general result, which directly follows from the properties of individual transfer functions, contradicts notions of ubiquitous complexity by showing input-dependent signal transmission inactivation.

Calcium and cAMP levels interact to determine attraction versus repulsion in axon guidance

Correct guidance of axons to their targets depends on an intricate network of signaling molecules in the growth cone. Calcium and cAMP are two key regulators of whether axons are attracted or repelled by molecular gradients, but how these molecules interact to determine guidance responses remains unclear. Here, we constructed a mathematical model for the relevant signaling network, which explained a large range of previous biological data and made predictions for when axons will be attracted or repelled. We then confirmed these predictions experimentally, in particular showing that while small increases in cAMP levels promote attraction large increases do not, and that under some circumstances reducing cAMP levels promotes attraction. Together, these results show that a relatively simple mathematical model can quantitatively predict guidance decisions across a wide range of conditions, and that calcium and cAMP levels play a more complex role in these decisions than previously determined.

CaMKII activation and dynamics are independent of the holoenzyme structure: an infinite subunit holoenzyme approximation

The combinatorial explosion produced by the multi-state, multi-subunit character of CaMKII has made analysis and modeling of this key signaling protein a significant challenge. Using rule-based and particle-based approaches, we construct exact models of CaMKII holoenzyme dynamics and study these models as a function of the number of subunits per holoenzyme, N. Without phosphatases the dynamics of activation are independent of the holoenzyme structure unless phosphorylation significantly alters the kinase activity of a subunit. With phosphatases the model is independent of holoenzyme size for N > 6. We introduce an infinite subunit holoenzyme approximation (ISHA), which simplifies the modeling by eliminating the combinatorial complexities encountered in any finite holoenzyme model. The ISHA is an excellent approximation to the full system over a broad range of physiologically relevant parameters. Finally, we demonstrate that the ISHA reproduces the behavior of exact models during synaptic plasticity protocols, which justifies its use as a module in large models of synaptic plasticity.

Modeling signal transduction leading to synaptic plasticity: evaluation and comparison of five models

An essential phenomenon of the functional brain is synaptic plasticity which is associated with changes in the strength of synapses between neurons. These changes are affected by both extracellular and intracellular mechanisms. For example, intracellular phosphorylation-dephosphorylation cycles have been shown to possess a special role in synaptic plasticity. We, here, provide the first computational comparison of models for synaptic plasticity by evaluating five models describing postsynaptic signal transduction networks. Our simulation results show that some of the models change their behavior completely due to varying total concentrations of protein kinase and phosphatase. Furthermore, the responses of the models vary when models are compared to each other. Based on our study, we conclude that there is a need for a general setup to objectively compare the models and an urgent demand for the minimum criteria that a computational model for synaptic plasticity needs to meet.