Closed-loop transcranial magnetic stimulation of real-time EEG based on the AR mode method

The applied principles of EEG analysis methods in neuroscience and clinical neurology

Electroencephalography (EEG) is a non-invasive measurement method for brain activity. Due to its safety, high resolution, and hypersensitivity to dynamic changes in brain neural signals, EEG has aroused much interest in scientific research and medical fields. This article reviews the types of EEG signals, multiple EEG signal analysis methods, and the application of relevant methods in the neuroscience field and for diagnosing neurological diseases. First, three types of EEG signals, including time-invariant EEG, accurate event-related EEG, and random event-related EEG, are introduced. Second, five main directions for the methods of EEG analysis, including power spectrum analysis, time-frequency analysis, connectivity analysis, source localization methods, and machine learning methods, are described in the main section, along with different sub-methods and effect evaluations for solving the same problem. Finally, the application scenarios of different EEG analysis methods are emphasized, and the advantages and disadvantages of similar methods are distinguished. This article is expected to assist researchers in selecting suitable EEG analysis methods based on their research objectives, provide references for subsequent research, and summarize current issues and prospects for the future.

Identifiability Analysis and Noninvasive Online Estimation of the First-Order Neural Activation Dynamics in the Brain With Closed-Loop Transcranial Magnetic Stimulation

BACKGROUND

Neurons demonstrate very distinct nonlinear activation dynamics, influenced by the neuron type, morphology, ion channel expression, and various other factors. The measurement of the activation dynamics can identify the neural target of stimulation and detect deviations, e.g., for diagnosis. This paper describes a tool for closed-loop sequential parameter estimation (SPE) of the activation dynamics through transcranial magnetic stimulation (TMS). The proposed SPE method operates in real time, selects ideal stimulus parameters, detects and processes the response, and concurrently estimates the input-output (IO) curve and the first-order approximation of the activated neural target.

OBJECTIVE

To develop a method for concurrent SPE of the first-order activation dynamics and IO curve with closed-loop TMS.

METHOD

First, identifiability of an integrated model of the first-order neural activation dynamics and IO curve is assessed, demonstrating that at least two IO curves need to be acquired with different pulse widths. Then, a two-stage SPE method is proposed. It estimates the IO curve by using Fisher information matrix (FIM) optimization in the first stage and subsequently estimates the membrane time constant as well as the coupling gain in the second stage. The procedure continues in a sequential manner until a stopping rule is satisfied.

RESULTS

The results of 73 simulation cases confirm the satisfactory estimation of the membrane time constant and coupling gain with average absolute relative errors (AREs) of 6.2% and 5.3%, respectively, with an average of 344 pulses (172 pulses for each IO curve or pulse width). The method estimates the IO curves' lower and upper plateaus, mid-point, and slope with average AREs of 0.2%, 0.7%, 0.9%, and 14.5%, respectively. The conventional time constant estimation method based on the strength-duration (S-D) curve leads to 33.3% ARE, which is 27.0% larger than 6.2% ARE obtained through the proposed real-time FIM-based SPE method in this paper.

CONCLUSIONS

SPE of the activation dynamics requires acquiring at least two IO curves with different pulse widths, which needs a controllable TMS (cTMS) device with adjustable pulse duration.

SIGNIFICANCE

The proposed SPE method enhances the cTMS functionality, which can contribute novel insights in research and clinical studies.

Closed-loop optimal and automatic tuning of pulse amplitude and width in EMG-guided controllable transcranial magnetic stimulation

This paper proposes an efficient algorithm for automatic and optimal tuning of pulse amplitude and width for sequential parameter estimation (SPE) of the neural membrane time constant and input-output (IO) curve parameters in closed-loop electromyography-guided (EMG-guided) controllable transcranial magnetic stimulation (cTMS). The proposed SPE is performed by administering a train of optimally tuned TMS pulses and updating the estimations until a stopping rule is satisfied or the maximum number of pulses is reached. The pulse amplitude is computed by the Fisher information maximization. The pulse width is chosen by maximizing a normalized depolarization factor, which is defined to separate the optimization and tuning of the pulse amplitude and width. The normalized depolarization factor maximization identifies the critical pulse width, which is an important parameter in the identifiability analysis, without any prior neurophysiological or anatomical knowledge of the neural membrane. The effectiveness of the proposed algorithm is evaluated through simulation. The results confirm satisfactory estimation of the membrane time constant and IO curve parameters for the simulation case. By defining the stopping rule based on the satisfaction of the convergence criterion with tolerance of 0.01 for 5 consecutive times for all parameters, the IO curve parameters are estimated with 52 TMS pulses, with absolute relative estimation errors (AREs) of less than 7%. The membrane time constant is estimated with 0.67% ARE, and the pulse width value tends to the critical pulse width with 0.16% ARE with 52 TMS pulses. The results confirm that the pulse width and amplitude can be tuned optimally and automatically to estimate the membrane time constant and IO curve parameters in real-time with closed-loop EMG-guided cTMS.

Optimization of phase prediction for brain-state dependent stimulation: a grid-search approach

Objective.Sources of heterogeneity in non-invasive brain stimulation literature can be numerous, with underlying brain states and protocol differences at the top of the list. Yet, incoherent results from brain-state-dependent stimulation experiments suggest that there are further factors adding to the variance. Hypothesizing that different signal processing pipelines might be partly responsible for heterogeneity; we investigated their effects on brain-state forecasting approaches.Approach.A grid-search was used to determine the fastest and most-accurate combination of preprocessing parameters and phase-forecasting algorithms. The grid-search was applied on a synthetic dataset and validated on electroencephalographic (EEG) data from a healthy (n= 18) and stroke (n= 31) cohort.Main results.Differences in processing pipelines led to different results; the grid-search chosen pipelines significantly increased the accuracy of published forecasting methods. The accuracy achieved in healthy was comparably high in stroke patients.Significance.This systematic offline analysis highlights the importance of the specific EEG processing and forecasting pipelines used for online state-dependent setups where precision in phase prediction is critical. Moreover, successful results in the stroke cohort pave the way to test state-dependent interventional treatment approaches.

Transcranial magnetic stimulation treatment in Alzheimer's disease: a meta-analysis of its efficacy as a function of protocol characteristics and degree of personalization

Alzheimer's disease (AD) represents the most common type of neurodegenerative disorder. Although our knowledge on the causes of AD remains limited and no curative treatments are available, several interventions have been proposed in trying to improve patients' symptomatology. Among those, transcranial magnetic stimulation (TMS) has been shown a promising, safe and noninvasive intervention to improve global cognitive functioning. Nevertheless, we currently lack agreement between research studies on the optimal stimulation protocol yielding the highest efficacy in these patients. To answer this query, we conducted a systematic literature search in PubMed, PsycINFO and Scopus databases and meta-analysis of studies published in the last 10 years (2010-2021) according to the Preferred Reporting Items for Systematic Reviews and Meta-Analyses (PRISMA) guidelines. Differently from prior published meta-analytic work, we investigated whether protocols that considered participants-specific neuroimaging scans for the selection of individualized stimulation targets held more successful outcomes compared to those relying on a generalized targeting selection criteria. We then compared the effect sizes of subsets of studies based on additional protocol characteristics (frequency, duration of intervention, number of stimulation sites, use of concomitant cognitive training and patients' educational level). Our results confirm TMS efficacy in improving global cognitive functioning in mild-to-moderate AD patients, but also highlight the flaws of current protocols characteristics, including a possible lack of sufficient personalization in stimulation protocols.

A formalism for sequential estimation of neural membrane time constant and input--output curve towards selective and closed-loop transcranial magnetic stimulation

OBJECTIVE

To obtain a formalism for real-time concurrent sequential estimation of neural membrane time constant and input--output (IO) curve with transcranial magnetic stimulation (TMS).

APPROACH

First, the neural membrane response and depolarization factor, which leads to motor evoked potentials (MEPs) with TMS are analytically computed and discussed. Then, an integrated model is developed which combines the neural membrane time constant and input--output curve. Identifiability of the proposed integrated model is discussed. A condition is derived, which assures estimation of the proposed integrated model. Finally, sequential parameter estimation (SPE) of the neural membrane time constant and IO curve is described through closed-loop optimal sampling and open-loop uniform sampling TMS. Without loss of generality, this paper focuses on a specific case of commercialized TMS pulse shapes. The proposed formalism and SPE method are directly applicable to other pulse shapes.

MAIN RESULTS

The results confirm satisfactory estimation of the membrane time constant and IO curve parameters. By defining a stopping rule based on five times consecutive convergence of the estimation parameters with a tolerances of 0.01, the membrane time constant and IO curve parameters are estimated with 82 TMS pulses with absolute relative estimation errors (AREs) of less than 4% with the optimal sampling SPE method. At this point, the uniform sampling SPE method leads to AREs up to 16%. The uniform sampling method does not satisfy the stopping rule due to the large estimation variations.

SIGNIFICANCE

This paper provides a tool for real-time closed-loop SPE of the neural time constant and IO curve, which can contribute novel insights in TMS studies. SPE of the membrane time constant enables selective stimulation, which can be used for advanced brain research, precision medicine and personalized medicine.

TAP: targeting and analysis pipeline for optimization and verification of coil placement in transcranial magnetic stimulation

Objective.Transcranial magnetic stimulation (TMS) can modulate brain function via an electric field (E-field) induced in a brain region of interest (ROI). The ROI E-field can be computationally maximized and set to match a specific reference using individualized head models to find the optimal coil placement and stimulus intensity. However, the available software lacks many practical features for prospective planning of TMS interventions and retrospective evaluation of the experimental targeting accuracy.Approach.The TMS targeting and analysis pipeline (TAP) software uses an MRI/fMRI-derived brain target to optimize coil placement considering experimental parameters such as the subject's hair thickness and coil placement restrictions. The coil placement optimization is implemented in SimNIBS 3.2, for which an additional graphical user interface (TargetingNavigator) is provided to visualize/adjust procedural parameters. The coil optimization process also computes the E-field at the target, allowing the selection of the TMS device intensity setting to achieve specific E-field strengths. The optimized coil placement information is prepared for neuronavigation software, which supports targeting during the TMS procedure. The neuronavigation system can record the coil placement during the experiment, and these data can be processed in TAP to quantify the accuracy of the experimental TMS coil placement and induced E-field.Main results.TAP was demonstrated in a study consisting of three repetitive TMS sessions in five subjects. TMS was delivered by an experienced operator under neuronavigation with the computationally optimized coil placement. Analysis of the experimental accuracy from the recorded neuronavigation data indicated coil location and orientation deviations up to about 2 mm and 2°, respectively, resulting in an 8% median decrease in the target E-field magnitude compared to the optimal placement.Significance.TAP supports navigated TMS with a variety of features for rigorous and reproducible stimulation delivery, including planning and evaluation of coil placement and intensity selection for E-field-based dosing.