Multiscale stochastic modeling of the tetrapartite synapse combining finite-volume geometry, Monte Carlo and Markov methods with Brownian diffusion to simulate vesicle release, neurotransmitter spread and postsynaptic channel dynamics calibrated to experimental data.