Variability of excitatory currents due to single-channel noise, receptor number and morphological heterogeneity

Patch clamp recordings of excitatory postsynaptic currents (EPSCs) in central neurons reveal large fluctuations in amplitudes and decay times of AMPA-receptor-mediated EPSCs. By using Monte Carlo simulations of synaptic transmission in brainstem interneurons, we tested several hypothesis that could account for the observed variability. The coefficient of variation (CV) of 0.5 for miniature amplitudes cannot be explained by fluctuations in vesicle content or receptor distribution, but is traced to variations in receptor number, which is estimated as 77+/-39 receptors per bouton. As the variability of rise times reflects fluctuations in size of the post-synaptic density and heterogeneity of the receptor distribution, the relatively small CV=0.37 of experimentally determined values points to a homogeneous arrangement of receptors. Within our model the large variability of decay times (CV=0.49) can only be explained by fluctuations in the transmitter time course (mean residence times of 0.4+/-0.13 ms), presumably resulting from heterogeneities in synaptic morphology. Hence, our simulations indicate that different noise sources control the variability of amplitudes, rise and decay times. In particular, the distribution of decay times yields information about the synaptic transmission process, which cannot be obtained from other observables.     

Modulation of glutamate mobility reveals the mechanism underlying slow-rising AMPAR EPSCs and the diffusion coefficient in the synaptic cleft

Fast- and slow-rising AMPA receptor-mediated EPSCs occur at central synapses. Fast-rising EPSCs are thought to be mediated by rapid local release of glutamate. However, two controversial mechanisms have been proposed to underlie slow-rising EPSCs: prolonged local release of transmitter via a fusion pore, and spillover of transmitter released rapidly from distant sites. We have investigated the mechanism underlying slow-rising EPSCs and the diffusion coefficient of glutamate in the synaptic cleft (Dglut) at cerebellar mossy fiber-granule cell synapses using a combination of diffusion modeling and patch-clamp recording. Simulations show that modulating Dglut has different effects on the peak amplitudes and time courses of EPSCs mediated by these two mechanisms. Slowing diffusion with the macromolecule dextran slowed slow-rising EPSCs and had little effect on their amplitude, indicating that glutamate spillover underlies these currents. Our results also suggest that under control conditions Dglut is approximately 3-fold lower than in free solution.     

A Brownian model of glutamate diffusion in excitatory synapses of hippocampus

We simulated the diffusion of glutamate, following the release of a single vesicle from a pre-synaptic terminal, in the synaptic cleft by using a Brownian diffusion model based on Langevin equations. The synaptic concentration time course and the time course of quantal excitatory post-synaptic current have been analyzed. The results showed that they depend on the number of receptors located at post-synaptic membrane. Their time course are dependent both on the total number of the post-synaptic receptors and on the eccentricity of the pre-synaptic glutamate vesicle.     

Spillover of glutamate onto synaptic AMPA receptors enhances fast transmission at a cerebellar synapse

Diffusion of glutamate from the synaptic cleft can activate high-affinity receptors, but is not thought to contribute to fast AMPA receptor-mediated transmission. Here, we show that single AMPA receptor EPSCs at the cerebellar mossy fiber-granule cell connection are mediated by both direct release of glutamate and rapid diffusion of glutamate from neighboring synapses. Immunogold localization revealed that AMPA receptors are located exclusively in postsynaptic densities, indicating that spillover of glutamate occurs between synaptic contacts. Spillover currents contributed half the synaptic charge and exhibited little trial-to-trial variability. We propose that spillover of glutamate improves transmission efficacy by both increasing the amplitude and duration of the EPSP and reducing fluctuations arising from the probabilistic nature of transmitter release.     

Glutamate escape from a tortuous synaptic cleft of the hippocampal mossy fibre synapse

The time course of neurotransmitter in the synaptic cleft contributes substantially to the fast kinetics of synaptic signalling. Hippocampal mossy fibres (MFs), a well-characterised excitatory pathway from dentate granule cells to the hippocampus proper, form large glutamatergic synapses at branched spiny structures in CA3 pyramidal cell dendrites. To what extent transmission at these synapses is affected by retarded glutamate clearance from the large tortuous synaptic cleft is not known. Here, we propose a simple geometrical approximation representing the 'typical' geometry of thorny excrescences that form the tortuous cleft interface at a MF synapse. We then employ Monte Carlo simulations to monitor movements of 3000 individual glutamate molecules released within the cleft. The results predict that, in the absence of neuronal glutamate transporters, it should take approximately 10 ms for 50% and 60-70 ms for 90% of glutamate molecules to escape the MF synapse.     

Control of Cleft Glutamate Concentration and Glutamate Spill-Out by Perisynaptic Glia: Uptake and Diffusion Barriers

Most glutamatergic synapses in the mammalian central nervous system are covered by thin astroglial processes that exert a dual action on synaptically released glutamate: they form physical barriers that oppose diffusion and they carry specific transporters that remove glutamate from the extracellular space. The present study was undertaken to investigate the dual action of glia by means of computer simulation. A realistic synapse model based on electron microscope data and Monte Carlo algorithms were used for this purpose. Results show (1) that physical obstacles formed by glial processes delay glutamate exit from the cleft and (2) that this effect is efficiently counteracted by glutamate uptake. Thus, depending on transporter densities, the presence of perisynaptic glia may result in increased or decreased glutamate transient in the synaptic cleft. Changes in temporal profiles of cleft glutamate concentration induced by glia differentially impact the response of the various synaptic and perisynaptic receptor subtypes. In particular, GluN2B- and GluN2C-NMDA receptor responses are strongly modified while GluN2A-NMDA receptor responses are almost unaffected. Thus, variations in glial transporter expression may allow differential tuning of NMDA receptors according to their subunit composition. In addition, simulation data suggest that the sink effect generated by transporters accumulation in the vicinity of the release site is the main mechanism limiting glutamate spill-out. Physical obstacles formed by glial processes play a comparatively minor role.     

Analytical description of the activation of multi-state receptors by continuous neurotransmitter signals at brain synapses.

Chemical synaptic transmission is a fundamental component of interneuronal communications in the central nervous system (CNS). Discharge of a presynaptic vesicle containing a few thousand molecules (a quantum) of neurotransmitter into the synaptic cleft generates a transmitter concentration signal that drives postsynaptic ion-channel receptors. These receptors exhibit multiple states, with state transition kinetics dependent on neurotransmitter concentration. Here, a novel and simple analytical approach for describing gating of multi-state receptors by signals with complex continuous time courses is used to describe the generation of glutamate-mediated quantal postsynaptic responses at brain synapses. The neurotransmitter signal, experienced by multi-state N-methyl-D-aspartate (NMDA)- and L-alpha-amino-3-hydroxy-5-methyl-4-isoxazolepropionate (AMPA)-type glutamate receptors at specific points in a synaptic cleft, is approximated by a series of step functions of different intensity and duration and used to drive a Markovian, multi-state kinetic scheme that describes receptor gating. Occupancy vectors at any point in time can be computed interatively from the occupancy vectors at the times of steps in transmitter concentration. Multi-state kinetic schemes for both the low-affinity AMPA subtype of glutamate receptor and for the high-affinity NMDA subtype are considered, and expected NMDA and AMPA components of synaptic currents are calculated. The amplitude of quantal responses mediated by postsynaptic receptor clusters having specific spatial distributions relative to foci of quantal neurotransmitter release is then calculated and related to the displacement between the center of the postsynaptic receptor cluster and the focus of synaptic vesicle discharge. Using this approach we show that the spatial relation between the focus of release and the center of the postsynaptic receptor cluster affects synaptic efficacy. We also show how variation in this relation contributes to variation in synaptic current amplitudes.     

Monte Carlo Simulation of Release of Vesicular Content in Neuroendocrine Cells

The release of transmitter from the vesicle, its diffusion through the fusion pore, and the cleft and its interaction with the carbon electrode were simulated using the Monte Carlo method. According to the simulation the transmitter release is largely determined by geometric factors – the ratio of the fusion pore cross-sectional and vesicular areas, if the diffusion constant is as in the aqueous solution – but the speed of transmitter dissociation from the gel matrix plays an important role during the rise phase of release. Transmitter is not depleted near the entrance to the fusion pore and there is no cleft-to-vesicle feedback, but the depletion becomes evident if the diffusion constant is reduced, especially if the pore is wide. In general, the time course of amperometric currents closely resembles the time course of the simulated transmitter concentration in the cleft and the time course of release. Surprisingly, even a tenfold change of the electrode efficiency has only a marginal effect on the amplitude or the time course of amperometric currents. Greater electrode efficiency however lowers the cleft concentration, but only if the cleft is narrow. As the cleft widens the current amplitudes diminish and rise times lengthen, but the decay times are less affected. Moreover, the amplitude dependence of the rise and decay times becomes steeper as the cleft widens and/or as the release kinetics slows. Finally, lower diffusion constant of transmitter in the narrow cleft does not further prolong the amperometric currents, whose slow time course reflects slow release kinetics.     

Asymmetry of glia near central synapses favors presynaptically directed glutamate escape

Recent findings demonstrate that synaptically released excitatory neurotransmitter glutamate activates receptors outside the immediate synaptic cleft and that the extent of such extrasynaptic actions is regulated by the high affinity glutamate uptake. The bulk of glutamate transporter systems are evenly distributed in the synaptic neuropil, and it is generally assumed that glutamate escaping the cleft affects pre- and postsynaptic receptors to a similar degree. To test whether this is indeed the case, we use quantitative electron microscopy and establish the stochastic pattern of glial occurrence in the three-dimensional (3D) vicinity of two common types of excitatory central synapses, stratum radiatum synapses in hippocampus and parallel fiber synapses in cerebellum. We find that the occurrence of glia postsynaptically is strikingly higher (3-4-fold) than presynaptically, in both types of synapses. To address the functional consequences of this asymmetry, we simulate diffusion and transport of synaptically released glutamate in these two brain areas using a detailed 3D compartmental model of the extracellular space with glutamate transporters arranged unevenly, in accordance with the obtained experimental data. The results predict that glutamate escaping the synaptic cleft is 2-4 times more likely to activate presynaptic compared to postsynaptic receptors. Simulations also show that postsynaptic neuronal transporters (EAAT4 type) at dendritic spines of cerebellar Purkinje cells exaggerate this asymmetry further. Our data suggest that the perisynaptic environment of these common central synapses favors fast presynaptic feedback in the information flow while preserving the specificity of the postsynaptic input.     

The effects of geometrical parameters on synaptic transmission: a Monte Carlo simulation study

Monte Carlo simulations of transmitter diffusion and its interactions with postsynaptic receptors have been used to study properties of quantal responses at central synapses. Fast synaptic responses characteristic of those recorded at glycinergic junctions on the teleost Mauthner cell (time to peak approximately 0.3-0.4 ms and decay time constant approximately 3-6 ms) served as the initial reference, and smaller contacts with fewer postsynaptic receptors were also modeled. Consistent with experimental findings, diffusion, simulated using a random walk algorithm and assuming a diffusion coefficient of 0.5-1.0 x 10(-5) cm2 s(-1), was sufficiently fast to account for transmitter removal from the synaptic cleft. Transmitter-receptor interactions were modeled as a two-step binding process, with the double-bound state having opened and closed conformations. Addition of a third binding step only slightly decreased response amplitude but significantly slowed both its rising and decay phases. The model allowed us to assess the sources of response variability and the likelihood of postsynaptic saturation as functions of multiple kinetic and spatial parameters. The method of nonstationary fluctuation analysis, typically used to estimate the number of functional channels at a synapse and single channel current, proved unreliable, presumably because the receptors in the postsynaptic matrix are not uniformly exposed to the same profile of transmitter concentration. Thus, the time course of the probability of channel opening most likely varies among receptors. Finally, possible substrates for phenomena of synaptic plasticity, such as long-term potentiation, were explored, including the diameter of the contact zone, defined by the region of pre- and postsynaptic apposition, the number and distribution of the receptors, and the degree of vesicle filling. Surprisingly, response amplitude is quite sensitive to the size of the receptor-free annulus surrounding the receptor cluster, such that expansion of the contact zone could produce an appreciable increase in quantal size, normally attributed to either the presence of more receptors or the release of more transmitter molecules.     

Variability of neurotransmitter concentration and nonsaturation of postsynaptic AMPA receptors at synapses in hippocampal cultures and slices

To understand the elementary unit of synaptic communication between CNS neurons, one must know what causes the variability of quantal postsynaptic currents and whether unitary packets of transmitter saturate postsynaptic receptors. We studied single excitatory synapses between hippocampal neurons in culture. Focal glutamate application at individual postsynaptic sites evoked currents (I(glu)) with little variability compared with quantal excitatory postsynaptic currents (EPSCs). The maximal I(glu) was >2-fold larger than the median EPSC. Thus, variations in [glu]cleft are the main source of variability in EPSC size, and glutamate receptors are generally far from saturation during quantal transmission. This conclusion was verified by molecular antagonism experiments in hippocampal cultures and slices. The general lack of glutamate receptor saturation leaves room for increases in [glu]cleft as a mechanism for synaptic plasticity.     

A Monte Carlo model reveals independent signaling at central glutamatergic synapses

We have developed a biophysically realistic model of receptor activation at an idealized central glutamatergic synapse that uses Monte Carlo techniques to simulate the stochastic nature of transmission following release of a single synaptic vesicle. For the a synapse with 80 AMPA and 20 NMDA receptors, a single quantum, with 3000 glutamate molecules, opened approximately 3 NMDARs and 20 AMPARs. The number of open receptors varied directly with the total number of receptors, and the fraction of open receptors did not depend on the ratio of co-localized AMPARs and NMDARs. Variability decreased with increases in either total receptor number or quantal size, and differences between the variability of AMPAR and NMDAR responses were due solely to unequal numbers of receptors at the synapse. Despite NMDARs having a much higher affinity for glutamate than AMPARs, quantal release resulted in similar occupancy levels in both receptor types. Receptor activation increased with number of transmitter molecules released or total receptor number, whereas occupancy levels were only dependent on quantal size. Tortuous diffusion spaces reduced the extent of spillover and the activation of extrasynaptic receptors. These results support the conclusion that signaling is spatially independent within and between central glutamatergic synapses.     

Extracellular glutamate diffusion determines the occupancy of glutamate receptors at CA1 synapses in the hippocampus

Following exocytosis at excitatory synapses in the brain, glutamate binds to several subtypes of postsynaptic receptors. The degree of occupancy of AMPA and NMDA receptors at hippocampal synapses is, however, not known. One approach to estimate receptor occupancy is to examine quantal amplitude fluctuations of postsynaptic signals in hippocampal neurons studied in vitro. The results of such experiments suggest that NMDA receptors at CA1 synapses are activated not only by glutamate released from the immediately apposed presynaptic terminals, but also by glutamate spillover from neighbouring terminals. Numerical simulations point to the extracellular diffusion coefficient as a critical parameter that determines the extent of activation of receptors positioned at different distances from the release site. We have shown that raising the viscosity of the extracellular medium can modulate the diffusion coefficient, providing an experimental tool to investigate the role of diffusion in activation of synaptic and extrasynaptic receptors. Whether intersynaptic cross-talk mediated by NMDA receptors occurs in vivo remains to be determined. The theoretical and experimental approaches described here also promise to shed light on the roles of metabotropic and kainate receptors, which often occur in an extrasynaptic distribution, and are therefore positioned to sense glutamate escaping from the synaptic cleft.     

Agrin and Synaptic Laminin Are Required to Maintain Adult Neuromuscular Junctions

As synapses form and mature the synaptic partners produce organizing molecules that regulate each other’s differentiation and ensure precise apposition of pre- and post-synaptic specializations. At the skeletal neuromuscular junction (NMJ), these molecules include agrin, a nerve-derived organizer of postsynaptic differentiation, and synaptic laminins, muscle-derived organizers of presynaptic differentiation. Both become concentrated in the synaptic cleft as the NMJ develops and are retained in adulthood. Here, we used mutant mice to ask whether these organizers are also required for synaptic maintenance. Deletion of agrin from a subset of adult motor neurons resulted in the loss of acetylcholine receptors and other components of the postsynaptic apparatus and synaptic cleft. Nerve terminals also atrophied and eventually withdrew from muscle fibers. On the other hand, mice lacking the presynaptic organizer laminin-α4 retained most of the synaptic cleft components but exhibited synaptic alterations reminiscent of those observed in aged animals. Although we detected no marked decrease in laminin or agrin levels at aged NMJs, we observed alterations in the distribution and organization of these synaptic cleft components suggesting that such changes could contribute to age-related synaptic disassembly. Together, these results demonstrate that pre- and post-synaptic organizers actively function to maintain the structure and function of adult NMJs.     

Dynamics of postsynaptic glutamate receptor targeting

Targeting of glutamate receptors to synapses is an important event in both developing and mature neurons. Glutamate receptors are delivered to nascent synapses during synaptogenesis and to existing synapses during activity-dependent synaptic strengthening. Increasing evidence suggests that glutamate receptors are inserted into the plasma membrane before they accumulate at the synapse. Lateral diffusion of receptors occurs at both synaptic and non-synaptic membranes, and glutamate receptors can exchange rapidly between synaptic and extrasynaptic sites. In addition, recent studies show that postsynaptic scaffold molecules can be highly mobile. The dynamic nature of the synapse suggests that many mechanisms might be involved in regulating synapse formation and synaptic plasticity.     

Presynaptic control of quantal size: kinetic mechanisms and implications for synaptic transmission and plasticity

Although the strength of quantal synaptic transmission is jointly controlled by pre- and post-synaptic mechanisms, the presynaptic mechanisms remain substantially less well characterized. Recent studies reveal that a single package of neurotransmitter is generally insufficient to activate all available postsynaptic receptors, whereas the sum of transmitter from multiple vesicles can result in receptor saturation. Thus, depending upon the number of vesicles released, a given synaptic pathway might be either 'reliable' or 'unreliable'. A lack of receptor saturation in turn makes it possible to modify quantal size by altering the flux of transmitter through the synaptic cleft. Studies are now illuminating several new mechanisms behind the regulation of this transmitter flux--characteristics that control how transmitter is loaded into vesicles, how it is released and the manner by which it interacts with postsynaptic receptors.     

Transmitter concentration profiles in the synaptic cleft: an analytical model of release and diffusion.

A three-dimensional model for release and diffusion of glutamate in the synaptic cleft was developed and solved analytically. The model consists of a source function describing transmitter release from the vesicle and a diffusion function describing the spread of transmitter in the cleft. Concentration profiles of transmitter at the postsynaptic side were calculated for different transmitter concentrations in a vesicle, release scenarios, and diffusion coefficients. From the concentration profiles the receptor occupancy could be determined using alpha-amino-3-hydroxy-5-methyl-4-isoxazolepropionic acid receptor kinetics. It turned out that saturation of receptors and sufficiently fast currents could only be obtained if the diffusion coefficient was one order of magnitude lower than generally assumed, and if the postsynaptic receptors formed clusters with a diameter of roughly 100 nm directly opposite the release sites. Under these circumstances the gradient of the transmitter concentration at the postsynaptic membrane outside the receptor clusters was steep, with minimal cross-talk among neighboring receptor clusters. These findings suggest that for each release site a corresponding receptor aggregate exists, subdividing an individual synapse into independent functional subunits without the need for specific lateral diffusion barriers.     

Kinetics of Synaptic Transmission at Ribbon Synapses of Rods and Cones

The ribbon synapse is a specialized structure that allows photoreceptors to sustain the continuous release of vesicles for hours upon hours and years upon years but also respond rapidly to momentary changes in illumination. Light responses of cones are faster than those of rods and, mirroring this difference, synaptic transmission from cones is also faster than transmission from rods. This review evaluates the various factors that regulate synaptic kinetics and contribute to kinetic differences between rod and cone synapses. Presynaptically, the release of glutamate-laden synaptic vesicles is regulated by properties of the synaptic proteins involved in exocytosis, influx of calcium through calcium channels, calcium release from intracellular stores, diffusion of calcium to the release site, calcium buffering, and extrusion of calcium from the cytoplasm. The rate of vesicle replenishment also limits the ability of the synapse to follow changes in release. Post-synaptic factors include properties of glutamate receptors, dynamics of glutamate diffusion through the cleft, and glutamate uptake by glutamate transporters. Thus, multiple synaptic mechanisms help to shape the responses of second-order horizontal and bipolar cells.