Calcium-dependent calcium decay explains STDP in a dynamic model of hippocampal synapses

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.

GluN2B-NMDAR subunit contribution on synaptic plasticity: A phenomenological model for CA3-CA1 synapses

Synaptic plasticity is believed to be a key mechanism underlying learning and memory. We developed a phenomenological N-methyl-D-aspartate (NMDA) receptor-based voltage-dependent synaptic plasticity model for synaptic modifications at hippocampal CA3-CA1 synapses on a hippocampal CA1 pyramidal neuron. The model incorporates the GluN2A-NMDA and GluN2B-NMDA receptor subunit-based functions and accounts for the synaptic strength dependence on the postsynaptic NMDA receptor composition and functioning without explicitly modeling the NMDA receptor-mediated intracellular calcium, a local trigger of synaptic plasticity. We embedded the model into a two-compartmental model of a hippocampal CA1 pyramidal cell and validated it against experimental data of spike-timing-dependent synaptic plasticity (STDP), high and low-frequency stimulation. The developed model predicts altered learning rules in synapses formed on the apical dendrites of the detailed compartmental model of CA1 pyramidal neuron in the presence of the GluN2B-NMDA receptor hypofunction and can be used in hippocampal networks to model learning in health and disease.

Synaptic Plasticity Is Predicted by Spatiotemporal Firing Rate Patterns and Robust to In Vivo-like Variability

Synaptic plasticity, the experience-induced change in connections between neurons, underlies learning and memory in the brain. Most of our understanding of synaptic plasticity derives from in vitro experiments with precisely repeated stimulus patterns; however, neurons exhibit significant variability in vivo during repeated experiences. Further, the spatial pattern of synaptic inputs to the dendritic tree influences synaptic plasticity, yet is not considered in most synaptic plasticity rules. Here, we investigate how spatiotemporal synaptic input patterns produce plasticity with in vivo-like conditions using a data-driven computational model with a plasticity rule based on calcium dynamics. Using in vivo spike train recordings as inputs to different size clusters of spines, we show that plasticity is strongly robust to trial-to-trial variability of spike timing. In addition, we derive general synaptic plasticity rules describing how spatiotemporal patterns of synaptic inputs control the magnitude and direction of plasticity. Synapses that strongly potentiated have greater firing rates and calcium concentration later in the trial, whereas strongly depressing synapses have hiring firing rates early in the trial. The neighboring synaptic activity influences the direction and magnitude of synaptic plasticity, with small clusters of spines producing the greatest increase in synaptic strength. Together, our results reveal that calcium dynamics can unify diverse plasticity rules and reveal how spatiotemporal firing rate patterns control synaptic plasticity.

Long-Term Effects of Platelet-Rich Fibrin on Fat Graft Survival and Their Optimal Mixing Ratio

BACKGROUND

Platelet-rich fibrin (PRF) can promote fat graft survival, but the reported mixing ratio of PRF to fat ranges from 1:25 to 1:2, lacking a clear standard for clinical application.

OBJECTIVES

The authors sought to explore the long-term effects of PRF on grafted fat and their optimal mixing ratio.

METHODS

Nude mice were randomly divided into a control group (receiving subcutaneous injection of fat granules) and 4 PRF groups (receiving subcutaneous injection of PRF and fat granules at volume ratios of 1:5, 1:10, 1:15, and 1:20, respectively). The graft samples (n = 12) were obtained in weeks 4, 8, and 12 to (1) calculate retention rates; (2) evaluate gene and protein expression of vascular endothelial growth factor A (VEGF-A), peroxisome proliferator-activated receptor-γ (PPAR-γ), type I collagen A1 (COL1-A1), and B-cell lymphoma-2 associated X protein (BAX); (3) perform hematoxylin and eosin, Masson's trichrome, α-smooth muscle action, and periplipin-1 stainings; and (4) count the microvessels and viable adipocytes.

RESULTS

Compared with the control group, PRF groups had higher retention rates, a higher gene/protein expression of VEGF-A, a lower gene/protein expression of COL1-A1 and BAX, less fibrosis, and more microvessels and viable adipocytes. Group 1:10 was superior to other groups in terms of retention rates and other evaluation indexes. The expression of PPAR-γ did not significantly differ among groups.

CONCLUSIONS

PRF may not play a long-term effect on adipogenesis, but it can still promote fat graft survival through facilitating vascularization, regulating collagen production, and inhibiting apoptosis. PRF can achieve the best promoting effect when the mixing ratio of PRF to fat is 1:10, which is recommended as the optimal ratio for clinical application.

Dendritic spine geometry and spine apparatus organization govern the spatiotemporal dynamics of calcium

Dendritic spines are small subcompartments that protrude from the dendrites of neurons and are important for signaling activity and synaptic communication. These subcompartments have been characterized to have different shapes. While it is known that these shapes are associated with spine function, the specific nature of these shape-function relationships is not well understood. In this work, we systematically investigated the relationship between the shape and size of both the spine head and spine apparatus, a specialized endoplasmic reticulum compartment within the spine head, in modulating rapid calcium dynamics using mathematical modeling. We developed a spatial multicompartment reaction-diffusion model of calcium dynamics in three dimensions with various flux sources, including N-methyl-D-aspartate receptors (NMDARs), voltage-sensitive calcium channels (VSCCs), and different ion pumps on the plasma membrane. Using this model, we make several important predictions. First, the volume to surface area ratio of the spine regulates calcium dynamics. Second, membrane fluxes impact calcium dynamics temporally and spatially in a nonlinear fashion. Finally, the spine apparatus can act as a physical buffer for calcium by acting as a sink and rescaling the calcium concentration. These predictions set the stage for future experimental investigations of calcium dynamics in dendritic spines.

Field-Driven Athermal Activation of Amorphous Metal Oxide Semiconductors for Flexible Programmable Logic Circuits and Neuromorphic Electronics

Despite extensive research, large-scale realization of metal-oxide electronics is still impeded by high-temperature fabrication, incompatible with flexible substrates. Ideally, an athermal treatment modifying the electronic structure of amorphous metal oxide semiconductors (AMOS) to generate sufficient carrier concentration would help mitigate such high-temperature requirements, enabling realization of high-performance electronics on flexible substrates. Here, a novel field-driven athermal activation of AMOS channels is demonstrated via an electrolyte-gating approach. Facilitating migration of charged oxygen species across the semiconductor-dielectric interface, this approach modulates the local electronic structure of the channel, generating sufficient carriers for charge transport and activating oxygen-compensated thin films. The thin-film transistors (TFTs) investigated here depict an enhancement of linear mobility from 51 to 105.25 cm² V^-1 s^-1 (ionic-gated) and from 8.09 to 14.49 cm² V^-1 s^-1 (back-gated), by creating additional oxygen vacancies. The accompanying stochiometric transformations, monitored via spectroscopic measurements (X-ray photoelectron spectroscopy) corroborate the detailed electrical (TFT, current evolution) parameter analyses, providing critical insights into the underlying oxygen-vacancy generation mechanism and clearly demonstrating field-induced activation as a promising alternative to conventional high-temperature annealing strategies. Facilitating on-demand active programing of the operation modes of transistors (enhancement vs depletion), this technique paves way for facile fabrication of logic circuits and neuromorphic transistors for bioinspired computing.

A Glutamatergic Spine Model to Enable Multi-Scale Modeling of Nonlinear Calcium Dynamics

In synapses, calcium is required for modulating synaptic transmission, plasticity, synaptogenesis, and synaptic pruning. The regulation of calcium dynamics within neurons involves cellular mechanisms such as synaptically activated channels and pumps, calcium buffers, and calcium sequestrating organelles. Many experimental studies tend to focus on only one or a small number of these mechanisms, as technical limitations make it difficult to observe all features at once. Computational modeling enables incorporation of many of these properties together, allowing for more complete and integrated studies. However, the scale of existing detailed models is often limited to synaptic and dendritic compartments as the computational burden rapidly increases when these models are integrated in cellular or network level simulations. In this article we present a computational model of calcium dynamics at the postsynaptic spine of a CA1 pyramidal neuron, as well as a methodology that enables its implementation in multi-scale, large-scale simulations. We first present a mechanistic model that includes individually validated models of various components involved in the regulation of calcium at the spine. We validated our mechanistic model by comparing simulated calcium levels to experimental data found in the literature. We performed additional simulations with the mechanistic model to determine how the simulated calcium activity varies with respect to presynaptic-postsynaptic stimulation intervals and spine distance from the soma. We then developed an input-output (IO) model that complements the mechanistic calcium model and provide a computationally efficient representation for use in larger scale modeling studies; we show the performance of the IO model compared to the mechanistic model in terms of accuracy and speed. The models presented here help achieve two objectives. First, the mechanistic model provides a comprehensive platform to describe spine calcium dynamics based on individual contributing factors. Second, the IO model is trained on the main dynamical features of the mechanistic model and enables nonlinear spine calcium modeling on the cell and network level simulation scales. Utilizing both model representations provide a multi-level perspective on calcium dynamics, originating from the molecular interactions at spines and propagating the effects to higher levels of activity involved in network behavior.

Modulation of Spike-Timing Dependent Plasticity: Towards the Inclusion of a Third Factor in Computational Models

In spike-timing dependent plasticity (STDP) change in synaptic strength depends on the timing of pre- vs. postsynaptic spiking activity. Since STDP is in compliance with Hebb's postulate, it is considered one of the major mechanisms of memory storage and recall. STDP comprises a system of two coincidence detectors with N-methyl-D-aspartate receptor (NMDAR) activation often posited as one of the main components. Numerous studies have unveiled a third component of this coincidence detection system, namely neuromodulation and glia activity shaping STDP. Even though dopaminergic control of STDP has most often been reported, acetylcholine, noradrenaline, nitric oxide (NO), brain-derived neurotrophic factor (BDNF) or gamma-aminobutyric acid (GABA) also has been shown to effectively modulate STDP. Furthermore, it has been demonstrated that astrocytes, via the release or uptake of glutamate, gate STDP expression. At the most fundamental level, the timing properties of STDP are expected to depend on the spatiotemporal dynamics of the underlying signaling pathways. However in most cases, due to technical limitations experiments grant only indirect access to these pathways. Computational models carefully constrained by experiments, allow for a better qualitative understanding of the molecular basis of STDP and its regulation by neuromodulators. Recently, computational models of calcium dynamics and signaling pathway molecules have started to explore STDP emergence in ex and in vivo-like conditions. These models are expected to reproduce better at least part of the complex modulation of STDP as an emergent property of the underlying molecular pathways. Elucidation of the mechanisms underlying STDP modulation and its consequences on network dynamics is of critical importance and will allow better understanding of the major mechanisms of memory storage and recall both in health and disease.

Pavlovian conditioning demonstrated with neuromorphic memristive devices

Pavlovian conditioning, a classical case of associative learning in a biological brain, is demonstrated using the Ni/Nb-SrTiO₃/Ti memristive device with intrinsic forgetting properties in the framework of the asymmetric spike-timing-dependent plasticity of synapses. Three basic features of the Pavlovian conditioning, namely, acquisition, extinction and recovery, are implemented in detail. The effects of the temporal relation between conditioned and unconditioned stimuli as well as the time interval between individual training trials on the Pavlovian conditioning are investigated. The resulting change of the response strength, the number of training trials necessary for acquisition and the number of extinction trials are illustrated. This work clearly demonstrates the hardware implementation of the brain function of the associative learning.

Inferring Neuronal Dynamics from Calcium Imaging Data Using Biophysical Models and Bayesian Inference

Calcium imaging has been used as a promising technique to monitor the dynamic activity of neuronal populations. However, the calcium trace is temporally smeared which restricts the extraction of quantities of interest such as spike trains of individual neurons. To address this issue, spike reconstruction algorithms have been introduced. One limitation of such reconstructions is that the underlying models are not informed about the biophysics of spike and burst generations. Such existing prior knowledge might be useful for constraining the possible solutions of spikes. Here we describe, in a novel Bayesian approach, how principled knowledge about neuronal dynamics can be employed to infer biophysical variables and parameters from fluorescence traces. By using both synthetic and in vitro recorded fluorescence traces, we demonstrate that the new approach is able to reconstruct different repetitive spiking and/or bursting patterns with accurate single spike resolution. Furthermore, we show that the high inference precision of the new approach is preserved even if the fluorescence trace is rather noisy or if the fluorescence transients show slow rise kinetics lasting several hundred milliseconds, and inhomogeneous rise and decay times. In addition, we discuss the use of the new approach for inferring parameter changes, e.g. due to a pharmacological intervention, as well as for inferring complex characteristics of immature neuronal circuits.

Calcium imaging of single neurons enables the indirect observation of neuronal dynamics, for example action potential firing. In contrast to the precise timing of spike trains, the calcium trace is temporally rather smeared and measured as a fluorescence trace. Consequently, several methods have been proposed to reconstruct spikes from calcium imaging data. However, a common feature of these methods is that they are not based on the biophysics of how neurons fire spikes and bursts. We propose to introduce well-established biophysical models to create a direct link between neuronal dynamics, e.g. the membrane potential, and fluorescence traces. Using both synthetic and experimental data, we show that this approach not only provides a robust and accurate spike reconstruction but also a reliable inference about the biophysically relevant parameters and variables. This enables novel ways of analyzing calcium imaging experiments in terms of the underlying biophysical quantities.