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.

Neural mass modelling for the masses: Democratising access to whole-brain biophysical modelling with FastDMF

A neural mass model with neuromodulation

Parameter Estimation in Brain Dynamics Models from Resting-State fMRI Data using Physics-Informed Neural Networks

On the intersection between data quality and dynamical modelling of large-scale fMRI signals

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

Biophysically inspired mean-field model of neuronal populations driven by ion exchange mechanisms

neurolib: A Simulation Framework for Whole-Brain Neural Mass Modeling

Physiologically informed dynamic causal modeling of fMRI data

Accurate and Efficient Simulation of Very High-Dimensional Neural Mass Models with Distributed-Delay Connectome Tensors

Data-driven modelling of brain activity using neural networks, Diffusion Maps, and the Koopman operator

Realistic head modeling of electromagnetic brain activity: An integrated Brainstorm pipeline from MRI data to the FEM solution

Modeling functional resting-state brain networks through neural message passing on the human connectome

Biomechanical modelling of brain atrophy through deep learning

RateML: A Code Generation Tool for Brain Network Models

A generative model of realistic brain cells with application to numerical simulation of diffusion-weighted MR signal

Increasing spectral DCM flexibility and speed by leveraging Julia's ModelingToolkit and automated differentiation

Multiple-shooting adjoint method for whole-brain dynamic causal modeling

Generalized Recurrent Neural Network accommodating Dynamic Causal Modeling for functional MRI analysis

Full Bayesian Modeling for fMRI Group Analysis

Deep coupled registration and segmentation of multimodal whole-brain images