Synaptic shot noise and conductance fluctuations affect the membrane voltage with equal significance

The subthresholdmembranevoltage of a neuron in active cortical tissue is a fluctuating quantity with a distribution that reflects the firing statistics of the presynaptic population. It was recently found that conductancebased synaptic drive can lead to distributions with a significant skew. Here it is demonstrated that the underlying shot noise caused by Poissonian spike arrival also skews the membrane distribution, but in the opposite sense. Using a perturbative method, we analyze the effects of shot noise on the distribution of synaptic conductances and calculate the consequent voltage distribution. To first order in the perturbation theory, the voltage distribution is a gaussian modulated by a prefactor that captures the skew. The gaussian component is identical to distributions derived using current-based models with an effective membrane time constant. The well-known effective-time-constant approximation can therefore be identified as the leading-order solution to the full conductance-based model. The higher-order modulatory prefactor containing the skew comprises terms due to both shot noise and conductance fluctuations. The diffusion approximation misses these shot-noise effects implying that analytical approaches such as the Fokker-Planck equation or simulation with filtered white noise cannot be used to improve on the gaussian approximation. It is further demonstrated that quantities used for fitting theory to experiment, such as the voltage mean and variance, are robust against these non-Gaussian effects. The effective-time-constant approximation is therefore relevant to experiment and provides a simple analytic base on which other pertinent biological details may be added

Spike shape and synaptic-amplitude distribution interact to set the high-frequency ring-rate response of neuronal populations

The dynamics of an ensemble of particles driven out of a potential well, with replacement, by the Poissonian arrival of amplitude-distributed shot noise is examined. A general formula for the high-frequency limit of the escape-rate susceptibility is derived. For certain choices of potential well and amplitude distribution the decay of the high-frequency susceptibility exhibits a nonuniversal exponent. This is a qualitatively different response to that predicted by the diffusion approximation. To provide an example the general framework is applied to a problem of current interest in the biophysics of neuronal voltage dynamics. It is shown that the firing-rate response of neurons to rapidly varying stimuli can be significantly enhanced depending on the ratio between the scale of excitatory postsynaptic potentials and the voltage range over which an action potential initiates. The result is robust to various choices of threshold definition and also to synaptic filtering at physiologically reasonable time scales

Adenosine A1 receptor activation mediates the developmental shift at layer 5 pyramidal cell synapses and is a determinant of mature synaptic strength

During the first postnatal month glutamatergic synapses between layer 5 pyramidal cells in the rodent neocortex switch from an immature state exhibiting high probability of neurotransmitter release, large unitary amplitude and synaptic depression to a mature state with decreased probability of release, smaller unitary amplitude and synaptic facilitation. Using paired recordings, we demonstrate that the developmental shift in release probability at synapses between rat somatosensory layer 5 thick-tufted pyramidal cells is due to a higher and more heterogeneous activation of presynaptic adenosine A1 receptors. Immature synapses under control conditions exhibited distributions of CV, failure rate and release probability that were almost coincident with the A1 receptor blocked condition; however, mature synapses under control conditions exhibited much broader distributions that spanned those of both the A1 receptor agonised and antagonised conditions. Immature and mature synapses expressed A1 receptors with no observable difference in functional efficacy and therefore the heterogeneous A1 receptor activation seen in the mature neocortex is due to increased adenosine concentrations that vary between synapses. Given the central role demonstrated for A1 receptor activation in determining synaptic amplitude and the statistics of transmission between mature layer 5 pyramidal cells, the emplacement of adenosine sources and sinks near the synaptic terminal could constitute a novel form of long-term synaptic plasticity

Transmission of temporally correlated spike trains through synapses with short-term depression

Short-term synaptic depression, caused by depletion of releasable neurotransmitter, modulates the strength of neuronal connections in a history-dependent manner. Quantifying the statistics of synaptic transmission requires stochastic models that link probabilistic neurotransmitter release with presynaptic spike-train statistics. Common approaches are to model the presynaptic spike train as either regular or a memory-less Poisson process: few analytical results are available that describe depressing synapses when the afferent spike train has more complex, temporally correlated statistics such as bursts. Here we present a series of analytical results—from vesicle release-site occupancy statistics, via neurotransmitter release, to the post-synaptic voltage mean and variance—for depressing synapses driven by correlated presynaptic spike trains. The class of presynaptic drive considered is that fully characterised by the inter-spike-interval distribution and encompasses a broad range of models used for neuronal circuit and network analyses, such as integrate-and-fire models with a complete post-spike reset and receiving sufficiently short-time correlated drive. We further demonstrate that the derived post-synaptic voltage mean and variance allow for a simple and accurate approximation of the firing rate of the post-synaptic neuron, using the exponential integrate-and-fire model as an example. These results extend the level of biological detail included in models of synaptic transmission and will allow for the incorporation of more complex and physiologically relevant firing patterns into future studies of neuronal networks

Short-Term Synaptic Plasticity Orchestrates the Response of Pyramidal Cells and Interneurons to Population Bursts

The synaptic drive from neuronal populations varies considerably over short time scales. Such changes in the pre-synaptic rate trigger many temporal processes absent under steady-state conditions. This paper examines the differential impact of pyramidal cell population bursts on post-synaptic pyramidal cells receiving depressing synapses, and on a class of interneuron that receives facilitating synapses. In experiment a significant shift of the order of one hundred milliseconds is seen between the response of these two cell classes to the same population burst. It is demonstrated here that such a temporal differentiation of the response can be explained by the synaptic and membrane properties without recourse to elaborate cortical wiring schemes. Experimental data is first used to construct models of the two types of dynamic synaptic response. A population-based approach is then followed to examine analytically the temporal synaptic filtering effects of the population burst for the two post-synaptic targets. The peak-to-peak delays seen in experiment can be captured by the model for experimentally realistic parameter ranges. It is further shown that the temporal separation of the response is communicated in the outgoing action potentials of the two post-synaptic cells: pyramidal cells fire at the beginning of the burst and the class of interneuron receiving facilitating synapses fires at the end of the burst. The functional role of such delays in the temporal organisation of activity in the cortical microcircuit is discusse

Bayesian Inference of Synaptic Quantal Parameters from Correlated Vesicle Release

Synaptic transmission is both history-dependent and stochastic, resulting in varying responses to presentations of the same presynaptic stimulus. This complicates attempts to infer synaptic parameters and has led to the proposal of a number of different strategies for their quantification. Recently Bayesian approaches have been applied to make more efficient use of the data collected in paired intracellular recordings. Methods have been developed that either provide a complete model of the distribution of amplitudes for isolated responses or approximate the amplitude distributions of a train of post-synaptic potentials, with correct short-term synaptic dynamics but neglecting correlations. In both cases the methods provided significantly improved inference of model parameters as compared to existing mean-variance fitting approaches. However, for synapses with high release probability, low vesicle number or relatively low restock rate and for data in which only one or few repeats of the same pattern are available, correlations between serial events can allow for the extraction of significantly more information from experiment: a more complete Bayesian approach would take this into account also. This has not been possible previously because of the technical difficulty in calculating the likelihood of amplitudes seen in correlated post-synaptic potential trains; however, recent theoretical advances have now rendered the likelihood calculation tractable for a broad class of synaptic dynamics models. Here we present a compact mathematical form for the likelihood in terms of a matrix product and demonstrate how marginals of the posterior provide information on covariance of parameter distributions. The associated computer code for Bayesian parameter inference for a variety of models of synaptic dynamics is provided in the supplementary material allowing for quantal and dynamical parameters to be readily inferred from experimental data sets

Linear and nonlinear integrate-and-fire neurons driven by synaptic shot noise with reversal potentials

The steady-state firing rate and firing-rate response of the leaky and exponential integrate-and-fire models receiving synaptic shot noise with excitatory and inhibitory reversal potentials is examined. For the particular case where the underlying synaptic conductances are exponentially distributed, it is shown that the master equation for a population of such model neurons can be reduced from an integrodifferential form to a more tractable set of three differential equations. The system is nevertheless more challenging analytically than for current-based synapses: where possible, analytical results are provided with an efficient numerical scheme and code provided for other quantities. The increased tractability of the framework developed supports an ongoing critical comparison between models in which synapses are treated with and without reversal potentials, such as recently in the context of networks with balanced excitatory and inhibitory conductances

Modelling microelectrode biosensors : free-flow calibration can substantially underestimate tissue concentrations

Microelectrode amperometric biosensors are widely used to measure concentrations of analytes in solution and tissue including acetylcholine, adenosine, glucose and glutamate. A great deal of experimental and modelling effort has been directed at quantifying the response of the biosensors themselves; however, the influence that the macroscopic tissue environment has on biosensor response has not been subjected to the same level of scrutiny. Here we identify an important issue in the way microelectrode biosensors are calibrated that is likely to have led to underestimations of analyte tissue concentrations. Concentration in tissue is typically determined by comparing the biosensor signal to that measured in free-flow calibration conditions. In a free-flow environment the concentration of the analyte at the outer surface of the biosensor can be considered constant. However, in tissue the analyte reaches the biosensor surface by diffusion through the extracellular space. Because the enzymes in the biosensor break down the analyte, a density gradient is set up resulting in a significantly lower concentration of analyte near the biosensor surface. This effect is compounded by the diminished volume fraction (porosity) and reduction in the diffusion coefficient due to obstructions (tortuosity) in tissue. We demonstrate this effect through modelling and experimentally verify our predictions in diffusive environments

Alpha-synuclein aggregates increase the conductance of substantia nigra dopamine neurons, an effect partly reversed by the KATP channel inhibitor glibenclamide

Dopaminergic neurons in the substantia nigra pars compacta (SNpc) form an important part of the basal ganglia circuitry, playing key roles in movement initiation and co-ordination. A hallmark of Parkinson’s disease (PD) is the degeneration of these SNpc dopaminergic neurons leading to akinesia, bradykinesia and tremor. There is gathering evidence that oligomeric alpha synuclein (α-syn) is one of the major pathological species in PD, with its deposition in Lewy bodies closely correlated with disease progression. However the precise mechanisms underlying the effects of oligomeric α-syn on dopaminergic neuron function have yet to be fully defined. Here we have combined electrophysiological recording and detailed analysis to characterise the time-dependent effects of α-syn aggregates (consisting of oligomers and possibly small fibrils) on the properties of SNpc dopaminergic neurons. The introduction of α-syn aggregates into single dopaminergic neurons via the patch electrode significantly reduced both the input resistance and the firing rate without changing the membrane potential. These effects occurred after 8-16 minutes of dialysis but did not occur with the monomeric form of α-syn. The effects of α-syn aggregates could be significantly reduced by pre-incubation with the ATP-sensitive potassium channel (KATP) inhibitor glibenclamide. This data suggests that accumulation of α-syn aggregates in dopaminergic neurons may chronically activate KATP channels leading to a significant loss of excitability and dopamine release

Experimentally verified parameter sets for modelling heterogeneous neocortical pyramidal-cell populations

Models of neocortical networks are increasingly including the diversity of excitatory and inhibitory neuronal classes. Significant variability in cellular properties are also seen within a nominal neuronal class and this heterogeneity can be expected to influence the population response and information processing in networks. Recent studies have examined the population and network effects of variability in a particular neuronal parameter with some plausibly chosen distribution. However, the empirical variability and covariance seen across multiple parameters are rarely included, partly due to the lack of data on parameter correlations in forms convenient for model construction. To addess this we quantify the heterogeneity within and between the neocortical pyramidal-cell classes in layers 2/3, 4, and the slender-tufted and thick-tufted pyramidal cells of layer 5 using a combination of intracellular recordings, single-neuron modelling and statistical analyses. From the response to both square-pulse and naturalistic fluctuating stimuli, we examined the class-dependent variance and covariance of electrophysiological parameters and identify the role of the h current in generating parameter correlations. A byproduct of the dynamic I-V method we employed is the straightforward extraction of reduced neuron models from experiment. Empirically these models took the refractory exponential integrate-and-fire form and provide an accurate fit to the perisomatic voltage responses of the diverse pyramidal-cell populations when the class-dependent statistics of the model parameters were respected. By quantifying the parameter statistics we obtained an algorithm which generates populations of model neurons, for each of the four pyramidal-cell classes, that adhere to experimentally observed marginal distributions and parameter correlations. As well as providing this tool, which we hope will be of use for exploring the effects of heterogeneity in neocortical networks, we also provide the code for the dynamic I-V method and make the full electrophysiological data set available