# (Chartove et al., 2019) Striatal FSI and SPN network DynaSim full simulation reproduction files.
This contains the complete [DynaSim](https://github.com/DynaSim/DynaSim) code files (in `dynasim`),
model mechanism files (in `dynasim/models/striatal_mechanisms`), and simulation runscripts
(in `runfiles`) needed for simulation of the striatum of:
A biophysical model of striatal microcircuits suggests delta/theta-rhythmically interleaved gamma
and beta oscillations mediate periodicity in motor control. Julia A. K. Chartove, Michelle M.
McCarthy, Benjamin R. Pittman-Polletta, Nancy J. Kopell. In preparation.
Adding the `dynasim` directory and all its subdirectories to your MATLAB path should enable you
to run the simulations in `runfiles`. Several plotting functions used in Dynasim require the MATLAB
Signal Processing Toolbox, so installing this toolbox is recommended if you want to use the built
in data visualizations.
The runfiles are labeled as follows:
FSI_1.m: This file simulates a single model fast-spiking interneuron at varying input levels, and
should be straightforward to run on a personal computer.
FSI_250_dynasim.m: This file simulates four model networks of 250 cells each:
sim1: Baseline dopamiergic tone, isolated FSI and SPN subnetworks (Figure 4i,5i)
sim2: High dopamiergic tone, isolated FSI and SPN subnetworks (Figure 4ii,5ii)
sim3: Baseline dopaminergic tone, full connectivity (Figure 4i,6i,7i)
sim4: High dopaminergic tone, full connectivity (Figure 4ii,6ii,7ii)
Each simulation runs for 4 simulated seconds at a high fixed-step resolution (Runge-Kutta, 0.01 ms)
before being downsampled to a time resolution of 0.1 ms. These simulations are set to run with a
default memory/RAM allowance of 64 GB and produce output files 300-400 MB in size. Running them may
require a parallel or high-performance scientific computing setup.
The files produced by golomb_100_dynasim.m can be used to replicate figures from the paper. In
`replicate_figures`, the files `make_fig4.m`, `make_fig5and6.m`, and `make_fig7.m` can be used to
generate the respective figures. After running FSI_250_dynasim.m, navigate to `full_network/data`
and run the following commands:
`make_fig1`,`make_fig2`,`make_fig3`, `make_fig_s1`, and `make_fig_s2` require summary data generated
over many runs and cannot be run based on the files in `runfiles`; instead, the data used to generate
them has been included in `replicate_figures/example_data`. If you are interested in running the
simulations used to produce these figures, they are included in `statistics`; however, the number of
simulations necessary to run is large enough to likely require a parallel computing cluster, and the
data generated requires some degree of cleaning by hand for use in figures, so this is not recommended.
The `statistics` folder is basically included for transparency's sake; you shouldn't need to use it
to replicate the main findings of the paper. If you would like a detailed explanation of how these
data were used to prepare figures, contact me at chartove 'at' bu 'dot' edu.
FSI models used in this simulation were based on
Golomb, D., Donner, K., Shacham, L., Shlosberg, D., Amitai, Y., & Hansel, D. (2007). Mechanisms of
firing patterns in fast-spiking cortical interneurons. PLoS Computational Biology, 3(8), 1498–1512.
and SPN models were based on
McCarthy, M. M., Moore-Kochlacs, C., Gu, X., Boyden, E. S., Han, X., & Kopell, N. (2011). Striatal
origin of the pathologic beta oscillations in Parkinson’s disease. Proceedings of the National
Academy of Sciences of the United States of America, 108(28), 11620–11625.
This code diverges from the given equations of (Golomb et al., 2007) and (McCarthy et al., 2011) in
several respects, discussed in the text of (Chartove et al., 2019)
Also note that figures in the text have been reformatted by hand in Inkscape; while no data was
changed, plots will not appear identical to those in the text.
The included copy of DynaSim is using commit
Later releases of Dynasim may not be compatible and may result in errors if used. Simulations were
run on MATLAB version 2017b.