Event Abstract Back to Event Vesicle release statistics at graded potential synapses: modeling study Mantas Gabrielaitis1, 2, 3*, Nikolai M. Chapochnikov1, 2, 3, Tobias Moser2, 3 and Fred Wolf1, 3 1 Max Planck INstitute for Dynamics and Self-Organization, Nonlinear Dynamics, Germany 2 InnerEarLab, Department of Otolaryngology, University of Goettingen Medical School, Germany 3 Bernstein Center for Computational Neuroscience, University of Goettingen, Germany Signal transfer characteristics of chemical synapses are strongly shaped by the molecular organization of the presynaptic active zone. Graded potential synapses are ubiquitous among sensory receptor neurons such as auditory and vestibular hair cells, retinal photoreceptors, and electroreceptors [1]. Generally, these synapses feature high fidelity, precision, and broad dynamic range when encoding transient or prolonged stimuli [1]. Experimental data describing the fine molecular organization of these synapses is rapidly accumulating, e.g., [2,3]. Here, we study an idealized biophysical model of the active zone of a graded potential synapse to relate its topography and molecular kinetic properties to the vesicle release dynamics. The model includes N independent Ca2+ channels described by two state Markov models. These channels can be arbitrarily distributed within the active zone along with M arbitrarily distributed vesicle docking sites. Each docking site is modeled as a two state Markov chain, one state being empty and the other - ready to release a docked vesicle. Release occurs immediately after binding of a single Ca2+ ion to the vesicle's Ca2+ sensor. The rate of binding is proportional to the calcium concentration [Ca2+] at the release site. Because Ca2+ diffusion from a channel to the release sensor is fast compared to the channel open times, we approximate [Ca2+] changes as instantaneous upon Ca2+ channel opening and closing [4]. We describe the [Ca2+] profile by a linearized buffered Ca2+ diffusion approximation [4]. Replenishment of empty docking sites occurs at constant rate. In order to characterize transient and sustained responses of this model active zone, we analyze first release event latency after a voltage step and stationary release inter-event interval distributions. We obtain analytical expressions for these quantities for various active zone topographies. By using our formalism, we examine the influence of Ca2+ channel number, positioning, and gating kinetics, vesicle docking site number and positioning, kinetics of the Ca2+ release sensor, and vesicle replenishment rate. We confirm the analytical results by direct numerical simulations. We find that, in a wide range of situations, the mean and variance of [Ca2+], and population variance of how much each channel contributes to the total [Ca2+] on average largely determine the release dynamics for a given active zone topography. For all active zone topographies, increasing number of release slots shifts the release inter-event interval statistics towards Poissonian. Increasing vesicle docking site and accompanying Ca2+ channel numbers proportionally changes the first event latency in the same way as effectively scaling the channel number. We hope that our insights on this analytically treatable model can help in analyzing how the active zone topography shapes information encoding at graded potential synapses in first approximation.