Modeling of quantal neurotransmitter release kinetics in the presence of fixed and mobile calcium buffers

The Mechanism of α2 adrenoreceptor-dependent Modulation of Neurotransmitter Release at the Neuromuscular Junctions

α2-Adrenoreceptors (ARs) are main Gi-protein coupled autoreceptors in sympathetic nerve terminals and targets for dexmedetomidine (DEX), a widely used sedative. We hypothesize that α2-ARs are also potent regulators of neuromuscular transmission via G protein-gated inwardly rectifying potassium (GIRK) channels. Using extracellular microelectrode recording of postsynaptic potentials, we found DEX-induced inhibition of spontaneous and evoked neurotransmitter release as well as desynchronization of evoked exocytotic events in the mouse diaphragm neuromuscular junction. These effects were suppressed by SKF-86,466, a selective α2-AR antagonist. An activator of GIRK channels ML297 had the same effects on neurotransmitter release as DEX. By contrast, inhibition of GIRK channels with tertiapin-Q prevented the action of DEX on evoked neurotransmitter release, but not on spontaneous exocytosis. The synaptic vesicle exocytosis is strongly dependent on Ca2+ influx through voltage-gated Ca2+ channels (VGCCs), which can be negatively regulated via α2-AR - GIRK channel axis. Indeed, inhibition of P/Q-, L-, N- or R-type VGCCs prevented the inhibitory action of DEX on evoked neurotransmitter release; antagonists of P/Q- and N-type channels also suppressed the DEX-mediated desynchronization of evoked exocytotic events. Furthermore, inhibition of P/Q-, L- or N-type VGCCs precluded the frequency decrease of spontaneous exocytosis upon DEX application. Thus, α2-ARs acting via GIRK channels and VGCCs (mainly, P/Q- and N-types) exert inhibitory effect on the neuromuscular communication by attenuating and desynchronizing evoked exocytosis. In addition, α2-ARs can suppress spontaneous exocytosis through GIRK channel-independent, but VGCC-dependent pathway.

Simulations of active zone structure and function at mammalian NMJs predict that loss of calcium channels alone is not sufficient to replicate LEMS effects

Lambert-Eaton myasthenic syndrome (LEMS) is an autoimmune-mediated neuromuscular disease thought to be caused by autoantibodies against P/Q-type voltage-gated calcium channels (VGCCs), which attack and reduce the number of VGCCs within transmitter release sites (active zones; AZs) at the neuromuscular junction (NMJ), resulting in neuromuscular weakness. However, patients with LEMS also have antibodies to other neuronal proteins, and about 15% of patients with LEMS are seronegative for antibodies against VGCCs. We hypothesized that a reduction in the number of P/Q-type VGCCs alone is not sufficient to explain LEMS effects on transmitter release. Here, we used a computational model to study a variety of LEMS-mediated effects on AZ organization and transmitter release constrained by electron microscopic, pharmacological, immunohistochemical, voltage imaging, and electrophysiological observations. We show that models of healthy AZs can be modified to predict the transmitter release and short-term facilitation characteristics of LEMS and that in addition to a decrease in the number of AZ VGCCs, disruption in the organization of AZ proteins, a reduction in AZ number, a reduction in the amount of synaptotagmin, and the compensatory expression of L-type channels outside the remaining AZs are important contributors to LEMS-mediated effects on transmitter release. Furthermore, our models predict that antibody-mediated removal of synaptotagmin in combination with disruption in AZ organization alone could mimic LEMS effects without the removal of VGCCs (a seronegative model). Overall, our results suggest that LEMS pathophysiology may be caused by a collection of pathological alterations to AZs at the NMJ, rather than by a simple loss of VGCCs.NEW & NOTEWORTHY We used a computational model of the active zone (AZ) in the mammalian neuromuscular junction to investigate Lambert-Eaton myasthenic syndrome (LEMS) pathophysiology. This model suggests that disruptions in presynaptic active zone organization and protein content (particularly synaptotagmin), beyond the simple removal of presynaptic calcium channels, play an important role in LEMS pathophysiology.

Presynaptic Acetylcholine Receptors Modulate the Time Course of Action Potential-Evoked Acetylcholine Quanta Secretion at Neuromuscular Junctions

For effective transmission of excitation in neuromuscular junctions, the postsynaptic response amplitude must exceed a critical level of depolarization to trigger action potential spreading along the muscle-fiber membrane. At the presynaptic level, the end-plate potential amplitude depends not only on the acetylcholine quanta number released from the nerve terminals in response to the nerve impulse but also on a degree of synchronicity of quanta releases. The time course of stimulus-phasic synchronous quanta secretion is modulated by many extra- and intracellular factors. One of the pathways to regulate the neurosecretion kinetics of acetylcholine quanta is an activation of presynaptic autoreceptors. This review discusses the contribution of acetylcholine presynaptic receptors to the control of the kinetics of evoked acetylcholine release from nerve terminals at the neuromuscular junctions. The timing characteristics of neurotransmitter release is nowadays considered an essential factor determining the plasticity and efficacy of synaptic transmission.

Bayesian analysis of the kinetics of quantal transmitter secretion at the neuromuscular junction

The timing of transmitter release from nerve endings is considered nowadays as one of the factors determining the plasticity and efficacy of synaptic transmission. In the neuromuscular junction, the moments of release of individual acetylcholine quanta are related to the synaptic delays of uniquantal endplate currents recorded under conditions of lowered extracellular calcium. Using Bayesian modelling, we performed a statistical analysis of synaptic delays in mouse neuromuscular junction with different patterns of rhythmic nerve stimulation and when the entry of calcium ions into the nerve terminal was modified. We have obtained a statistical model of the release timing which is represented as the summation of two independent statistical distributions. The first of these is the exponentially modified Gaussian distribution. The mixture of normal and exponential components in this distribution can be interpreted as a two-stage mechanism of early and late periods of phasic synchronous secretion. The parameters of this distribution depend on both the stimulation frequency of the motor nerve and the calcium ions' entry conditions. The second distribution was modelled as quasi-uniform, with parameters independent of nerve stimulation frequency and calcium entry. Two different probability density functions for the distribution of synaptic delays suggest at least two independent processes controlling the time course of secretion, one of them potentially involving two stages. The relative contribution of these processes to the total number of mediator quanta released depends differently on the motor nerve stimulation pattern and on calcium ion entry into nerve endings.

Estimation of presynaptic calcium currents and endogenous calcium buffers at the frog neuromuscular junction with two different calcium fluorescent dyes

At the frog neuromuscular junction, under physiological conditions, the direct measurement of calcium currents and of the concentration of intracellular calcium buffers—which determine the kinetics of calcium concentration and neurotransmitter release from the nerve terminal—has hitherto been technically impossible. With the aim of quantifying both Ca²⁺ currents and the intracellular calcium buffers, we measured fluorescence signals from nerve terminals loaded with the low-affinity calcium dye Magnesium Green or the high-affinity dye Oregon Green BAPTA-1, simultaneously with microelectrode recordings of nerve-action potentials and end-plate currents. The action-potential-induced fluorescence signals in the nerve terminals developed much more slowly than the postsynaptic response. To clarify the reasons for this observation and to define a spatiotemporal profile of intracellular calcium and of the concentration of mobile and fixed calcium buffers, mathematical modeling was employed. The best approximations of the experimental calcium transients for both calcium dyes were obtained when the calcium current had an amplitude of 1.6 ± 0.08 pA and a half-decay time of 1.2 ± 0.06 ms, and when the concentrations of mobile and fixed calcium buffers were 250 ± 13 μM and 8 ± 0.4 mM, respectively. High concentrations of endogenous buffers define the time course of calcium transients after an action potential in the axoplasm, and may modify synaptic plasticity.

Nanodomain coupling explains Ca²⁺ independence of transmitter release time course at a fast central synapse

A puzzling property of synaptic transmission, originally established at the neuromuscular junction, is that the time course of transmitter release is independent of the extracellular Ca²⁺ concentration ([Ca²⁺]_o), whereas the rate of release is highly [Ca²⁺]_o-dependent. Here, we examine the time course of release at inhibitory basket cell-Purkinje cell synapses and show that it is independent of [Ca²⁺]_o. Modeling of Ca²⁺-dependent transmitter release suggests that the invariant time course of release critically depends on tight coupling between Ca²⁺ channels and release sensors. Experiments with exogenous Ca²⁺ chelators reveal that channel-sensor coupling at basket cell-Purkinje cell synapses is very tight, with a mean distance of 10–20 nm. Thus, tight channel-sensor coupling provides a mechanistic explanation for the apparent [Ca²⁺]_o independence of the time course of release.

DOI:http://dx.doi.org/10.7554/eLife.04057.001

The nervous system sends information around the body in the form of electrical signals that travel through cells called neurons. However, these electrical signals cannot cross the synapses between neurons. Instead, the information is carried across the synapse by molecules called neurotransmitters.

Calcium ions control the release of neurotransmitters. There is a high concentration of calcium ions outside the neuron but they are not able to pass through the cell membrane under normal conditions. However, when an electrical impulse reaches the synapse, ion channels in the membrane open and allow calcium ions to enter the cell. Once inside, the ions activate the release of neurotransmitters by binding to proteins called release sensors.

Several experiments on the release of neurotransmitters have studied the synapses between neurons and muscle fibers. These studies found that the higher the concentration of ions outside the neuron, the higher the rate at which the neurotransmitters were released. However, the timing of release—the length of time over which the neurotransmitters were released—did not depend on the concentration of calcium ions.

Arai and Jonas have now studied neurotransmitter release at a synapse in a region of the brain called the cerebellum. These experiments also found that the timing of the release did not depend on the ion concentration, suggesting that this may be a general property of neurotransmitter release.

To find out more, Arai and Jonas created a mathematical model of neurotransmitter release. This model suggests that for the timing of release to remain the same, the ion channel and the release sensor must be located close together in the presynaptic terminal. If they are not close together, the timing of release becomes blurred and more dependent on the external calcium concentration.

Further experiments confirm the prediction of the model by showing that the calcium channels and the release sensors in these synapses are very close together. The next challenge will be to find out whether the conclusions are also valid for other synapses where the calcium channels and release sensors are further apart.

DOI:http://dx.doi.org/10.7554/eLife.04057.002

An excess-calcium-binding-site model predicts neurotransmitter release at the neuromuscular junction

Despite decades of intense experimental studies, we still lack a detailed understanding of synaptic function. Fortunately, using computational approaches, we can obtain important new insights into the inner workings of these important neural systems. Here, we report the development of a spatially realistic computational model of an entire frog active zone in which we constrained model parameters with experimental data, and then used Monte Carlo simulation methods to predict the Ca(2+)-binding stoichiometry and dynamics that underlie neurotransmitter release. Our model reveals that 20-40 independent Ca(2+)-binding sites on synaptic vesicles, only a fraction of which need to bind Ca(2+) to trigger fusion, are sufficient to predict physiological release. Our excess-Ca(2+)-binding-site model has many functional advantages, agrees with recent data on synaptotagmin copy number, and is the first (to our knowledge) to link detailed physiological observations with the molecular machinery of Ca(2+)-triggered exocytosis. In addition, our model provides detailed microscopic insight into the underlying Ca(2+) dynamics during synapse activation.

Kinetics of neurotransmitter release in neuromuscular synapses of newborn and adult rats

The kinetics of the phasic synchronous and delayed asynchronous release of acetylcholine quanta was studied at the neuromuscular junctions of aging rats from infant to mature animals at various frequencies of rhythmic stimulation of the motor nerve. We found that in infants 6 (P6) and 10 (P10) days after birth a strongly asynchronous phase of quantal release was observed, along with a reduced number of quanta compared to the synapses of adults. The rise time and decay of uni-quantal end-plate currents were significantly longer in infant synapses. The presynaptic immunostaining revealed that the area of the synapses in infants was significantly (up to six times) smaller than in mature junctions. The intensity of delayed asynchronous release in infants increased with the frequency of stimulation more than in adults. A blockade of the ryanodine receptors, which can contribute to the formation of delayed asynchronous release, had no effect on the kinetics of delayed secretion in the infants unlike synapses of adults. Therefore, high degree of asynchrony of quantal release in infants is not associated with the activity of ryanodine receptors and with the liberation of calcium ions from intracellular calcium stores.

Redox-sensitive synchronizing action of adenosine on transmitter release at the neuromuscular junction

The kinetics of neurotransmitter release was recognized recently as an important contributor to synaptic efficiency. Since adenosine is the ubiquitous modulator of presynaptic release in peripheral and central synapses, in the current project we studied the action of this purine on the timing of acetylcholine quantal release from motor nerve terminals in the skeletal muscle. Using extracellular recording from frog neuromuscular junction we tested the action of adenosine on the latencies of single quantal events in the pro-oxidant and antioxidant conditions. We found that adenosine, in addition to previously known inhibitory action on release probability, also synchronized release by removing quantal events with long latencies. This action of adenosine on release timing was abolished by oxidants whereas in the presence of the antioxidant the synchronizing action of adenosine was further enhanced. Interestingly, unlike the timing of release, the inhibitory action of adenosine on release probability was redox-independent. Modulation of release timing by adenosine was mediated by purinergic A1 receptors as it was eliminated by the specific A1 antagonist 8-cyclopentyl-1,3-dipropylxanthine (DPCPX) and mimicked by the specific A1 agonist N(6)-cyclopentyl-adenosine. Consistent with data obtained from dispersion of single quantal events, adenosine also reduced the rise-time of multiquantal synaptic currents. The latter effect was reproduced in the model based on synchronizing effect of adenosine on release timing. Thus, adenosine which is generated at the neuromuscular junction from the breakdown of the co-transmitter ATP induces the synchronization of quantal events. The effect of adenosine on release timing should preserve the fidelity of synaptic transmission via "cost-effective" use of less transmitter quanta. Our findings also revealed important crosstalk between purinergic and redox modulation of synaptic processes which could take place in the elderly or in neuromuscular diseases associated with oxidative stress like lateral amyotrophic sclerosis.

Combined computational and experimental approaches to understanding the Ca(2+) regulatory network in neurons

Ca(2+) is a ubiquitous signaling ion that regulates a variety of neuronal functions by binding to and altering the state of effector proteins. Spatial relationships and temporal dynamics of Ca(2+) elevations determine many cellular responses of neurons to chemical and electrical stimulation. There is a wealth of information regarding the properties and distribution of Ca(2+) channels, pumps, exchangers, and buffers that participate in Ca(2+) regulation. At the same time, new imaging techniques permit characterization of evoked Ca(2+) signals with increasing spatial and temporal resolution. However, understanding the mechanistic link between functional properties of Ca(2+) handling proteins and the stimulus-evoked Ca(2+) signals they orchestrate requires consideration of the way Ca(2+) handling mechanisms operate together as a system in native cells. A wide array of biophysical modeling approaches is available for studying this problem and can be used in a variety of ways. Models can be useful to explain the behavior of complex systems, to evaluate the role of individual Ca(2+) handling mechanisms, to extract valuable parameters, and to generate predictions that can be validated experimentally. In this review, we discuss recent advances in understanding the underlying mechanisms of Ca(2+) signaling in neurons via mathematical modeling. We emphasize the value of developing realistic models based on experimentally validated descriptions of Ca(2+) transport and buffering that can be tested and refined through new experiments to develop increasingly accurate biophysical descriptions of Ca(2+) signaling in neurons.

Ca2+ dependence of the binomial parameters p and n at the mouse neuromuscular junction

The Ca(2+) dependence of synaptic quantal release is generally thought to be restricted to probability of vesicular release. However, some studies have suggested that the number of release sites (n) at the neuromuscular junction (NMJ) is also Ca(2+) dependent. In this study, we recorded endplate currents over a wide range of extracellular Ca(2+) concentrations and found the expected Ca(2+) dependency of release. A graphical technique was used to estimate p (probability of release) and n using standard binomial assumptions. The results suggested n was Ca(2+) dependent. The data were simulated using compound binomial statistics with variable n (Ca(2+) dependent) or fixed n (Ca(2+) independent). With fixed n, successful simulation of increasing Ca(2+) required that p increase abruptly at some sites from very low to high values. Successful simulation with variable n required the introduction of previously silent release sites (p = 0) with high values of p. Thus the success of both simulations required abrupt, large increases of p at a subset of release sites with initially low or zero p. Estimates of the time course of release obtained by deconvolving evoked endplate currents with average miniature endplate currents decreased slightly as Ca(2+) increased, thus arguing against sequential release of multiple quanta at higher Ca(2+) levels. Our results suggest that the apparent Ca(2+) dependence of n at the NMJ can be explained by an underlying Ca(2+) dependence of a spatially variable p such that p increases abruptly at a subset of sites as Ca(2+) is increased.