Metamodelling of a two-population spiking neural network

Abstract In computational neuroscience, hypotheses are often formulated as bottom-up mechanistic models of the systems in question, consisting of differential equations that can be numerically integrated forward in time. Candidate models can then be validated by comparison against experimental data. The model outputs of neural network models depend on both neuron parameters, connectivity parameters and other model inputs. Successful model fitting requires sufficient exploration of the model parameter space, which can be computationally demanding. Additionally, identifying degeneracy in the parameters, i.e. different combinations of parameter values that produce similar outputs, is of interest, as they define the subset of models consistent with the data. In this computational study, we apply metamodels to a two-population recurrent spiking network of point-neurons, the so-called Brunel network. Metamodels are data-driven approximations to more complex models with more desirable computational properties, which can be run considerably faster than the original model. Specifically, we apply two different metamodelling techniques, Masked Autoregressive Flows and Deep Gaussian Process Regression, to estimate the power spectra of two different signals; the population spiking activities and the local field potential. We find that the metamodels are able to accurately model the power spectra in the asynchronous irregular regime. Using the metamodels, we invert the network, and estimate the probability distributions over parameters that can give rise to observed simulation outputs. We find that these distributions correctly identify parameter combinations that give similar model outputs, and that the parameters are constrained to different degrees by observing the two different signals. Author summary In computational neuroscience, mechanistic models are used to simulate networks of neurons. These models exhibit complex dynamics, and the parameters of the neurons and connections between neurons shape the model’s behaviour. Due to the model complexity, running the simulations and fitting the model to experimental data can be computationally demanding. In this study, we apply metamodelling techniques, data-driven approximations that are much faster to run, to two different signals generated by a two-population recurrent network model, the population spiking activities and the local field potential (LFP). Further, we invert the metamodels, and demonstrate that it can reliably find the different combinations of parameters that can give rise to an observed simulation output. We compare the metamodels on both the forward and inverse problem, and investigate to what degree the parameters are constrained by the two different signals.