Abstract:Biophysically detailed neural models are a powerful technique to study neural dynamics in health and disease with a growing number of established and openly available models. A major challenge in the use of such models is that parameter inference is an inherently difficult and unsolved problem. Identifying unique parameter distributions that can account for observed neural dynamics, and differences across experimental conditions, is essential to their meaningful use. Recently, simulation based inference (SBI) has been proposed as an approach to perform Bayesian inference to estimate parameters in detailed neural models. SBI overcomes the challenge of not having access to a likelihood function, which has severely limited inference methods in such models, by leveraging advances in deep learning to perform density estimation. While the substantial methodological advancements offered by SBI are promising, their use in large scale biophysically detailed models is challenging and methods for doing so have not been established, particularly when inferring parameters that can account for time series waveforms. We provide guidelines and considerations on how SBI can be applied to estimate time series waveforms in biophysically detailed neural models starting with a simplified example and extending to specific applications to common MEG/EEG waveforms using the the large scale neural modeling framework of the Human Neocortical Neurosolver. Specifically, we describe how to estimate and compare results from example oscillatory and event related potential simulations. We also describe how diagnostics can be used to assess the quality and uniqueness of the posterior estimates. The methods described provide a principled foundation to guide future applications of SBI in a wide variety of applications that use detailed models to study neural dynamics. A central problem in computational neural modeling is estimating model parameters that can account for observed activity patterns. While several techniques exist to perform parameter inference in special classes of abstract neural models, there are comparatively few approaches for large scale biophysically detailed neural models. In this work, we describe challenges and solutions in applying a deep learning based statistical framework to estimate parameters in a biophysically detailed large scale neural model, and emphasize the particular difficulties in estimating parameters for time series data. Our example uses a multi-scale model designed to connect human MEG/EEG recordings to the underlying cell and circuit level generators. Our approach allows for crucial insight into how cell-level properties interact to produce measured neural activity, and provides guidelines for diagnosing the quality of the estimate and uniqueness of predictions for different MEG/EEG biomarkers.

Probabilistic Inference on Virtual Brain Models of Disorders

Simulation-Based Inference for Whole-Brain Network Modeling of Epilepsy using Deep Neural Density Estimators

Principles and Operation of Virtual Brain Twins

Dynamic Causal Modeling in Probabilistic Programming Languages

Methods and considerations for estimating parameters in biophysically detailed neural models with simulation based inference

On the influence of prior information evaluated by fully Bayesian criteria in a personalized whole-brain model of epilepsy spread

Network inference in a stochastic multi-population neural mass model via approximate Bayesian computation

Fifty shades of The Virtual Brain: Converging optimal working points yield biologically plausible electrophysiological and imaging features

Virtual brain twins: from basic neuroscience to clinical use

Modeling brain dynamics after tumor resection using The Virtual Brain

Simulation-based inference for efficient identification of generative models in computational connectomics

Parameter estimation in a whole-brain network model of epilepsy: Comparison of parallel global optimization solvers

Probabilistic prediction of neurological disorders with a statistical assessment of neuroimaging data modalities

Learning Large Number of Local Statistical Models Via Variational Bayesian Inference for Brain Voxel Classification in Magnetic Resonance Images

Model-based whole-brain perturbational landscape of neurodegenerative diseases

Whole-Brain Network Models: From Physics to Bedside

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

Subject-specific features of excitation/inhibition profiles in neurodegenerative diseases

Bayesian Models of Functional Connectomics and Behavior

Predicting Brain Activity Using a Bayesian Spatial Model.

A Deep Probabilistic Spatiotemporal Framework for Dynamic Graph Representation Learning with Application to Brain Disorder Identification