The model response was recorded either in voltage clamp or in current-clamp at the soma, and the voltage-dependent ionic conductances were modified in some cases to test their impact on BPR generation. The membrane voltage, the calcium concentration, the ionic currents from each model section were recorded and saved in a nested MATLAB file. The model response properties were analyzed using MATLAB routines (MathWorks, Natick, MA, USA) and custom made Python scripts, using data recorded during each simulations. Each simulation data file contained information about the membrane voltage, ionic channels and synaptic currents recorded in specific location throughout the models.