Abstract
A three-dimensional model of the reaction-diffusion processes of a neurotransmitter and its ligand receptor in a disk shaped volume is proposed which represents the transmission process of acetylcholine in the synaptic cleft in the neuromuscular junction. The behavior of the reaction-diffusion system is described by a three-dimensional diffusion equation with nonlinear reaction terms due to the rate processes of acetylcholine with the receptor. A new stable and accurate numerical method is used to solve the equations with Neumann boundaries in cylindrical coordinates. The simulation analysis agrees with experimental measurements of end-plate current, and agrees well with the results of the conformational state of the acetylcholine receptor as a function of time and acetylcholine concentration of earlier investigations with a smaller error compared to experiments. Asymmetric emission of acetylcholine in the synaptic cleft and the subsequent effects on open receptor population is simulated. Sensitivity of the open receptor dynamics to the changes in the diffusion parameters and neuromuscular junction volume is investigated. The effects of anisotropic diffusion and non-symmetric emission of transmitter at the presynaptic membrane is simulated.