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

Computational models of neurotransmission at cerebellar synapses unveil the impact on network computation

The neuroscientific field benefits from the conjoint evolution of experimental and computational techniques, allowing for the reconstruction and simulation of complex models of neurons and synapses. Chemical synapses are characterized by presynaptic vesicle cycling, neurotransmitter diffusion, and postsynaptic receptor activation, which eventually lead to postsynaptic currents and subsequent membrane potential changes. These mechanisms have been accurately modeled for different synapses and receptor types (AMPA, NMDA, and GABA) of the cerebellar cortical network, allowing simulation of their impact on computation. Of special relevance is short-term synaptic plasticity, which generates spatiotemporal filtering in local microcircuits and controls burst transmission and information flow through the network. Here, we present how data-driven computational models recapitulate the properties of neurotransmission at cerebellar synapses. The simulation of microcircuit models is starting to reveal how diverse synaptic mechanisms shape the spatiotemporal profiles of circuit activity and computation.

Computational modeling of trans-synaptic nanocolumns, a modulator of synaptic transmission

Understanding synaptic transmission is of crucial importance in neuroscience. The spatial organization of receptors, vesicle release properties and neurotransmitter molecule diffusion can strongly influence features of synaptic currents. Newly discovered structures coined trans-synaptic nanocolumns were shown to align presynaptic vesicles release sites and postsynaptic receptors. However, how these structures, spanning a few tens of nanometers, shape synaptic signaling remains little understood. Given the difficulty to probe submicroscopic structures experimentally, computer modeling is a useful approach to investigate the possible functional impacts and role of nanocolumns. In our in silico model, as has been experimentally observed, a nanocolumn is characterized by a tight distribution of postsynaptic receptors aligned with the presynaptic vesicle release site and by the presence of trans-synaptic molecules which can modulate neurotransmitter molecule diffusion. In this work, we found that nanocolumns can play an important role in reinforcing synaptic current mostly when the presynaptic vesicle contains a small number of neurotransmitter molecules. Our work proposes a new methodology to investigate in silico how the existence of trans-synaptic nanocolumns, the nanometric organization of the synapse and the lateral diffusion of receptors shape the features of the synaptic current such as its amplitude and kinetics.

Heterogeneity of glutamatergic synapses: cellular mechanisms and network consequences

Chemical synapses are commonly known as a structurally and functionally highly diverse class of cell-cell contacts specialized to mediate communication between neurons. They represent the smallest "computational" unit of the brain and are typically divided into excitatory and inhibitory as well as modulatory categories. These categories are subdivided into diverse types, each representing a different structure-function repertoire that in turn are thought to endow neuronal networks with distinct computational properties. The diversity of structure and function found among a given category of synapses is referred to as heterogeneity. The main building blocks for this heterogeneity are synaptic vesicles, the active zone, the synaptic cleft, the postsynaptic density, and glial processes associated with the synapse. Each of these five structural modules entails a distinct repertoire of functions, and their combination specifies the range of functional heterogeneity at mammalian excitatory synapses, which are the focus of this review. We describe synapse heterogeneity that is manifested on different levels of complexity ranging from the cellular morphology of the pre- and postsynaptic cells toward the expression of different protein isoforms at individual release sites. We attempt to define the range of structural building blocks that are used to vary the basic functional repertoire of excitatory synaptic contacts and discuss sources and general mechanisms of synapse heterogeneity. Finally, we explore the possible impact of synapse heterogeneity on neuronal network function.

The glutamatergic synapse: a complex machinery for information processing

Being the most abundant synaptic type, the glutamatergic synapse is responsible for the larger part of the brain's information processing. Despite the conceptual simplicity of the basic mechanism of synaptic transmission, the glutamatergic synapse shows a large variation in the response to the presynaptic release of the neurotransmitter. This variability is observed not only among different synapses but also in the same single synapse. The synaptic response variability is due to several mechanisms of control of the information transferred among the neurons and suggests that the glutamatergic synapse is not a simple bridge for the transfer of information but plays an important role in its elaboration and management. The control of the synaptic information is operated at pre, post, and extrasynaptic sites in a sort of cooperation between the pre and postsynaptic neurons which also involves the activity of other neurons. The interaction between the different mechanisms of control is extremely complicated and its complete functionality is far from being fully understood. The present review, although not exhaustively, is intended to outline the most important of these mechanisms and their complexity, the understanding of which will be among the most intriguing challenges of future neuroscience.

Synaptic environment and extrasynaptic glutamate signals: The quest continues

Behaviour of a mammal relies on the brain's excitatory circuits equipped with glutamatergic synapses. In most cases, glutamate escaping from the synaptic cleft is rapidly buffered and taken up by high-affinity transporters expressed by nearby perisynaptic astroglial processes (PAPs). The spatial relationship between glutamatergic synapses and PAPs thus plays a crucial role in understanding glutamate signalling actions, yet its intricate features can only be fully appreciated using methods that operate beyond the diffraction limit of light. Here, we examine principal aspects pertaining to the receptor actions of glutamate, inside and outside the synaptic cleft in the brain, where the organisation of synaptic micro-physiology and micro-environment play a critical part. In what conditions and how far glutamate can escape the synaptic cleft activating its target receptors outside the immediate synapse has long been the subject of debate. Evidence is also emerging that neuronal activity- and astroglia-dependent glutamate spillover actions could be important across the spectrum of cognitive functions This article is part of the special issue on 'Glutamate Receptors - The Glutamatergic Synapse'.

Quantification of trans-synaptic protein alignment: A data analysis case for single-molecule localization microscopy

Nanoscale distribution of proteins and their relative positioning within a defined subcellular region are key to their physiological functions. Thanks to the super-resolution imaging methods, especially single-molecule localization microscopy (SMLM), mapping the three-dimensional distribution of multiple proteins has been easier and more efficient than ever. Nevertheless, in spite of the many tools available for efficient localization detection and image rendering, it has been a challenge to quantitatively analyze the 3D distribution and relative positioning of proteins in these SMLM data. Here, using heterogeneously distributed synaptic proteins as examples, we describe in detail a series of analytical methods including detection of nanoscale density clusters, quantification of the trans-synaptic alignment between these protein densities, and automatic en face projection and averaging. These analyses were performed within customized Matlab routines and we make the full scripts available. The concepts behind these analytical methods and the scripts can be adapted for quantitative analysis of spatial organization of other macromolecular complexes.

Role of Heterogeneous Macromolecular Crowding and Geometrical Irregularity at Central Excitatory Synapses in Shaping Synaptic Transmission

Besides the geometrical tortousity due to the extrasynaptic structures, macromolecular crowding and geometrical irregularities constituting the cleft composition at central excitatory synapses has a major and direct role in retarding the glutamate diffusion within the cleft space. However, the cleft composition may not only coarsely reduce the overall diffusivity of the glutamate but may also lead to substantial spatial variation in the diffusivity across the cleft space. Decrease in the overall diffusivity of the glutamate may have straightforward consequences to the glutamate transients in the cleft. However, how spatial variation in the diffusivity may further affect glutamate transients is an intriguing aspect. Therefore, to understand the role of cleft heterogeneity, the present study adopts a novel approach of glutamate diffusion which considers a gamma statistical distribution of the diffusion coefficient of glutamate (Dglut) across the cleft space, such that its moments discernibly capture the dual impacts of the cleft composition, and further applies the framework of superstatistics. The findings reveal a power law behavior in the glutamate transients, akin to the long-range anomalous subdiffusion, which leads to slower decay profile of cleft glutamate at higher intensity of cleft heterogeneity. Moreover, increase in the cleft heterogeneity is seen to eventually cause slower-rising excitatory postsynaptic currents with higher amplitudes, lesser noise, and prolonged duration of charge transfer across the postsynaptic membrane. Further, with regard to the conventional standard diffusion approach, the study suggests that the effective Dglut essentially derives from the median of the Dglut distribution and does not necessarily need to be the mean Dglut. Together, the findings indicate a strong implication of cleft heterogeneity to the metabolically cost-effective tuning of synaptic response during the phenomenon of plasticity at individual synapses and also provide an additional factor of variability in transmission across identical synapses.

Efficient integration of synaptic events by NMDA receptors in three-dimensional neuropil

Sustained activation of NMDA receptors (NMDARs) plays an important role in controlling activity of neural circuits in the brain. However, whether this activation reflects the ambient level of excitatory neurotransmitter glutamate in brain tissue or whether it depends mainly on local synaptic discharges remains poorly understood. To shed light on the underlying biophysics here we developed and explored a detailed Monte Carlo model of a realistic three-dimensional neuropil fragment containing 54 excitatory synapses. To trace individual molecules and their individual receptor interactions on this scale, we have designed and implemented a dedicated computer cluster and the appropriate software environment. Our simulations have suggested that sparse synaptic discharges are 20-30 times more efficient than nonsynaptic (stationary, leaky) supply of glutamate in controlling sustained NMDAR occupancy in the brain. This mechanism could explain how the brain circuits provide substantial background activation of NMDARs while maintaining a negligible ambient glutamate level in the extracellular space. Thus the background NMDAR occupancy, rather than the background glutamate level, is likely to reflect the ongoing activity in local excitatory networks.

Model of an excitatory synapse based on stochastic processes

We present a mathematical model of a biological synapse based on stochastic processes to establish the temporal behavior of the postsynaptic potential following a quantal synaptic transmission. This potential form is the basis of the neural code. We suppose that the release of neurotransmitters in the synaptic cleft follows a Poisson process, and that they diffuse according to integrated Ornstein-Uhlenbeck processes in 3-D with random initial positions and velocities. The diffusion occurs in an isotropic environment between two infinite parallel planes representing the pre- and postsynaptic membrane. We state that the presynaptic membrane is perfectly reflecting and that the other is perfectly absorbing. The activation of the receptors polarizes the postsynaptic membrane according to a parallel RC circuit scheme. We present the results obtained by simulations according to a Gillespie algorithm and we show that our model exhibits realistic postsynaptic behaviors from a simple quantal occurrence.

Swim pacemaker response to bath applied neurotransmitters in the cubozoan Tripedalia cystophora

The four rhopalia of cubomedusae are integrated parts of the central nervous system carrying their many eyes and thought to be the centres of visual information processing. Rhopalial pacemakers control locomotion through a complex neural signal transmitted to the ring nerve and the signal frequency is modulated by the visual input. Since electrical synapses have never been found in the cubozoan nervous system all signals are thought to be transmitted across chemical synapses, and so far information about the neurotransmitters involved are based on immunocytochemical or behavioural data. Here we present the first direct physiological evidence for the types of neurotransmitters involved in sensory information processing in the rhopalial nervous system. FMRFamide, serotonin and dopamine are shown to have inhibitory effect on the pacemaker frequency. There are some indications that the fast acting acetylcholine and glycine have an initial effect and then rapidly desensitise. Other tested neuroactive compounds (GABA, glutamate, and taurine) could not be shown to have a significant effect.

Impact of synaptic neurotransmitter concentration time course on the kinetics and pharmacological modulation of inhibitory synaptic currents

The time course of synaptic currents is a crucial determinant of rapid signaling between neurons. Traditionally, the mechanisms underlying the shape of synaptic signals are classified as pre- and post-synaptic. Over the last two decades, an extensive body of evidence indicated that synaptic signals are critically shaped by the neurotransmitter time course which encompasses several phenomena including pre- and post-synaptic ones. The agonist transient depends on neurotransmitter release mechanisms, diffusion within the synaptic cleft, spill-over to the extra-synaptic space, uptake, and binding to post-synaptic receptors. Most estimates indicate that the neurotransmitter transient is very brief, lasting between one hundred up to several hundreds of microseconds, implying that post-synaptic activation is characterized by a high degree of non-equilibrium. Moreover, pharmacological studies provide evidence that the kinetics of agonist transient plays a crucial role in setting the susceptibility of synaptic currents to modulation by a variety of compounds of physiological or clinical relevance. More recently, the role of the neurotransmitter time course has been emphasized by studies carried out on brain slice models that revealed a striking, cell-dependent variability of synaptic agonist waveforms ranging from rapid pulses to slow volume transmission. In the present paper we review the advances on studies addressing the impact of synaptic neurotransmitter transient on kinetics and pharmacological modulation of synaptic currents at inhibitory synapses.

Molecular diffusion model of neurotransmitter homeostasis around synapses supporting gradients

Neurotransmitter homeostasis in and around a synapse involves complex random processes such as diffusion, molecular binding, and uptake by glial transporters. A three-dimensional stochastic diffusion model of a synapse was developed to provide molecular-level details of neurotransmitter homeostasis not predicted by alternative models based on continuum approaches. The development was illustrated through an example case cortico-accumbens synapse that successfully integrated neuroadaptations observed after chronic cocaine. By incorporating cystine-glutamate exchanger as a nonsynaptic release site for glutamate, the stochastic model was used to quantify the relative contributions of synaptic and nonsynaptic sources to extracellular concentration and to estimate molecular influx rates into the perisynapse. A perturbation analysis showed that among the parameters considered, variation in surface density of glial transporters had the largest effect on glutamate concentrations. The stochastic diffusion model of the example synapse was further generalized to characterize glial morphology by studying the role of diffusion path length in supporting neurotransmitter gradients and isolating the synapse. For the same set of parameters, diffusion path length was found to be proportional to the gradient supported.

Realistic modelling of receptor activation in hippocampal excitatory synapses: analysis of multivesicular release, release location, temperature and synaptic cross-talk

Chemically mediated synaptic transmission results from fusion of synaptic vesicles with the presynaptic plasma membrane, subsequent release of the vesicular content into the cleft and binding to postsynaptic receptors. Previous modelling studies of excitatory neurotransmitter glutamate were based on simplified geometries failing to account for the biologically realistic synaptic environment, in particular, the presence of astrocytes, the geometry of extracellular space, and the neurotransmitter uptake mechanism. Using 3-dimensional reconstructions of hippocampal glutamatergic synapses including the surrounding astrocytic processes we have developed a biologically realistic model to analyse receptor activation in different conditions. We used the finite element method to simulate glutamate release, analyse glutamate diffusion following single and multiple vesicle release and binding at the postsynaptic site to AMPA and NMDA receptors. We demonstrate that: (1) the transmitter diffusion is highly temperature-sensitive; (2) release conditions and geometry more specifically affect AMPARs than NMDARs; (3) the sensitivities of AMPARs and NMDARs to simultaneous vesicular release are different; (4) in the case of multivesicle neurotransmitter release with variable delays, the binding of glutamate to AMPARs is additive up to 1 ms after the release, then becomes independent, but to NMDARs the binding is additive up to 33 ms; (5) the number of AMPARs varies more than the number of NMDRs in response to the input firing patterns; (6) the presence of astrocytes effectively blocks synaptic cross-talk; and (7) synaptic cross-talk, mediated by NMDARs but not AMPARs, is only possible after quasi-simultaneous multivesicular release at physiological temperature (35 degrees C) without intervening astrocytes, but not at 25 degrees C. Our simulations demonstrate the importance of temperature and ultrastructural synaptic environment in synaptic transmission and synaptic cross-talk.

Computational model of extracellular glutamate in the nucleus accumbens incorporates neuroadaptations by chronic cocaine

Chronic cocaine administration causes instability in extracellular glutamate in the nucleus accumbens that is thought to contribute to the vulnerability to relapse. A computational framework was developed to model glutamate in the extracellular space, including synaptic and nonsynaptic glutamate release, glutamate elimination by glutamate transporters and diffusion, and negative feedback on synaptic release via metabotropic glutamate receptors (mGluR2/3). This framework was used to optimize the geometry of the glial sheath surrounding excitatory synapses, and by inserting physiological values, accounted for known stable extracellular, extrasynaptic concentrations of glutamate measured by microdialysis and glutamatergic tone on mGluR2/3. By using experimental values for cocaine-induced reductions in cystine-glutamate exchange and mGluR2/3 signaling, and by predicting the down-regulation of glutamate transporters, the computational model successfully represented the experimentally observed increase in glutamate that is seen in rats during cocaine-seeking. This model provides a mathematical framework for describing how pharmacological or pathological conditions influence glutamate transmission measured by microdialysis.

The undertow of sulfur metabolism on glutamatergic neurotransmission

Metabolic interdependence between specialized cells in an organ represents a strategy for energy economy by requiring expression of only a subset of pathway genes in a given cell type. In brain, sulfur metabolism exemplifies this principle of metabolic cooperation between glial and neuronal cells and furnishes three key reagents: S-adenosylmethionine, glutathione and taurine. The pathways for glutathione and taurine syntheses depend on metabolic integration between astrocytes and neurons and intersect with the glutamine-glutamate cycle, which underlies glutamatergic synaptic transmission and requires cooperation between these cell types. We propose that underlying waves of glutamate clearance by astrocytes are activation of cystine import and taurine efflux that result, respectively, from a shared transporter and an increase in solute concentration that triggers osmoregulatory responses.

Parametric spectral analysis of nonstationary fluctuations of excitatory synaptic currents

We assessed on Monte-Carlo simulated excitatory post-synaptic currents the ability of autoregressive (AR)-model fitting to evaluate their fluctuations. AR-model fitting consists of a linear filter describing the process that generates the fluctuations when driven with a white noise. Its fluctuations provide a filtered version of the signal and have a spectral density depending on the properties of the linear filter. When the spectra of the non-stationary fluctuations of excitatory post-synaptic currents were estimated by fitting AR-models to the segments of current fluctuations, assumed to be stationary and independent, the parameter and spectral estimates were scattered. The scatter was much reduced if the time-variant AR-models were fitted using stochastic adaptive estimators (Kalman, recursive least squares and least mean squares). The ability of time-variant AR-models to accurately fit the current fluctuations was monitored by comparing the fluctuations with predicted fluctuations, and by evaluating the model-learning rate. The median frequency of current fluctuations, which could be rapidly tracked and estimated from the individual quantal events (either Monte-Carlo simulated or recorded from pyramidal neurons of rat hippocampus), rose during the rise phase, before declining to a lower steady-state level during the decay phase of quantal event, whereas the variance showed a broad peak. The closing rate of AMPA channels directly affects the steady-state median frequency, whereas the transient peak can be modulated by a variety of factors-number of molecules released, ability of glutamate molecules to re-enter the synaptic cleft, diffusion constant of glutamate in the cleft and opening rate of AMPA channels. In each case, the effect on the amplitude and decay time of mEPSCs and on the current fluctuations differs. Each factor thus leaves its own kinetic fingerprint arguing that the contribution of such factors can be inferred from the combined kinetic properties of individual mEPSCs.

Computational modeling of paroxysmal depolarization shifts in neurons induced by the glutamate release from astrocytes

Recent experimental studies have shown that astrocytes respond to external stimuli with a transient increase of the intracellular calcium concentration or can exhibit self-sustained spontaneous activity. Both evoked and spontaneous astrocytic calcium oscillations are accompanied by exocytosis of glutamate caged in astrocytes leading to paroxysmal depolarization shifts (PDS) in neighboring neurons. Here, we present a simple mathematical model of the interaction between astrocytes and neurons that is able to numerically reproduce the experimental results concerning the initiation of the PDS. The timing of glutamate release from the astrocyte is studied by means of a combined modeling of a vesicle cycle and the dynamics of SNARE-proteins. The neuronal slow inward currents (SICs), induced by the astrocytic glutamate and leading to PDS, are modeled via the activation of presynaptic glutamate receptors. The dependence of the bidirectional communication between neurons and astrocytes on the concentration of glutamate transporters is analyzed, as well. Our numerical results are in line with experimental findings showing that astrocyte can induce synchronous PDSs in neighboring neurons, resulting in a transient synchronous spiking activity.

Electrophysiological description of mechanisms determining synaptic transmission and its modulation

Signal integration in neurons is a complex process that depends on e.g. the kinetics of synaptic currents, distribution of synaptic connections as well as passive and excitatory membrane properties. The time course of synaptic currents is largely determined by the kinetics of the postsynaptic receptors and the time course of synaptic neurotransmitter concentration. The analysis of current responses to rapid agonist applications provides the means to study the ligand-gated receptor gating but experimentally based estimation of neurotransmitter transient at central synapses was an important challenge during the last decade. Both theoretical as well as experimentally based approaches indicated that synaptic agonist transient is very brief, implying that the activation of postsynaptic receptors occurs in conditions of extreme non-equilibrium. Such a dynamic pattern of activation of postsynaptic receptors has a crucial impact not only on the kinetics of synaptic currents but also on their susceptibility to pharmacological modulation.

Long-term exposure to histamine induces the expression of an embryonic-like motor pattern in an adult nervous system

Neuromodulatory inputs play important roles in shaping the outputs of neural networks. While the actions of neuromodulatory substances over the short term (seconds, minutes) have been examined in detail, far less is known about the possible longer-term (hours) effects of these substances. To investigate this issue, we used the stomatogastric nervous system (STNS) of the lobster to examine the short- and long-term effects of histamine on rhythmic network activity. The application of histamine to the entire STNS had strong inhibitory effects on all three of the STNS networks, observable within minutes. In contrast, longer-term (> 1 h) application of histamine induced the expression of a single, unified rhythm involving neurons from all three networks. Selective application of histamine to different regions of the STNS demonstrated that a unified rhythm arises following the long-term application of histamine to the commissural ganglia (CoGs; modulatory centres), but not the stomatogastric ganglion (site of neural networks). Strikingly, the single rhythm observed following the long-term application of histamine to the CoGs exhibits many similarities with the single rhythm expressed by the embryonic STNS. Together, these results demonstrate that histamine has markedly different short- and long-term effects on network activity; short-term effects arising through direct actions on the networks and long-term effects mediated by actions on modulatory neurons. Furthermore, they indicate that histamine is able to induce the expression of an embryonic-like rhythm in an adult system, suggesting that long-term actions of histamine may play key roles in the development of the STNS networks.

Glutamate-induced exocytosis of glutamate from astrocytes

Recent studies indicate that astrocytes can play a much more active role in neuronal circuits than previously believed, by releasing neurotransmitters such as glutamate and ATP. Here we report that local application of glutamate or glutamine synthetase inhibitors induces astrocytic release of glutamate, which activates a slowly decaying transient inward current (SIC) in CA1 pyramidal neurons and a transient inward current in astrocytes in hippocampal slices. The occurrence of SICs was accompanied by an appearance of large vesicles around the puffing pipette. The frequency of SICs was positively correlated with [glutamate]o. EM imaging of anti-glial fibrillary acid protein-labeled astrocytes showed glutamate-induced large astrocytic vesicles. Imaging of FM 1-43 fluorescence using two-photon laser scanning microscopy detected glutamate-induced formation and fusion of large vesicles identified as FM 1-43-negative structures. Fusion of large vesicles, monitored by collapse of vesicles with a high intensity FM 1-43 stain in the vesicular membrane, coincided with SICs. Glutamate induced two types of large vesicles with high and low intravesicular [Ca2+]. The high [Ca2+] vesicle plays a major role in astrocytic release of glutamate. Vesicular fusion was blocked by infusing the Ca2+ chelator, 1,2-bis(2-aminophenoxy)ethane-N,N,N',N'-tetraacetic acid, or the SNARE blocker, tetanus toxin, suggesting Ca2+- and SNARE-dependent fusion. Infusion of the vesicular glutamate transport inhibitor, Rose Bengal, reduced astrocytic glutamate release, suggesting the involvement of vesicular glutamate transports in vesicular transport of glutamate. Our results demonstrate that local [glutamate]o increases induce formation and exocytotic fusion of glutamate-containing large astrocytic vesicles. These large vesicles could play important roles in the feedback control of neuronal circuits and epileptic seizures.

Modeling of slow glutamate diffusion and AMPA receptor activation in the cerebellar glomerulus

Synaptic conductances are influenced markedly by the geometry of the space surrounding the synapse since the transient glutamate concentration in the synaptic cleft is determined by this geometry. Our paper is an attempt to understand the reasons for slow glutamate diffusion in the cerebellar glomerulus, a structure situated around the enlarged mossy fiber terminal in the cerebellum and surrounded by a glial sheath. For this purpose, analytical expressions for glutamate diffusion in the glomerulus were considered in models with two-, three-, and fractional two-three-dimensional (2D-3D) geometry with an absorbing boundary. The time course of average glutamate concentration in the synaptic cleft of the mossy fiber-granule cell connection was calculated for both direct release of glutamate from the same synaptic unit, and for cumulative spillover of glutamate from neighboring release sites. Several kinetic schemes were examined, and the parameters of the diffusion models were estimated by identifying theoretical activation of AMPA receptors with direct release and spillover components of published experimental AMPA receptor-mediated EPSCs. For model selection, the correspondence of simulated paired-pulse ratio and EPSC increase after prevention of desensitization to experimental values were also taken into consideration. Our results suggest at least a 7- to 10-fold lower apparent diffusion coefficient of glutamate in the porous medium of the glomerulus than in water. The modeling of glutamate diffusion in the 2D-3D geometry gives the best fit of experimental EPSCs. We show that it could be only partly explained by normal diffusion of glutamate in the complex geometry of the glomerulus. We assume that anomalous diffusion of glutamate occurs in the glomerulus. A good match of experimental estimations and theoretical parameters, obtained in the simulations that use an approximation of anomalous diffusion by a solution for fractional Brownian motion, confirms our assumption.

Developmental Changes of GABA Synaptic Transient in Cerebellar Granule Cells

The time course of synaptic currents is largely determined by the microscopic gating of the postsynaptic receptors and the temporal profile of the synaptic neurotransmitter concentration. Although several lines of evidence indicate that developmental changes of GABAergic synaptic current time course are clearly correlated with a switch in postsynaptic receptors, much less is known about the modification of GABA release during development. To address this issue, we studied the sensitivity of miniature inhibitory postsynaptic currents (mIPSCs) to a quickly dissociating competitive antagonist, 1,2,5,6-tetrahydropyridine-4-yl)methylphosphinic acid (TPMPA), in neurons cultured for 6 to 8 days in vitro (DIV) ("young") and for 12 to 14 DIV ("old"). mIPSCs recorded in young neurons were significantly more resistant to the block by TPMPA. This observation was interpreted as a consequence of a more efficient displacement of TPMPA from GABA(A) receptors caused by a stronger GABA release in young neurons. The change in mIPSC sensitivity to TPMPA during development was not affected by the deletion of alpha(1) subunit, supporting its presynaptic origin. The effects of a second quickly dissociating antagonist, SR-95103 [2-(carboxy-3'-propyl)-3-amino-4-methyl-6-phenylpyridazinium chloride], on young, old, and alpha(1) -/- neurons were qualitatively the same as those obtained with TPMPA. Moreover, the analysis of current responses to ultrafast GABA applications showed that the unbinding rates of TPMPA in DIV 6 to 8 and in DIV 12 to 14 neurons are not significantly different, ruling out the postsynaptic mechanism of differential TPMPA action. Thus, we provide evidence that presynaptic GABA uniquantal release is developmentally regulated.

Presynaptic source of quantal size variability at GABAergic synapses in rat hippocampal neurons in culture

The variability of quantal size depends on both presynaptic (profile of the neurotransmitter concentration in the cleft) and postsynaptic (number and gating properties of postsynaptic receptors) factors. Here we have examined the possibility that at nonsaturated synapses in cultured hippocampal neurons, changes in both the transmitter concentration peak and its clearance from the synaptic cleft may influence the variability of spontaneous miniature synaptic GABAergic currents (mIPSCs). We found that, in contrast to the slow-off GABAA receptor antagonist bicuculline, fast-off competitive antagonists such as SR-95103 and TPMPA differentially blocked small and large mIPSCs. In the presence of flurazepam, a drug believed to increase the affinity of GABA for GABAAR, small mIPSCs were enhanced more efficiently than large events. Moreover, the addition of dextran, which increases the viscosity of the extracellular fluid, preferentially increased small mIPSCs with respect to large ones. These observations suggest that changes in the concentration peak and the speed of GABA clearance in the cleft may be an important source of synaptic variability. The study of the correlation between peak amplitude and kinetics of mIPSCs allowed determination of the relative contribution of transmitter peak concentration vs. time of GABA clearance. Small synaptic responses were associated with fast onset and decay kinetics while large amplitude currents were associated with slow kinetics, indicating a crucial role for GABA synaptic clearance in variability of mIPSCs. By using model simulations we were able to estimate the range of variability of both the concentration and the speed of clearance of the GABA transient in the synaptic cleft.

Dynamism of GABAA receptor activation shapes the “personality” of inhibitory synapses

The kinetics of synaptic currents is largely determined by the postsynaptic receptor gating and the concentration time course of synaptic neurotransmitter. While the analysis of current responses to rapid agonist application provides the means to study the ligand-gated receptor gating, no direct tools are available to measure the neurotransmitter transient at GABAergic and glutamatergic synapses. Several lines of evidence indicate that the synaptic agonist transient is very brief suggesting that the activation of postsynaptic receptors occurs in conditions of extreme non-equilibrium. Such a dynamic pattern of activation has a crucial impact not only on the kinetics of synaptic currents but also on their susceptibility to pharmacological modulation. Thus, changes in the synaptic agonist waveform due to, for example modulation of the release machinery or uptake system may considerably alter both kinetics and pharmacology of synaptic currents. The use of modifiers of GABA(A) receptor gating and low-affinity antagonists provides a tool to estimate the time course of the agonist transient revealing that synaptic neurotransmitter is not saturating and that the agonist clearance occurs at a sub-millisecond time scale. It is proposed that dynamic conditions of synaptic receptor activation assure a broad spectrum of performance rendering the synapse extremely susceptible to a variety of modulatory processes.

Properties of Quantal Transmission at CA1 Synapses

We have used Monte Carlo simulations to understand the generation of quantal responses at the single active zones of CA1 synapses. We constructed a model of AMPA channel activation that accounts for the responses to controlled glutamate application and a model of glutamate diffusion in the synaptic cleft. With no further adjustments to these models, we simulated the response to the release of glutamate from a single vesicle. The predicted response closely matches the rise time of observed responses, which recent measurements show is much faster (<100 μs) than previously thought. The simulations show that initial channel opening is driven by a brief (<100 μs) glutamate spike near the site of vesicle fusion, producing a hotspot of channel activation (diameter: ∼250 nm) smaller than many synapses. Quantal size therefore depends more strongly on the density of channels than their number, a finding that has important implications for measuring synaptic strength. Recent measurements allow estimation of AMPA receptor density at CA1 synapses. Using this value, our simulations correctly predicts a quantal amplitude of ∼10 pA. We have also analyzed the properties of excitatory postsynaptic currents (EPSCs) generated by the multivesicular release that can occur during evoked responses. We find that summation is nearly linear and that the existence of multiple narrow peaks in amplitude histograms can be accounted for. It has been unclear how to reconcile the existence of these narrow peaks, which indicate that the variation of quantal amplitude is small (CV < 0.2) with the highly variable amplitude of miniature EPSCs (mEPSCs; CV ∼ 0.6). According to one theory, mEPSC variability arises from variation in vesicle glutamate content. However, both our modeling results and recent experimental results indicate that this view cannot account for the observed rise time/amplitude correlation of mEPSCs. In contrast, this correlation and the high mEPSC variability can be accounted for if some mEPSCs are generated by two or more vesicles released with small temporal jitter. We conclude that a broad range of results can be accounted for by simple principles: quantal amplitude (∼10 pA) is stereotyped, some mEPSCs are multivesicular at moderate and large synapses, and evoked responses are generated by quasi-linear summation of multiple quanta.

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.

Wavelet analysis of nonstationary fluctuations of Monte Carlo-simulated excitatory postsynaptic currents

Tracking spectral changes of rapidly varying signals is a demanding task. In this study, we explore on Monte Carlo-simulated glutamate-activated AMPA patch and synaptic currents whether a wavelet analysis offers such a possibility. Unlike Fourier methods that determine only the frequency content of a signal, the wavelet analysis determines both the frequency and the time. This is owing to the nature of the basis functions, which are infinite for Fourier transforms (sines and cosines are infinite), but are finite for wavelet analysis (wavelets are localized waves). In agreement with previous reports, the frequency of the stationary patch current fluctuations is higher for larger currents, whereas the mean-variance plots are parabolic. The spectra of the current fluctuations and mean-variance plots are close to the theoretically predicted values. The median frequency of the synaptic and nonstationary patch currents is, however, time dependent, though at the peak of synaptic currents, the median frequency is insensitive to the number of glutamate molecules released. Such time dependence demonstrates that the "composite spectra" of the current fluctuations gathered over the whole duration of synaptic currents cannot be used to assess the mean open time or effective mean open time of AMPA channels. The current (patch or synaptic) versus median frequency plots show hysteresis. The median frequency is thus not a simple reflection of the overall receptor saturation levels and is greater during the rise phase for the same saturation level. The hysteresis is due to the higher occupancy of the doubly bound state during the rise phase and not due to the spatial spread of the saturation disk, which remains remarkably constant. Albeit time dependent, the variance of the synaptic and nonstationary patch currents can be accurately determined. Nevertheless the evaluation of the number of AMPA channels and their single current from the mean-variance plots of patch or synaptic currents is not highly accurate owing to the varying number of the activatable AMPA channels caused by desensitization. The spatial nonuniformity of open, bound, and desensitized AMPA channels, and the time dependence and spatial nonuniformity of the glutamate concentration in the synaptic cleft, further reduce the accuracy of estimates of the number of AMPA channels from synaptic currents. In conclusion, wavelet analysis of nonstationary fluctuations of patch and synaptic currents expands our ability to determine accurately the variance and frequency of current fluctuations, demonstrates the limits of applicability of techniques currently used to evaluate the single channel current and number of AMPA channels, and offers new insights into the mechanisms involved in the generation of unitary quantal events at excitatory central synapses.

Saturation in excitatory synapses of hippocampus investigated by computer simulations

The standard view of the synaptic function in excitatory synapses has been deeply questioned by recent experimental data on hippocampal glutamate synapses both for possible receptor nonsaturation and for larger and non-Gaussian peak amplitude fluctuations. Our previous investigations of the mechanisms involved in the variability of the response of hippocampal glutamatergic synapses, carried out by computer simulation of simple Brownian models of glutamate diffusion, furnished initial evidence about their presynaptic character. A new, refined model, reported here, assumes a collision volume for the glutamate molecule and a more realistic description of receptors and their binding dynamics. Based on this model, conditions for AMPA and NMDA receptor saturation have been investigated and new miniature (or quantal) EPSC parameters have been computed. The results corroborate the hypothesis that the lack of AMPA and NMDA receptor saturation and the EPSC stochastic variability are attributable to the small volume of glutamatergic synaptic vesicles and hence to the small number of glutamate molecules diffusing in the cleft after a vesicle release. The investigations better characterize some not well-known elements of the synaptic structure, such as the fusion pore, and provide useful information on AMPA receptor dynamics. Indeed, a nice fit between computed EPSCs and some miniature EPSCs in recent experimental literature allowed for the computation of new transition time values among the different AMPA receptor states through a trial-and-error optimization procedure. Moreover, the model has been used to evaluate two hypotheses on the genesis of the long-term potentiation phenomenon.

Simulation modeling of ligand receptor interactions at non-equilibrium conditions: processing of noisy inputs by ionotropic receptors

The first event in signal transduction at a synapse is the binding of transmitters to receptors. Because of rapidly changing transmitter levels this binding is unlikely to occur at equilibrium. We describe a mathematical approach that models complex receptor interactions in which the timing and amplitude of transmitter release are noisy. We show that exact solutions for simple bimolecular interactions and receptor transitions can be used to model complex reaction schemes by expressing them in sets of difference equations. Results from the difference equation method to describe binding and channel opening at extended time points compare well with standard solutions using ordinary differential equations. Because it is applicable to noisy systems we used the difference method to investigate the information processing capabilities of GABA receptors and predict how pharmacological agents may modify these properties. As previously demonstrated, the response to a single pulse of GABA is prolonged through entry into a desensitized state. During trains of stimuli the signal to noise ratio can change, and even increase progressively, but the overall transmitted fidelity of the signal decreases with increased driving frequency. The GABA modulator chlorpromazine (primarily affects agonist on and off rates) is predicated to increase receptor signal to noise ratio at all frequencies whereas pregnenolone sulfate (affects receptor desensitization) completely inhibits information transfer.

Electric fields of synaptic currents could influence diffusion of charged neurotransmitter molecules

Rapid activation of synaptic receptor-channels evokes an ion current that flows through the narrow synaptic cleft; this exerts a significant voltage drop and therefore strong electric field (10(4) V/m range) directed towards the current sinks in the cleft. To what extent this field affects fast diffusion of charged neurotransmitter molecules is not known. We draw a theoretical framework for this complex electrodiffusion phenomenon and establish the basic relationships between the synaptic current and the time course of neurotransmitter in the cleft. The analyses predict that excitatory currents could significantly accelerate the dispersion of negatively charged molecules from the cleft while attracting the positively charged molecules towards the current sinks. This previously unrecognized mechanism should affect the kinetics of synaptic receptor currents, thus contributing to fast synaptic signaling in the brain.

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.

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.

Mechanisms shaping fast excitatory postsynaptic currents in the central nervous system

How different factors contribute to determine the time course of the basic element of fast glutamate-mediated excitatory postsynaptic currents (mEPSCs) in the central nervous system has been a focus of interest of neurobiologists for some years. In spite of intensive investigations, these mechanisms are not well understood. In this review, basic hypotheses are summarized, and a new hypothesis is proposed, which holds that desensitization of AMPA receptors plays a major role in shaping the time course of fast mEPSCs. According to the new hypothesis, desensitization shortens the time course of mEPSCs largely by reducing the buffering of glutamate molecules by AMPA receptors. The hypothesis accounts for numerous findings on fast mEPSCs and is expected to be equally fruitful as a framework for further experimental and theoretical investigations.

The role of perisynaptic glial sheaths in glutamate spillover and extracellular Ca(2+) depletion

Recent findings suggest that rapid activation of extrasynaptic receptors and transient depletion of extracellular Ca(2+) may represent an important component of glutamatergic synaptic transmission. These phenomena imply a previously unrecognized role for synaptic glial sheaths: to retard extracellular diffusion in the synaptic vicinity. The present study is an attempt to assess the extent and physiological implications of this retardation using a detailed compartmental model of the typical synaptic environment. The model allows reconstruction of a partial (asymmetric) glial sheath covered with transporter molecules, which gives a more realistic representation of the vicinity of central synapses. Simulations show to what extent, in conditions compatible with physiology, the occupancy of synaptic receptors and the depletion of Ca(2+) in the cleft increase with increased glial coverage. The impact of glial sheaths on synaptic transmission is shown to become greater with smaller synapses and with slower kinetics of perisynaptic ion transients. At a calyceal synapse, a profound temporal filtering of fast Ca(2+) influx is found, and similar phenomena are predicted to occur following simultaneous activation of multiple synapses in the neuropil. The results provide a quantitative guidance for interpretation of physiological experiments that address fast transients of neurotransmitters and small ions in the brain tissue.

Glutamate uptake

Brain tissue has a remarkable ability to accumulate glutamate. This ability is due to glutamate transporter proteins present in the plasma membranes of both glial cells and neurons. The transporter proteins represent the only (significant) mechanism for removal of glutamate from the extracellular fluid and their importance for the long-term maintenance of low and non-toxic concentrations of glutamate is now well documented. In addition to this simple, but essential glutamate removal role, the glutamate transporters appear to have more sophisticated functions in the modulation of neurotransmission. They may modify the time course of synaptic events, the extent and pattern of activation and desensitization of receptors outside the synaptic cleft and at neighboring synapses (intersynaptic cross-talk). Further, the glutamate transporters provide glutamate for synthesis of e.g. GABA, glutathione and protein, and for energy production. They also play roles in peripheral organs and tissues (e.g. bone, heart, intestine, kidneys, pancreas and placenta). Glutamate uptake appears to be modulated on virtually all possible levels, i.e. DNA transcription, mRNA splicing and degradation, protein synthesis and targeting, and actual amino acid transport activity and associated ion channel activities. A variety of soluble compounds (e.g. glutamate, cytokines and growth factors) influence glutamate transporter expression and activities. Neither the normal functioning of glutamatergic synapses nor the pathogenesis of major neurological diseases (e.g. cerebral ischemia, hypoglycemia, amyotrophic lateral sclerosis, Alzheimer's disease, traumatic brain injury, epilepsy and schizophrenia) as well as non-neurological diseases (e.g. osteoporosis) can be properly understood unless more is learned about these transporter proteins. Like glutamate itself, glutamate transporters are somehow involved in almost all aspects of normal and abnormal brain activity.

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.

Quantal release of serotonin

We have studied the origin of quantal variability for small synaptic vesicles (SSVs) and large dense-cored vesicles (LDCVs). As a model, we used serotonergic Retzius neurons of leech that allow for combined amperometrical and morphological analyses of quantal transmitter release. We find that the transmitter amount released by a SSV varies proportionally to the volume of the vesicle, suggesting that serotonin is stored at a constant intravesicular concentration and is completely discharged during exocytosis. Transmitter discharge from LDCVs shows a higher degree of variability than is expected from their size distribution, and bulk release from LDCVs is slower than release from SSVs. On average, differences in the transmitter amount released from SSVs and LDCVs are proportional to the size differences of the organelles, suggesting that transmitter is stored at similar concentrations in SSVs and LDCVs.

Effect of voltage drop within the synaptic cleft on the current and voltage generated at a single synapse

In a model of a single synapse with a circular contact zone and a single concentric zone containing receptor-gated channels, we studied the dependence of the synaptic current on the synaptic cleft width and on the relative size of the receptor zone. During synaptic excitation, the extracellular current entered the cleft and flowed into the postsynaptic cell through receptor channels distributed homogeneously over the receptor zone. The membrane potential and channel currents were smaller toward the cleft center if compared to the cleft edges. This radial gradient was due to the voltage drop produced by the synaptic current on the cleft resistance. The total synaptic current conducted by the same number of open channels was sensitive to changes in the receptor zone radius and the cleft width. We conclude that synaptic geometry may affect synaptic currents by defining the volume resistor of the cleft. The in-series connection of the resistances of the intracleft medium and the receptor channels plays the role of the synaptic voltage divider. This voltage dividing effect should be taken into account when the conductance of single channels or synaptic contacts is estimated from experimental measurements of voltage-current relationships.

Altering the concentration of GABA in the synaptic cleft potentiates miniature IPSCs in rat occipital cortex

We have tested the effect of dextran (40 kDa, 5%) on miniature IPSCs (mIPSCs) recorded in layer V cortical pyramidal cells. This compound increases the amplitude of mIPSCs at room and physiological temperatures by 15%, leaving their duration unaffected at room temperature and slightly increased at physiological temperature. The amplitude increase is attributable to an increase in the number of receptors bound by GABA during synaptic transmission, as shown by the occlusion between the effects of dextran and zolpidem on mIPSC amplitude at room temperature. As dextran presumably enhances the concentration and dwell time of GABA in the synaptic cleft, these results demonstrate that the postsynaptic GABAA receptors are not saturated at room and physiological temperatures.

Electrodiffusion of synaptic receptors: a mechanism to modify synaptic efficacy?

We analysed physical forces that act on synaptic receptor-channels following the release of neurotransmitter. These forces are: 1) electrostatic interaction between receptors, 2) stochastic Brownian diffusion in the membrane, 3) transient electric field force generated by currents through open channels, 4) viscous drag force elicited by the flowing molecules and 5) strong in-membrane friction. By considering alpha-amino-3-hydroxy-5-methyl-4-isoxazolepropionic acid (AMPA) type receptors, we show that, depending on the size and electrophoretic charge of the extracellular receptor domain, release of an excitatory neurotransmitter (glutamate) can induce receptor clustering towards the release site on a fast time scale (8-100 ms). This clustering progresses whenever repetitive synaptic activation exceeds a critical frequency (20-100 s(-1), depending on the currents through individual channels). As a result, a higher proportion of the receptors is exposed to higher glutamate levels. This should increase by 50-100% the peak synaptic current induced by the same amount of released neurotransmitter. In order for this mechanism to contribute to long-term changes of synaptic efficacy, we consider the possibility that the in-membrane motility of the AMPA receptors is transiently increased during synaptic activity, e. g., through the breakage of receptor anchors in the postsynaptic membrane due to activation of N-methyl-d-aspartic acid receptors.

Clustering of the group III metabotropic glutamate receptor 4a at parallel fiber synaptic terminals in the rat cerebellar molecular layer

We report in this study with a pre-embedding immunogold method, the clustering of the group III metabotropic glutamate receptor 4a (mGluR4a) along the presynaptic membrane of parallel fiber synaptic terminals in the cerebellar molecular layer. The mGluR4a clusters were homogeneously distributed and interspaced by about 60 nm. These results suggest a particular arrangement of mGluR4a which might help to a rapid and effective activation of this receptor by glutamate.

Building a bilayer model of the neuromuscular synapse

Progress over the past 10 years has made it possible to construct a simple model of neurotransmitter release. Currently, some models use artificially formed vesicles to represent synaptic vesicles and a planar lipid bilayer as a presynaptic membrane. Fusion of vesicles with the bilayer is via channel proteins in the vesicle membrane and an osmotic gradient. In this paper; a framework is presented for the successful construction of a more complete model of synaptic transmission. This model includes real synaptic vesicles that fuse with a planar bilayer. The bilayer contains acetylcholine receptor (AChR) channels which function as autoreceptors in the membrane. Vesicle fusion is initiated following a Ca2+ flux through voltage-gated Ca2+ channels. Key steps in the plan are validated by mathematical modeling. Specifically, the probability that a reconstituted AChR channel opens following the release of ACh from a fusing vesicle, is calculated as a function of time, quantal content, and number of reconstituted AChRs. Experimentally obtainable parameters for construction of a working synapse are given. The inevitable construction of a full working model will mean that the minimal structures necessary for synaptic transmission are identified. This will open the door in determining regulatory and modulatory factors of transmitter release.