Geometric tortuosity at invaginating rod synapses slows glutamate diffusion and shapes synaptic responses: insights from anatomically realistic Monte Carlo simulations

At the first synapse in the vertebrate retina, rod photoreceptor terminals form deep invaginations occupied by multiple second-order rod bipolar and horizontal cell (RBP and HC) dendrites. Synaptic vesicles are released into this invagination at multiple sites beneath an elongated presynaptic ribbon. We investigated the impact of this complex architecture on the diffusion of synaptic glutamate and activity of postsynaptic receptors. We obtained serial electron micrographs of mouse retina and reconstructed four rod terminals along with their postsynaptic RBP and HC dendrites. We incorporated these structures into an anatomically realistic Monte Carlo simulation of neurotransmitter diffusion and receptor activation. We compared passive diffusion of glutamate in these realistic structures to existing, geometrically simplified models of the synapse and found that glutamate exits anatomically realistic synapses ten times more slowly than previously predicted. By comparing simulations with electrophysiological recordings, we modeled synaptic activation of EAAT5 glutamate transporters in rods, AMPA receptors on HC dendrites, and metabotropic glutamate receptors (mGluR6) on RRBP dendrites. Our simulations suggested that ∼3,000 EAAT5 transporters populate the rod presynaptic membrane and that, while uptake by surrounding glial Müller cells retrieves much of the glutamate released by rods, binding and uptake by EAAT5 influences RBP response kinetics. The relatively long lifetime of glutamate within the cleft allows mGluR6 on RBP dendrites to temporally integrate the steady stream of vesicles released at this synapse in darkness. Glutamate’s tortuous diffusional path through realistic synaptic geometry confers quantal variability, as release from nearby ribbon sites exerts larger effects on RBP and HC receptors than release from more distant sites. While greater integration may allow slower sustained release rates, added quantal variability complicates the challenging task of detecting brief decreases in release produced by rod light responses at scotopic threshold.