Code for "Gap junctions desynchronize a neural circuit to stabilize insect flight"

This code produces the figures F3B, F3C, F3D, F3E, F3H, and F3J in the manuscript "Gap junctions desynchronize a neural circuit to stabilize insect flight". The published article can now be found at DOI: 10.1038/s41586-023-06099-0

Overview:

A small circuit is simulated, which consists of five motoneurons that innervate the dorsal longitudinal wing depressor muscle of Drosophila melanogaster. Each motoneuron is described by a single-compartment model based on Berger et al. (2015). Five identical single-neuron models are coupled by linear non-rectifying gap junction currents. Then, channel conductances, input currents, and coupling strengths are varied. The resulting outputs are compared to experimentally observed firing patterns. Please see the main manuscript for details. 

Additional information:

Instructions for installation and running the code can be found in the README. The data folder contains experimentally recorded spike trains for the comparison between models and experiments in Figures F3D and F3H. Please note that the data for the other figures have been uploaded separately and can be found at DOI: 10.5281/zenodo.7737730