- Poster presentation
- Open Access
Is the Langevin phase equation an efficient model for stochastic limit cycle oscillators in real neurons?
BMC Neurosciencevolume 10, Article number: P236 (2009)
The Langevin phase equation dϕ/dt = 1 + Z(ϕ)(G(t) + σs(t)), where ϕ is the phase, which is disturbed by a perturbation G and Langevin force s of intensity σ, and Z is the phase response curve (PRC), has been deemed to be a good model for stochastic limit cycle oscillators , and it has been extensively used in theoretical neuroscience as a model neural oscillator . Inspired by the theoretical research, experimental researchers have measured PRCs, but none of them have identified the Langevin phase equation for real neurons directly. In fact, biological experiments have yet to show whether this equation is a good model for neural oscillators.
Here, we demonstrate that the Langevin phase equation is a good model for neural oscillators in rat hippocampal CA1 pyramidal neurons, through two steps in the same neuron. (A) Estimation step: We estimated the parameters of the Langevin phase equation, i.e., a PRC and the intensity of the Langevin force, from physiological noisy data (phase shifts disturbed by one-shot rectangle pulses 20 pA, 5 msec) by using the MAP estimation algorithm, from our previous study . Figure 1 shows the estimated PRC (a) and the hyper parameters (b). The estimated PRC was fully positive, which indicates that hippocampal CA1 pyramidal neurons could be classified as Type I excitability neurons, while the estimated Langevin force intensity suggests that the inherent noise in the recording neuron is 9.0 pA. (B) Prediction step: We injected two different periodic perturbations, consisting successive rectangular pulses. One consisted of 20 pA, 5 msec pulses (the same values used in the PRC estimation), the other 10 pA, 20 msec pulses. We obtained histograms of phase differences between spikes and pulses for each perturbation. Then we ascertained whether the Fokker Planck equation (FPE) derived from the Langevin phase equation with the estimated PRC and Langevin force intensity could predict these histograms.
Results and discussion
The distribution derived from FPE was in good agreement with the experimental histogram (Fig. 2), even when the shape of the perturbation was not the same as the one used in the PRC measurement (Fig. 2b). This result suggests that the Langevin phase equation can describe the stochastic dynamics of neural oscillators, no matter what perturbation disturbs them.
Nakao H, Teramae J, Ermentrout GB: Comment on "phase reduction of stochastic limit cycle oscillators". [http://arxiv.org/abs/0812.3205]
Ermentrout GB, Galan RF, Urban NN: Relating neural dynamics to neural coding. Phys Rev Lett. 2007, 99: 248103-10.1103/PhysRevLett.99.248103.
Ota K, Omori T, Aonishi T: MAP estimation algorithm for phase reponse curves based on analysis of the observation process. J Comput Neurosci.