Abstract:A bstract Connectome-based neural mass modelling is the emerging computational neuroscience paradigm for simulating large-scale network dynamics observed in whole-brain activity measurements such as fMRI, M/EEG, and related techniques. Estimating physiological parameters by fitting these models to empirical data is challenging however, due to large network sizes, often physiologically detailed fast-timescale system equations, and the need for long (e.g. tens of minutes) simulation runs. Here we introduce a novel approach to connectome-based neural mass model parameter estimation by employing optimization tools developed for deep learning. We cast the system of differential equations representing both neural and haemodynamic activity dynamics as a deep neural network, implemented within a widely used machine learning programming environment (PyTorch). This allows us to use robust industry-standard optimization algorithms, automatic differentiation for computation of gradients, and other useful functionality. The approach is demonstrated using a connectome-based network with nodal dynamics specified by the two-state RWW mean-field neural mass model equations, which we use here as a model of fMRI-measured activity and correlation fluctuations. Additional optimization constraints are explored and prove fruitful, including restricting the model to domains of parameter space near a bifurcation point that yield metastable dynamics. Using these techniques, we first show robust recovery of physiological model parameters in synthetic data and then, as a proof-of-principle, apply the framework to modelling of empirical resting-state fMRI data from the Human Connectome Project database. For resting state activity, the system can be understood as a deep net that receives uncorrelated noise on its input layer, which is transformed into network-wide modelled functional connectivity on its output layer. This is consistent with the prevailing conception in theoretical neuroscience of resting-state functional connectivity patterns as an emergent phenomenon that is driven by (effectively) random activity fluctuations, which are then in turn spatiotemporally filtered by anatomical connectivity and local neural dynamics.

Efficient Reconstruction of Neural Mass Dynamics Modeled by Linear-Threshold Networks

A Robust Evolutionary Optimisation Approach For Parameterising A Neural Mass Model

Deep Learning-Based Parameter Estimation for Neurophysiological Models of Neuroimaging Data

Binless Kernel Machine: Modeling Spike Train Transformation for Cognitive Neural Prostheses

Efficient state and parameter estimation for high-dimensional nonlinear system identification with application to MEG brain network modeling

Robust Reconstruction Of Physiological Parameters From Dynamic Pet Data

Bayesian Approaches to Nonlinear Network Reconstruction

A Self-Organizing State-Space-Model Approach for Parameter Estimation in Hodgkin-Huxley-Type Models of Single Neurons

Is the brain macroscopically linear? A system identification of resting state dynamics

Reconstruction of sparse recurrent connectivity and inputs from the nonlinear dynamics of neuronal networks

Reconstruction of Neuronal Input Through Modeling Single-Neuron Dynamics and Computations

Reconstructing dynamics of complex systems from noisy time series with hidden variables

Contraction-Based Identification of a Neuron Model with Nonlinear Parameterization via Synchronization

Approximating Nonlinear Functions With Latent Boundaries in Low-Rank Excitatory-Inhibitory Spiking Networks

Linear-Threshold Network Models for Describing and Analyzing Brain Dynamics

A State Space Approach for Piecewise-Linear Recurrent Neural Networks for Reconstructing Nonlinear Dynamics from Neural Measurements

Deep inverse modeling reveals dynamic-dependent invariances in neural circuit mechanisms

Identifying nonlinear dynamical systems via generative recurrent neural networks with applications to fMRI

Reconstruction of recurrent synaptic connectivity of thousands of neurons from simulated spiking activity

A scalable generative model for dynamical system reconstruction from neuroimaging data