A Leadfield-Free Optimization Framework for Transcranially Applied Electric Currents.

Background: Transcranial Electrical Stimulation (TES), Temporal Interference Stimulation (TIS), Electroconvulsive Therapy (ECT) and Tumor Treating Fields (TTFields) are based on the application of electric current patterns to the brain.

Objective: The optimal electrode positions, shapes and alignments for generating a desired current pattern in the brain vary between persons due to anatomical variability. The aim is to develop a flexible and efficient computational approach to determine individually optimal montages based on electric field simulations.

Fast EEG/MEG BEM-based forward problem solution for high-resolution head models.

A fast BEM (boundary element method) based approach is developed to solve an EEG/MEG forward problem for a modern high-resolution head model. The method utilizes a charge-based BEM accelerated by the fast multipole method (BEM-FMM) with an adaptive mesh pre-refinement method (called b-refinement) close to the singular dipole source(s). No costly matrix-filling or direct solution steps typical for the standard BEM are required; the method generates on-skin voltages as well as MEG magnetic fields for high-resolution head models within 90 s after initial model assembly using a regular workstation.

Enabling electric field model of microscopically realistic brain.

Background: Modeling brain stimulation at the microscopic scale may reveal new paradigms for various stimulation modalities.

Objective: We present the largest map to date of extracellular electric field distributions within a layer L2/L3 mouse primary visual cortex brain sample. This was enabled by the automated analysis of serial section electron microscopy images with improved handling of image defects, covering a volume of 250 × 140 × 90 μm³.

Improving EEG Forward Modeling Using High-Resolution Five-Layer BEM-FMM Head Models: Effect on Source Reconstruction Accuracy.

Electroencephalographic (EEG) source localization is a fundamental tool for clinical diagnoses and brain-computer interfaces. We investigate the impact of model complexity on reconstruction accuracy by comparing the widely used three-layer boundary element method (BEM) as an inverse method against a five-layer BEM accelerated by the fast multipole method (BEM-FMM) and coupled with adaptive mesh refinement (AMR) as forward solver. Modern BEM-FMM with AMR can solve high-resolution multi-tissue models efficiently and accurately.

Statistical method accounts for microscopic electric field distortions around neurons when simulating activation thresholds.

Notwithstanding advances in computational models of neuromodulation, there are mismatches between simulated and experimental activation thresholds. Transcranial Magnetic Stimulation (TMS) of the primary motor cortex generates motor evoked potentials (MEPs). At the threshold of MEP generation, whole-head models predict macroscopic (at millimeter scale) electric fields (50-70 V/m) which are considerably below conventionally simulated cortical neuron thresholds (200-300 V/m).

Fast and Accurate EEG/MEG BEM-Based Forward Problem Solution for High-Resolution Head Models.

A BEM (boundary element method) based approach is developed to accurately solve an EEG/MEG forward problem for a modern high-resolution head model in approximately 60 seconds using a common workstation. The method utilizes a charge-based BEM with fast multipole acceleration (BEM-FMM) and a "smart" mesh pre-refinement (called -refinement) close to the singular source(s). No costly matrix-filling or direct solution steps typical for the standard BEM are required; the method generates on-skin voltages as well as MEG magnetic fields for high-resolution head models in approximately 60 seconds after initial model assembly.

Accuracy of dipole source reconstruction in the 3-layer BEM model against the 5-layer BEM-FMM model.

Objective: To compare cortical dipole fitting spatial accuracy between the widely used yet highly simplified 3-layer and modern more realistic 5-layer BEM-FMM models with and without (AMR) methods.

Methods: We generate simulated noiseless 256-channel EEG data from 5-layer (7-compartment) meshes of 15 subjects from the Connectome Young Adult dataset. For each subject, we test four dipole positions, three sets of conductivity values, and two types of head segmentation.

Enabling Electric Field Model of Microscopically Realistic Brain.

Unlabelled: Modeling brain stimulation at the microscopic scale may reveal new paradigms for various stimulation modalities. We present the largest map to date of extracellular electric field distributions within a layer L2/L3 mouse primary visual cortex brain sample. This was enabled by the automated analysis of serial section electron microscopy images with improved handling of image defects, covering a volume of 250 × 140 × 90 μm .

PyRates-A code-generation tool for modeling dynamical systems in biology and beyond.

The mathematical study of real-world dynamical systems relies on models composed of differential equations. Numerical methods for solving and analyzing differential equation systems are essential when complex biological problems have to be studied, such as the spreading of a virus, the evolution of competing species in an ecosystem, or the dynamics of neurons in the brain. Here we present PyRates, a Python-based software for modeling and analyzing differential equation systems via numerical methods.

Modeling the effects of transcranial magnetic stimulation on spatial attention.

. Transcranial magnetic stimulation (TMS) has been widely used to modulate brain activity in healthy and diseased brains, but the underlying mechanisms are not fully understood. Previous research leveraged biophysical modeling of the induced electric field (E-field) to map causal structure-function relationships in the primary motor cortex.

Laminar neural dynamics of auditory evoked responses: Computational modeling of local field potentials in auditory cortex of non-human primates.

Evoked neural responses to sensory stimuli have been extensively investigated in humans and animal models both to enhance our understanding of brain function and to aid in clinical diagnosis of neurological and neuropsychiatric conditions. Recording and imaging techniques such as electroencephalography (EEG), magnetoencephalography (MEG), local field potentials (LFPs), and calcium imaging provide complementary information about different aspects of brain activity at different spatial and temporal scales. Modeling and simulations provide a way to integrate these different types of information to clarify underlying neural mechanisms.

Spatiotemporal phase slip patterns for visual evoked potentials, covert object naming tasks, and insight moments extracted from 256 channel EEG recordings.

Phase slips arise from state transitions of the coordinated activity of cortical neurons which can be extracted from the EEG data. The phase slip rates (PSRs) were studied from the high-density (256 channel) EEG data, sampled at 16.384 kHz, of five adult subjects during covert visual object naming tasks.

Precise motor mapping with transcranial magnetic stimulation.

We describe a routine to precisely localize cortical muscle representations within the primary motor cortex with transcranial magnetic stimulation (TMS) based on the functional relation between induced electric fields at the cortical level and peripheral muscle activation (motor-evoked potentials; MEPs). Besides providing insights into structure-function relationships, this routine lays the foundation for TMS dosing metrics based on subject-specific cortical electric field thresholds. MEPs for different coil positions and orientations are combined with electric field modeling, exploiting the causal nature of neuronal activation to pinpoint the cortical origin of the MEPs.

White matter brain structure predicts language performance and learning success.

Individual differences in the ability to process language have long been discussed. Much of the neural basis of these, however, is yet unknown. Here we investigated the relationship between long-range white matter connectivity of the brain, as revealed by diffusion tractography, and the ability to process syntactically complex sentences in the participants' native language as well as the improvement thereof by multiday training.

Rate and oscillatory switching dynamics of a multilayer visual microcircuit model.

The neocortex is organized around layered microcircuits consisting of a variety of excitatory and inhibitory neuronal types which perform rate- and oscillation-based computations. Using modeling, we show that both superficial and deep layers of the primary mouse visual cortex implement two ultrasensitive and bistable switches built on mutual inhibitory connectivity motives between somatostatin, parvalbumin, and vasoactive intestinal polypeptide cells. The switches toggle pyramidal neurons between high and low firing rate states that are synchronized across layers through translaminar connectivity.

Comparison of the performance and reliability between improved sampling strategies for polynomial chaos expansion.

As uncertainty and sensitivity analysis of complex models grows ever more important, the difficulty of their timely realizations highlights a need for more efficient numerical operations. Non-intrusive Polynomial Chaos methods are highly efficient and accurate methods of mapping input-output relationships to investigate complex models. There is substantial potential to increase the efficacy of the method regarding the selected sampling scheme.

Cross-frequency phase-amplitude coupling in repetitive movements in patients with Parkinson's disease.

Bradykinesia is a cardinal motor symptom in Parkinson's disease (PD), the pathophysiology of which is not fully understood. We analyzed the role of cross-frequency coupling of oscillatory cortical activity in motor impairment in patients with PD and healthy controls. High-density EEG signals were recorded during various motor activities and at rest.

Modelling the effect of ephaptic coupling on spike propagation in peripheral nerve fibres.

Experimental and theoretical studies have shown that ephaptic coupling leads to the synchronisation and slowing down of spikes propagating along the axons within peripheral nerve bundles. However, the main focus thus far has been on a small number of identical axons, whereas realistic peripheral nerve bundles contain numerous axons with different diameters. Here, we present a computationally efficient spike propagation model, which captures the essential features of propagating spikes and their ephaptic interaction, and facilitates the theoretical investigation of spike volleys in large, heterogeneous fibre bundles.

The effect of meninges on the electric fields in TES and TMS. Numerical modeling with adaptive mesh refinement.

Background: When modeling transcranial electrical stimulation (TES) and transcranial magnetic stimulation (TMS) in the brain, the meninges - dura, arachnoid, and pia mater - are often neglected due to high computational costs.

Objective: We investigate the impact of the meningeal layers on the cortical electric field in TES and TMS while considering the headreco segmentation as the base model.

Method: We use T1/T2 MRI data from 16 subjects and apply the boundary element fast multipole method with adaptive mesh refinement, which enables us to accurately solve this problem and establish method convergence at reasonable computational cost.

Resting-State fMRI Functional Connectivity Strength Predicts Local Activity Change in the Dorsal Cingulate Cortex: A Multi-Target Focused rTMS Study.

Previous resting state functional magnetic resonance imaging (RS-fMRI) studies suggested that repetitive transcranial magnetic stimulation (rTMS) can modulate local activity in distant areas via functional connectivity (FC). A brain region has more than one connection with the superficial cortical areas. The current study proposed a multi-target focused rTMS protocol for indirectly stimulating a deep region, and to investigate 1) whether FC strength between stimulation targets (right middle frontal gyrus [rMFG] and right inferior parietal lobule [rIPL]) and effective region (dorsal anterior cingulate cortex [dACC]) can predict local activity changes of dACC and 2) whether multiple stimulation targets can focus on the dACC via FC.

Efficient high-resolution TMS mapping of the human motor cortex by nonlinear regression.

Transcranial magnetic stimulation (TMS) is a powerful tool to investigate causal structure-function relationships in the human brain. However, a precise delineation of the effectively stimulated neuronal populations is notoriously impeded by the widespread and complex distribution of the induced electric field. Here, we propose a method that allows rapid and feasible cortical localization at the individual subject level.

On the Role of Arkypallidal and Prototypical Neurons for Phase Transitions in the External Pallidum.

The external pallidum (globus pallidus pars externa [GPe]) plays a central role for basal ganglia functions and dynamics and, consequently, has been included in most computational studies of the basal ganglia. These studies considered the GPe as a homogeneous neural population. However, experimental studies have shown that the GPe contains at least two distinct cell types (prototypical and arkypallidal cells).

Functional brain plasticity during L1 training on complex sentences: Changes in gamma-band oscillatory activity.

The adult human brain remains plastic even after puberty. However, whether first language (L1) training in adults can alter the language network is yet largely unknown. Thus, we conducted a longitudinal training experiment on syntactically complex German sentence comprehension.