Precision data-driven modeling of cortical dynamics reveals idiosyncratic mechanisms underlying canonical oscillations

Task-free brain activity affords unique insight into the functional structure of brain network dynamics and is a strong marker of individual differences. In this work, we present an algorithmic optimization framework that makes it possible to directly invert and parameterize brain-wide dynamical-systems models involving hundreds of interacting brain areas, from single-subject time-series recordings. This technique provides a powerful neurocomputational tool for interrogating mechanisms underlying individual brain dynamics ("precision brain models") and making quantitative predictions. We extensively validate the models' performance in forecasting future brain activity and predicting individual variability in key M/EEG markers. Lastly, we demonstrate the power of our technique in resolving individual differences in the generation of alpha and beta-frequency oscillations. We characterize subjects based upon model attractor topology and a dynamical-systems mechanism by which these topologies generate individual variation in the expression of alpha vs. beta rhythms. We trace these phenomena back to global variation in excitation-inhibition balance, highlighting the explanatory power of our framework in generating mechanistic insights.

Brain model state space reconstruction using an LSTM neural network

Objective
Kalman filtering has previously been applied to track neural model states and parameters, particularly at the scale relevant to electroencephalography (EEG). However, this approach lacks a reliable method to determine the initial filter conditions and assumes that the distribution of states remains Gaussian. This study presents an alternative, data-driven method to track the states and parameters of neural mass models (NMMs) from EEG recordings using deep learning techniques, specifically a Long Short-Term Memory (LSTM) neural network.
Approach
An LSTM filter was trained on simulated EEG data generated by a neural mass model using a wide range of parameters. With an appopriately customised loss function, the LSTM filter can learn the behaviour of NMMs. As a result, it can output the state vector and parameters of NMMs given observation data as the input.
Main Results
Test results using simulated data yielded correlations with R squared of around 0.99 and verified that the method is robust to noise and can be more accurate than a nonlinear Kalman filter when the initial conditions of the Kalman filter are not accurate. As an example of real-world application, the LSTM filter was also applied to real EEG data that included epileptic seizures, and revealed changes in connectivity strength parameters at the beginnings of seizures.
Significance
Tracking the state vector and parameters of mathematical brain models is of great importance in the area of brain modelling, monitoring, imaging and control. This approach has no need to specify the initial state vector and parameters, which is very difficult to do in practice because many of the variables being estimated cannot be measured directly in physiological experiments. This method may be applied using any neural mass model and, therefore, provides a general, novel, efficient approach to estimate brain model variables that are often difficult to measure.&#xD.

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

This paper introduces methods and a novel toolbox that efficiently integrates high-dimensional Neural Mass Models (NMMs) specified by two essential components. The first is the set of nonlinear Random Differential Equations (RDEs) of the dynamics of each neural mass. The second is the highly sparse three-dimensional Connectome Tensor (CT) that encodes the strength of the connections and the delays of information transfer along the axons of each connection. To date, simplistic assumptions prevail about delays in the CT, often assumed to be Dirac-delta functions. In reality, delays are distributed due to heterogeneous conduction velocities of the axons connecting neural masses. These distributed-delay CTs are challenging to model. Our approach implements these models by leveraging several innovations. Semi-analytical integration of RDEs is done with the Local Linearization (LL) scheme for each neural mass model, ensuring dynamical fidelity to the original continuous-time nonlinear dynamic. This semi-analytic LL integration is highly computationally-efficient. In addition, a tensor representation of the CT facilitates parallel computation. It also seamlessly allows modeling distributed delays CT with any level of complexity or realism. This ease of implementation includes models with distributed-delay CTs. Consequently, our algorithm scales linearly with the number of neural masses and the number of equations they are represented with, contrasting with more traditional methods that scale quadratically at best. To illustrate the toolbox's usefulness, we simulate a single Zetterberg-Jansen-Rit (ZJR) cortical column, a single thalmo-cortical unit, and a toy example comprising 1000 interconnected ZJR columns. These simulations demonstrate the consequences of modifying the CT, especially by introducing distributed delays. The examples illustrate the complexity of explaining EEG oscillations, e.g., split alpha peaks, since they only appear for distinct neural masses. We provide an open-source Script for the toolbox.

Intensity dependence of sub-harmonics in cortical response to photic stimulation

Objective. Periodic photic stimulation of human volunteers at 10 Hz is known to entrain their Electroencephalography (EEG) signals. This entrainment manifests as an increment in power at 10, 20, 30 Hz. We observed that this entrainment is accompanied by the emergence of sub-harmonics, but only at specific frequencies and higher intensities of the stimulating signal. Thereafter, we describe our results and explain them using the physiologically inspired Jansen and Rit Neural Mass Model (NMM).Approach. Four human volunteers were separately exposed to both high and low intensity 10 Hz and 6 Hz stimulation. A total of 4 experiments per subject were therefore performed. Simulations and bifurcation analysis of the NMM were carried out and compared with the experimental findings. <i> Main results. High intensity 10 Hz stimulation led to an increment in power at 5 Hz across all the 4 subjects. No increment of power was observed with low intensity stimulation. However, when the same protocol was repeated with a 6 Hz photic stimulation, neither high nor low intensity stimulation were found to cause a discernible change in power at 3 Hz. We found that the NMM was able to recapitulate these results. A further numerical analysis indicated that this arises from the underlying bifurcation structure of the NMM. <i> Significance. The excellent match between theory and experiment suggest that the bifurcation properties of the NMM are mirroring similar features possessed by the actual neural masses producing the EEG dynamics. Neural Mass Models could thus be valuable for understanding properties and pathologies of EEG dynamics, and may contribute to the engineering of brain-computer interface technologies.