Systems-level modeling of neuronal circuits for leech swimming

This paper describes a mathematical model of the neuronal central pattern generator (CPG) that controls the rhythmic body motion of the swimming leech. The systems approach is employed to capture the neuronal dynamics essential for generating coordinated oscillations of cell membrane potentials by a simple CPG architecture with a minimal number of parameters. Based on input/output data from physiological experiments, dynamical components (neurons and synaptic interactions) are first modeled individually and then integrated into a chain of nonlinear oscillators to form a CPG. We show through numerical simulations that the values of a few parameters can be estimated within physiologically reasonable ranges to achieve good fit of the data with respect to the phase, amplitude, and period. This parameter estimation leads to predictions regarding the synaptic coupling strength and intrinsic period gradient along the nerve cord, the latter of which agrees qualitatively with experimental observations.