Action potential-triggered release of neurotransmitters at chemical synapses forms the key basis of communication between two neurons. To quantify the stochastic dynamics of the number of neurotransmitters released, we investigate a model where neurotransmitter-filled vesicles attach to a finite number of docking sites in the axon terminal, and are subsequently released when the action potential arrives. We formulate the model as a Stochastic Hybrid System (SHS) that combines three key noise mechanisms: random arrival of action potentials, stochastic refilling of docking sites, and probabilistic release of docked vesicles. This SHS representation is used to derive exact analytical formulas for the mean and noise (as quantified by Fano factor) in the number of vesicles released per action potential. Interestingly, results show that in relevant parameter regimes, noise in the number of vesicles released is sub-Poissonian at low frequencies, super-Poissonian at intermediate frequencies, and approaches a Poisson limit at high frequencies.In contrast, noise in the number of neurotransmitters in the synaptic cleft is always super-Poissonian, but is lowest at intermediate frequencies.We further investigate changes in these noise properties for non-Poissonian arrival of action potentials, and when the probability of release is frequency dependent. In summary, these results provide the first glimpse into synaptic parameters not only determining the mean synaptic strength, but also shaping its stochastic dynamics that is critical for information transfer between neurons.rich literature regarding synaptic transmission as a deterministic process [1][2][3][4][5], increasing evidence points to diverse noise mechanisms at play during this process [6][7][8][9][10]. Counterintuitively, noise has been shown to sometimes enhance signaling and information transfer between neurons [11][12][13][14][15][16].Stochastic Hybrid Systems (SHS) constitute an important class of mathematical models that integrate discrete stochastic events with continuous dynamics. Given their generality and scope, SHS have been successfully used for modeling stochastic phenomena in a variety of biological processes [17][18][19][20][21][22][23][24][25][26][27][28][29][30][31][32][33], including neuronal dynamics [34,35]. Building up on our previous work [36], we consider the SHS framework for modeling the neurotransmitter dynamics. More specifically, the model consists of M ∈ {1, 2, . . .} docking site at the axon terminal (Fig. 1). Neurotransmitter-filled vesicles attach to these docking sites with a given probabilistic rate that is proportional to M − n, where n is the number of already docked vesicles. Action Potentials (APs) arrive at a given frequency, such that, the time between two successive APs follows an arbitrary positively-valued probability density function g. Each docked vesicle has a certain probability of release, and based on this probability, APs cause a fraction of the docked vesicles to release their neurotransmitter content into the synapt...