A model is constructed to study the phenomenon of bursting in cultured neuronal networks by considering the effects of glial release and the extrasynaptic receptors on neurons. In the frequently observed situations of synchronized bursting, the whole neuronal network can be described by a mean-field model. In this model, the dynamics of the synchronized network in the presence of glia is represented by an effective two-compartment neuron with stimulations on both the dendrite and soma. Numerical simulations of this model show that most of the experimental observations in bursting, in particular the high plateau and the slow repolarization, can be reproduced. Our findings suggest that the effects of glia release and extrasynaptic receptors, which are usually neglected in neuronal models, can become important in intense network activities. Furthermore, simulations of the model are also performed for the case of glia-suppressed cultures to compare with recent experimental results.
