Due to the heterogeneous macromolecular crowding and geometrical irregularity at central excitatory synapses, the diffusion coefficient of glutamate may exhibit spatial variation across the cleft space. To take into account the effect of emergent cleft heterogeneity on the generation of excitatory postsynaptic currents (EPSCs), a gamma statistical distribution of the glutamate diffusion coefficient is considered and, using the principle of superstatistics, the glutamate transients are computed as well as the activation of AMPA receptors is performed. This model demonstrates the numerical simulation of the Brownian diffusion of glutamate under distributed diffusion coefficient, the subsequent stochastic activation of AMPA receptors using Milstein-Nicoll scheme and modified Gillespie algorithm with minimum time-step correction, and the eventual stochastic profile of EPSC generation. The study is based on the CA1 synapses located at the dendrites of CA1 pyramidal neurons in the mammalian hippocampal region.
Model Type: Synapse
Simulation Environment: MATLAB
Implementer(s): Gupta, Rahul [gupta.sbt at gmail.com]