Utilising activity patterns of a complex biophysical network model to optimise intra-striatal deep brain stimulation

A large-scale biophysical network model for the isolated striatal body is developed to optimise potential intrastriatal deep brain stimulation applied to, e.g. obsessive-compulsive disorder. The model is based on modified Hodgkin-Huxley equations with small-world connectivity, while the spatial information about the positions of the neurons is taken from a detailed human atlas. The model produces neuronal spatiotemporal activity patterns segregating healthy from pathological conditions. Three biomarkers were used for the optimisation of stimulation protocols regarding stimulation frequency, amplitude and localisation: the mean activity of the entire network, the frequency spectrum of the entire network (rhythmicity) and a combination of the above two. By minimising the deviation of the aforementioned biomarkers from the normal state, we compute the optimal deep brain stimulation parameters, regarding position, amplitude and frequency. Our results suggest that in the DBS optimisation process, there is a clear trade-off between frequency synchronisation and overall network activity, which has also been observed during in vivo studies.