Michael Schirner

The Virtual Brain Integrates Computational Modeling and Multimodal Neuroimaging

Brain function is thought to emerge from the interactions among neuronal populations. Apart from traditional efforts to reproduce brain dynamics from the micro- to macroscopic scales, complementary approaches develop phenomenological models of lower complexity. Such macroscopic models typically generate only a few selected-ideally functionally relevant-aspects of the brain dynamics. Importantly, they often allow an understanding of the underlying mechanisms beyond computational reproduction. Adding detail to these models will widen their ability to reproduce a broader range of dynamic features of the brain. For instance, such models allow for the exploration of consequences of focal and distributed pathological changes in the system, enabling us to identify and develop approaches to counteract those unfavorable processes. Toward this end, The Virtual Brain (TVB) ( www.thevirtualbrain.org ), a neuroinformatics platform with a brain simulator that incorporates a range of neuronal models and dynamics at its core, has been developed. This integrated framework allows the model-based simulation, analysis, and inference of neurophysiological mechanisms over several brain scales that underlie the generation of macroscopic neuroimaging signals. In this article, we describe how TVB works, and we present the first proof of concept.

Vascular endothelial growth factor (VEGF) receptor II-derived peptides inhibit VEGF.

Vascular endothelial growth factor (VEGF) directly stimulates endothelial cell proliferation and migration via tyrosine kinase receptors of the split kinase domain family. It mediates vascular growth and angiogenesis in the embryo but also in the adult in a variety of physiological and pathological conditions. The potential binding site of VEGF with its receptor was identified using cellulose-bound overlapping peptides of the extracytosolic part of the human vascular endothelial growth factor receptor II (VEGFR II). Thus, a peptide originating from the third globular domain of the VEGFR II comprising residues 247RTELNVGIDFNWEYP261was revealed as contiguous sequence stretch, which bound125I-VEGF165. A systematic replacement with L-amino acids within the peptide representing the putative VEGF-binding site on VEGFR II indicates Asp255 as the hydrophilic key residue for binding. The dimerized peptide (RTELNVGIDFNWEYPAS)2K inhibits VEGF165 binding with an IC50 of 0.5 μm on extracellular VEGFR II fragments and 30 μm on human umbilical vein cells. VEGF165-stimulated autophosphorylation of VEGFR II as well as proliferation and migration of microvascular endothelial cells was inhibited by the monomeric peptide RTELNVGIDFNWEYPASK at a half-maximal concentration of 3–10, 0.1, and 0.1 μm, respectively. We conclude that transduction of the VEGF165 signal can be interrupted with a peptide derived from the third Ig-like domain of VEGFR II by blockade of VEGF165 binding to its receptor.

An automated pipeline for constructing personalized virtual brains from multimodal neuroimaging data.

Large amounts of multimodal neuroimaging data are acquired every year worldwide. In order to extract high-dimensional information for computational neuroscience applications standardized data fusion and efficient reduction into integrative data structures are required. Such self-consistent multimodal data sets can be used for computational brain modeling to constrain models with individual measurable features of the brain, such as done with The Virtual Brain (TVB). TVB is a simulation platform that uses empirical structural and functional data to build full brain models of individual humans. For convenient model construction, we developed a processing pipeline for structural, functional and diffusion-weighted magnetic resonance imaging (MRI) and optionally electroencephalography (EEG) data. The pipeline combines several state-of-the-art neuroinformatics tools to generate subject-specific cortical and subcortical parcellations, surface-tessellations, structural and functional connectomes, lead field matrices, electrical source activity estimates and region-wise aggregated blood oxygen level dependent (BOLD) functional MRI (fMRI) time-series. The output files of the pipeline can be directly uploaded to TVB to create and simulate individualized large-scale network models that incorporate intra- and intercortical interaction on the basis of cortical surface triangulations and white matter tractograpy. We detail the pitfalls of the individual processing streams and discuss ways of validation. With the pipeline we also introduce novel ways of estimating the transmission strengths of fiber tracts in whole-brain structural connectivity (SC) networks and compare the outcomes of different tractography or parcellation approaches. We tested the functionality of the pipeline on 50 multimodal data sets. In order to quantify the robustness of the connectome extraction part of the pipeline we computed several metrics that quantify its rescan reliability and compared them to other tractography approaches. Together with the pipeline we present several principles to guide future efforts to standardize brain model construction. The code of the pipeline and the fully processed data sets are made available to the public via The Virtual Brain website (thevirtualbrain.org) and via github (https://github.com/BrainModes/TVB-empirical-data-pipeline). Furthermore, the pipeline can be directly used with High Performance Computing (HPC) resources on the Neuroscience Gateway Portal (http://www.nsgportal.org) through a convenient web-interface.

Structural architecture supports functional organization in the human aging brain at a regionwise and network level

Functional interactions in the brain are constrained by the underlying anatomical architecture, and structural and functional networks share network features such as modularity. Accordingly, age‐related changes of structural connectivity (SC) may be paralleled by changes in functional connectivity (FC). We provide a detailed qualitative and quantitative characterization of the SC–FC coupling in human aging as inferred from resting‐state blood oxygen‐level dependent functional magnetic resonance imaging and diffusion‐weighted imaging in a sample of 47 adults with an age range of 18–82. We revealed that SC and FC decrease with age across most parts of the brain and there is a distinct age‐dependency of regionwise SC–FC coupling and network‐level SC–FC relations. A specific pattern of SC–FC coupling predicts age more reliably than does regionwise SC or FC alone (r = 0.73, 95% CI = [0.7093, 0.8522]). Hence, our data propose that regionwise SC–FC coupling can be used to characterize brain changes in aging. Hum Brain Mapp 37:2645–2661, 2016.

How do parcellation size and short-range connectivity affect dynamics in large-scale brain network models?

Recent efforts to model human brain activity on the scale of the whole brain rest on connectivity estimates of large-scale networks derived from diffusion magnetic resonance imaging (dMRI). This type of connectivity describes white matter fiber tracts. The number of short-range cortico-cortical white-matter connections is, however, underrepresented in such large-scale brain models. It is still unclear on the one hand, which scale of representation of white matter fibers is optimal to describe brain activity on a large-scale such as recorded with magneto- or electroencephalography (M/EEG) or functional magnetic resonance imaging (fMRI), and on the other hand, to which extent short-range connections that are typically local should be taken into account. In this article we quantified the effect of connectivity upon large-scale brain network dynamics by (i) systematically varying the number of brain regions before computing the connectivity matrix, and by (ii) adding generic short-range connections. We used dMRI data from the Human Connectome Project. We developed a suite of preprocessing modules called SCRIPTS to prepare these imaging data for The Virtual Brain, a neuroinformatics platform for large-scale brain modeling and simulations. We performed simulations under different connectivity conditions and quantified the spatiotemporal dynamics in terms of Shannon Entropy, dwell time and Principal Component Analysis. For the reconstructed connectivity, our results show that the major white matter fiber bundles play an important role in shaping slow dynamics in large-scale brain networks (e.g. in fMRI). Faster dynamics such as gamma oscillations (around 40 Hz) are sensitive to the short-range connectivity if transmission delays are considered.

Bridging multiple scales in the human brain using computational modelling

Brain dynamics span multiple spatial and temporal scales, from fast spiking neurons to slow fluctuations over distributed areas. No single experimental method links data across scales. Here, we bridge this gap using The Virtual Brain connectome-based modelling platform to integrate multimodal data with biophysical models and support neurophysiological inference. Simulated cell populations were linked with subject-specific white-matter connectivity estimates and driven by electroencephalography-derived electric source activity. The models were fit to subject-specific resting-state functional magnetic resonance imaging data, and overfitting was excluded using 5-fold cross-validation. Further evaluation of the models show how balancing excitation with feedback inhibition generates an inverse relationship between α-rhythms and population firing on a faster time scale and resting-state network oscillations on a slower time scale. Lastly, large-scale interactions in the model lead to the emergence of scale-free power-law spectra. Our novel findings underscore the integrative role for computational modelling to complement empirical studies.

Differentiation of Alzheimer's disease based on local and global parameters in personalized Virtual Brain models

Alzheimers disease (AD) is marked by cognitive dysfunction emerging from neuropathological processes impacting brain function. AD affects brain dynamics at the local level, such as changes in the balance of inhibitory and excitatory neuronal populations, as well as long-range changes to the global network. Individual differences in these changes as they relate to behaviour are poorly understood. Here, we use a multi-scale neurophysiological model, “The Virtual Brain (TVB)”, based on empirical multi-modal neuroimaging data, to study how local and global dynamics correlate with individual differences in cognition. In particular, we modeled individual resting-state functional activity of 124 individuals across the behavioural spectrum from healthy aging, to amnesic Mild Cognitive Impairment (MCI), to AD. The model parameters required to accurately simulate empirical functional brain imaging data correlated significantly with cognition, and exceeded the predictive capacity of empirical connectomes.

Subject specificity of the correlation between large-scale structural and functional connectivity

Structural connectivity (SC), the physical pathways connecting regions in the brain, and functional connectivity (FC), the temporal coactivations, are known to be tightly linked. However, the nature of this relationship is still not understood. In the present study, we examined this relation more closely in six separate human neuroimaging datasets with different acquisition and preprocessing methods. We show that using simple linear associations, the relation between an individual’s SC and FC is not subject specific for five of the datasets. Subject specificity of SC-FC fit is achieved only for one of the six datasets, the multimodal Glasser Human Connectome Project (HCP) parcellated dataset. We show that subject specificity of SC-FC correspondence is limited across datasets due to relatively small variability between subjects in SC compared with the larger variability in FC. Author Summary We present evidence that, in most standard datasets, the subject variation in structural connectivity (SC) may be too weak to be reflected in the functional connectivity (FC) variability. However, subject specificity of SC-FC can be captured via fine, multimodally parcellated data because of greater SC variability across subjects. Nonetheless, SC and FC each show a large component that is common across subjects, which sets limitations on the extent of SC-FC subject specificity. Implications of these findings for personalized medicine should be considered. Namely, attention to the quality of processing and parcellation methods is critical for furthering our understanding of the relationship between individual SC and FC.