Focused Ultrasound for Treatment of Movement Disorders: A Review of Non-Food and Drug Administration Approved Indications

INTRODUCTION

MRI-guided focused ultrasound (FUS) is an incisionless thermo-ablative procedure that may be used to treat medication-refractory movement disorders, with a growing number of potential anatomic targets and clinical applications. As of this article's publication, the only US Food and Drug Administration (FDA)-approved uses of FUS for movement disorders are thalamotomy for essential tremor (ET) and tremor-dominant Parkinson's Disease (PD), and pallidotomy for other cardinal symptoms of PD. We present a state-of-the-art review on all non-FDA approved indications of FUS for movement disorders, beyond the most well-described indications of ET and PD. Our objective was to summarize the safety and efficacy of FUS in this setting and provide a roadmap for future directions of FUS for movement disorders.

METHODS

A state-of-the-art review was conducted on use of FUS for non-FDA approved movement disorders. All movement disorders excluding FDA-approved uses for ET and PD were included.

RESULTS

A total of 25 studies on 172 patients were included. In patients with tremor plus dystonia syndromes (n = 6), ventralis intermediate nucleus of the thalamus (VIM)-FUS gave >50% tremor reduction, with no improvement in dystonia and worsened dystonia in 2/6 patients. Ventral-oralis complex (VO)-FUS gave >50% improvement for focal hand dystonia (n = 6) and 100% return to musical performance in musician's dystonia (n = 6). In patients with multiple sclerosis (MS) and tremor (n = 3), improvement in tremor was seen in 2 patients with a favorable skull density ratio; no MS disease change was noted after VIM-FUS. In patients with tremor and comorbid ataxia syndromes (n = 3), none were found to have worsened ataxia after VIM-FUS; all had clinically significant tremor improvement. Subthalamic nucleus (STN)-FUS for PD (n = 49) gave approximately 50% improvement in PD motor symptoms, with dystonia and mild dyskinesias as possible adverse effects. Cerebellothalamic tract (CTT-FUS) for ET (n = 42) gave 55-90% tremor improvement, with gait dysfunction as a rare persistent adverse effect. Pallidothalamic tract (PTT-FUS) for PD (n = 50) gave approximately 50% improvement in motor symptoms, with mild speech dysfunction as a possible adverse effect.

CONCLUSION

VIM-FUS appeared safe and effective for heterogenous tremor etiologies, and VO-FUS appeared most effective for isolated segmental dystonia. STN-FUS was effective for PD symptom reduction; postoperative dystonia and mild on-medication dyskinesias required medical management. Tractography-based targeting with CTT-FUS for ET and PTT-FUS for PD demonstrated promising early results. Larger prospective trials with long-term follow-up are needed to the evaluate the safety and efficacy non-FDA approved indications for FUS.

The seizure severity score: a quantitative tool for comparing seizures and their response to therapy

Objective.Epilepsy is a neurological disorder characterized by recurrent seizures which vary widely in severity, from clinically silent to prolonged convulsions. Measuring severity is crucial for guiding therapy, particularly when complete control is not possible. Seizure diaries, the current standard for guiding therapy, are insensitive to the duration of events or the propagation of seizure activity across the brain. We present a quantitative seizure severity score that incorporates electroencephalography (EEG) and clinical data and demonstrate how it can guide epilepsy therapies.Approach.We collected intracranial EEG and clinical semiology data from 54 epilepsy patients who had 256 seizures during invasive, in-hospital presurgical evaluation. We applied an absolute slope algorithm to EEG recordings to identify seizing channels. From this data, we developed a seizure severity score that combines seizure duration, spread, and semiology using non-negative matrix factorization. For validation, we assessed its correlation with independent measures of epilepsy burden: seizure types, epilepsy duration, a pharmacokinetic model of medication load, and response to epilepsy surgery. We investigated the association between the seizure severity score and preictal network features.Main results.The seizure severity score augmented clinical classification by objectively delineating seizure duration and spread from recordings in available electrodes. Lower preictal medication loads were associated with higher seizure severity scores (p= 0.018, 97.5% confidence interval = [-1.242, -0.116]) and lower pre-surgical severity was associated with better surgical outcome (p= 0.042). In 85% of patients with multiple seizure types, greater preictal change from baseline was associated with higher severity.Significance.We present a quantitative measure of seizure severity that includes EEG and clinical features, validated on gold standard in-patient recordings. We provide a framework for extending our tool's utility to ambulatory EEG devices, for linking it to seizure semiology measured by wearable sensors, and as a tool to advance data-driven epilepsy care.

Quantifying interictal intracranial EEG to predict focal epilepsy

Introduction

Intracranial EEG (IEEG) is used for 2 main purposes, to determine: (1) if epileptic networks are amenable to focal treatment and (2) where to intervene. Currently these questions are answered qualitatively and sometimes differently across centers. There is a need for objective, standardized methods to guide surgical decision making and to enable large scale data analysis across centers and prospective clinical trials.

Methods

We analyzed interictal data from 101 patients with drug resistant epilepsy who underwent presurgical evaluation with IEEG. We chose interictal data because of its potential to reduce the morbidity and cost associated with ictal recording. 65 patients had unifocal seizure onset on IEEG, and 36 were non-focal or multi-focal. We quantified the spatial dispersion of implanted electrodes and interictal IEEG abnormalities for each patient. We compared these measures against the "5 Sense Score (5SS)," a pre-implant estimate of the likelihood of focal seizure onset, and assessed their ability to predict the clinicians' choice of therapeutic intervention and the patient outcome.

Results

The spatial dispersion of IEEG electrodes predicted network focality with precision similar to the 5SS (AUC = 0.67), indicating that electrode placement accurately reflected pre-implant information. A cross-validated model combining the 5SS and the spatial dispersion of interictal IEEG abnormalities significantly improved this prediction (AUC = 0.79; p<0.05). The combined model predicted ultimate treatment strategy (surgery vs. device) with an AUC of 0.81 and post-surgical outcome at 2 years with an AUC of 0.70. The 5SS, interictal IEEG, and electrode placement were not correlated and provided complementary information.

Conclusions

Quantitative, interictal IEEG significantly improved upon pre-implant estimates of network focality and predicted treatment with precision approaching that of clinical experts. We present this study as an important step in building standardized, quantitative tools to guide epilepsy surgery.

Remote effects of temporal lobe epilepsy surgery: Long-term morphological changes after surgical resection

OBJECTIVE

Epilepsy surgery is an effective treatment for drug-resistant patients. However, how different surgical approaches affect long-term brain structure remains poorly characterized. Here, we present a semiautomated method for quantifying structural changes after epilepsy surgery and compare the remote structural effects of two approaches, anterior temporal lobectomy (ATL), and selective amygdalohippocampectomy (SAH).

METHODS

We studied 36 temporal lobe epilepsy patients who underwent resective surgery (ATL = 22, SAH = 14). All patients received same-scanner MR imaging preoperatively and postoperatively (mean 2 years). To analyze postoperative structural changes, we segmented the resection zone and modified the Advanced Normalization Tools (ANTs) longitudinal cortical pipeline to account for resections. We compared global and regional annualized cortical thinning between surgical treatments.

RESULTS

Across procedures, there was significant cortical thinning in the ipsilateral insula, fusiform, pericalcarine, and several temporal lobe regions outside the resection zone as well as the contralateral hippocampus. Additionally, increased postoperative cortical thickness was seen in the supramarginal gyrus. Patients treated with ATL exhibited greater annualized cortical thinning compared with SAH cases (ATL: -0.08 ± 0.11 mm per year, SAH: -0.01 ± 0.02 mm per year, t = 2.99, P = 0.006). There were focal postoperative differences between the two treatment groups in the ipsilateral insula (P = 0.039, corrected). Annualized cortical thinning rates correlated with preoperative cortical thickness (r = 0.60, P < 0.001) and had weaker associations with age at surgery (r = -0.33, P = 0.051) and disease duration (r = -0.42, P = 0.058).

SIGNIFICANCE

Our evidence suggests that selective procedures are associated with less cortical thinning and that earlier surgical intervention may reduce long-term impacts on brain structure.

Quantitative approaches to guide epilepsy surgery from intracranial EEG

Over the past 10 years, the drive to improve outcomes from epilepsy surgery has stimulated widespread interest in methods to quantitatively guide epilepsy surgery from intracranial EEG (iEEG). Many patients fail to achieve seizure freedom, in part due to the challenges in subjective iEEG interpretation. To address this clinical need, quantitative iEEG analytics have been developed using a variety of approaches, spanning studies of seizures, interictal periods, and their transitions, and encompass a range of techniques including electrographic signal analysis, dynamical systems modeling, machine learning and graph theory. Unfortunately, many methods fail to generalize to new data and are sensitive to differences in pathology and electrode placement. Here, we critically review selected literature on computational methods of identifying the epileptogenic zone from iEEG. We highlight shared methodological challenges common to many studies in this field and propose ways that they can be addressed. One fundamental common pitfall is a lack of open-source, high-quality data, which we specifically address by sharing a centralized high-quality, well-annotated, multicentre dataset consisting of >100 patients to support larger and more rigorous studies. Ultimately, we provide a road map to help these tools reach clinical trials and hope to improve the lives of future patients.

A pharmacokinetic model of antiseizure medication load to guide care in the epilepsy monitoring unit

OBJECTIVE

Evaluating patients with drug-resistant epilepsy often requires inducing seizures by tapering antiseizure medications (ASMs) in the epilepsy monitoring unit (EMU). The relationship between ASM taper strategy, seizure timing, and severity remains unclear. In this study, we developed and validated a pharmacokinetic model of total ASM load and tested its association with seizure occurrence and severity in the EMU.

METHODS

We studied 80 patients who underwent intracranial electroencephalographic recording for epilepsy surgery planning. We developed a first order pharmacokinetic model of the ASMs administered in the EMU to generate a continuous metric of overall ASM load. We then related modeled ASM load to seizure likelihood and severity. We determined the association between the rate of ASM load reduction, the length of hospital stay, and the probability of having a severe seizure. Finally, we used modeled ASM load to predict oncoming seizures.

RESULTS

Seizures occurred in the bottom 50th percentile of sampled ASM loads across the cohort (p < .0001, Wilcoxon signed-rank test), and seizures requiring rescue therapy occurred at lower ASM loads than seizures that did not require rescue therapy (logistic regression mixed effects model, odds ratio = .27, p = .01). Greater ASM decrease early in the EMU was not associated with an increased likelihood of having a severe seizure, nor with a shorter length of stay.

SIGNIFICANCE

A pharmacokinetic model can accurately estimate ASM levels for patients in the EMU. Lower modeled ASM levels are associated with increased seizure likelihood and seizure severity. We show that ASM load, rather than ASM taper speed, is associated with severe seizures. ASM modeling has the potential to help optimize taper strategy to minimize severe seizures while maximizing diagnostic yield.

Addressing spatial bias in intracranial EEG functional connectivity analyses for epilepsy surgical planning

Objective.To determine the effect of epilepsy on intracranial electroencephalography (EEG) functional connectivity, and the ability of functional connectivity to localize the seizure onset zone (SOZ), controlling for spatial biases.Approach.We analyzed intracranial EEG data from patients with drug-resistant epilepsy admitted for pre-surgical planning. We calculated intracranial EEG functional networks and determined whether changes in functional connectivity lateralized the SOZ using a spatial subsampling method to control for spatial bias. We developed a 'spatial null model' to localize the SOZ electrode using only spatial sampling information, ignoring EEG data. We compared the performance of this spatial null model against models incorporating EEG functional connectivity and interictal spike rates.Main results.About 110 patients were included in the study, although the number of patients differed across analyses. Controlling for spatial sampling, the average connectivity was lower in the SOZ region relative to the same anatomic region in the contralateral hemisphere. A model using intra-hemispheric connectivity accurately lateralized the SOZ (average accuracy 75.5%). A spatial null model incorporating spatial sampling information alone achieved moderate accuracy in classifying SOZ electrodes (mean AUC = 0.70, 95% CI 0.63-0.77). A model incorporating intracranial EEG functional connectivity and spike rate data further outperformed this spatial null model (AUC 0.78,p= 0.002 compared to spatial null model). However, a model incorporating functional connectivity without spike rate data did not significantly outperform the null model (AUC 0.72,p= 0.38).Significance.Intracranial EEG functional connectivity is reduced in the SOZ region, and interictal data predict SOZ electrode localization and laterality, however a predictive model incorporating functional connectivity without interictal spike rates did not significantly outperform a spatial null model. We propose constructing a spatial null model to provide an estimate of the pre-implant hypothesis of the SOZ, and to serve as a benchmark for further machine learning algorithms in order to avoid overestimating model performance because of electrode sampling alone.

Normative intracranial EEG maps epileptogenic tissues in focal epilepsy

Planning surgery for patients with medically refractory epilepsy often requires recording seizures using intracranial EEG. Quantitative measures derived from interictal intracranial EEG yield potentially appealing biomarkers to guide these surgical procedures; however, their utility is limited by the sparsity of electrode implantation as well as the normal confounds of spatiotemporally varying neural activity and connectivity. We propose that comparing intracranial EEG recordings to a normative atlas of intracranial EEG activity and connectivity can reliably map abnormal regions, identify targets for invasive treatment and increase our understanding of human epilepsy. Merging data from the Penn Epilepsy Center and a public database from the Montreal Neurological Institute, we aggregated interictal intracranial EEG retrospectively across 166 subjects comprising >5000 channels. For each channel, we calculated the normalized spectral power and coherence in each canonical frequency band. We constructed an intracranial EEG atlas by mapping the distribution of each feature across the brain and tested the atlas against data from novel patients by generating a z-score for each channel. We demonstrate that for seizure onset zones within the mesial temporal lobe, measures of connectivity abnormality provide greater distinguishing value than univariate measures of abnormal neural activity. We also find that patients with a longer diagnosis of epilepsy have greater abnormalities in connectivity. By integrating measures of both single-channel activity and inter-regional functional connectivity, we find a better accuracy in predicting the seizure onset zones versus normal brain (area under the curve = 0.77) compared with either group of features alone. We propose that aggregating normative intracranial EEG data across epilepsy centres into a normative atlas provides a rigorous, quantitative method to map epileptic networks and guide invasive therapy. We publicly share our data, infrastructure and methods, and propose an international framework for leveraging big data in surgical planning for refractory epilepsy.

Intracranial electroencephalographic biomarker predicts effective responsive neurostimulation for epilepsy prior to treatment

OBJECTIVE

Despite the overall success of responsive neurostimulation (RNS) therapy for drug-resistant focal epilepsy, clinical outcomes in individuals vary significantly and are hard to predict. Biomarkers that indicate the clinical efficacy of RNS-ideally before device implantation-are critically needed, but challenges include the intrinsic heterogeneity of the RNS patient population and variability in clinical management across epilepsy centers. The aim of this study is to use a multicenter dataset to evaluate a candidate biomarker from intracranial electroencephalographic (iEEG) recordings that predicts clinical outcome with subsequent RNS therapy.

METHODS

We assembled a federated dataset of iEEG recordings, collected prior to RNS implantation, from a retrospective cohort of 30 patients across three major epilepsy centers. Using ictal iEEG recordings, each center independently calculated network synchronizability, a candidate biomarker indicating the susceptibility of epileptic brain networks to RNS therapy.

RESULTS

Ictal measures of synchronizability in the high-γ band (95-105 Hz) significantly distinguish between good and poor RNS responders after at least 3 years of therapy under the current RNS therapy guidelines (area under the curve = .83). Additionally, ictal high-γ synchronizability is inversely associated with the degree of therapeutic response.

SIGNIFICANCE

This study provides a proof-of-concept roadmap for collaborative biomarker evaluation in federated data, where practical considerations impede full data sharing across centers. Our results suggest that network synchronizability can help predict therapeutic response to RNS therapy. With further validation, this biomarker could facilitate patient selection and help avert a costly, invasive intervention in patients who are unlikely to benefit.

OUP accepted manuscript

Objective

Materials and Methods

Results

Discussion and Conclusion

Electrocorticography and stereo EEG provide distinct measures of brain connectivity: implications for network models

Brain network models derived from graph theory have the potential to guide functional neurosurgery, and to improve rates of post-operative seizure freedom for patients with epilepsy. A barrier to applying these models clinically is that intracranial EEG electrode implantation strategies vary by centre, region and country, from cortical grid & strip electrodes (Electrocorticography), to purely stereotactic depth electrodes (Stereo EEG), to a mixture of both. To determine whether models derived from one type of study are broadly applicable to others, we investigate the differences in brain networks mapped by electrocorticography and stereo EEG in a cohort of patients who underwent surgery for temporal lobe epilepsy and achieved a favourable outcome. We show that networks derived from electrocorticography and stereo EEG define distinct relationships between resected and spared tissue, which may be driven by sampling bias of temporal depth electrodes in patients with predominantly cortical grids. We propose a method of correcting for the effect of internodal distance that is specific to electrode type and explore how additional methods for spatially correcting for sampling bias affect network models. Ultimately, we find that smaller surgical targets tend to have lower connectivity with respect to the surrounding network, challenging notions that abnormal connectivity in the epileptogenic zone is typically high. Our findings suggest that effectively applying computational models to localize epileptic networks requires accounting for the effects of spatial sampling, particularly when analysing both electrocorticography and stereo EEG recordings in the same cohort, and that future network studies of epilepsy surgery should also account for differences in focality between resection and ablation. We propose that these findings are broadly relevant to intracranial EEG network modelling in epilepsy and an important step in translating them clinically into patient care.

IRIS: A Modular Platform for Continuous Monitoring and Caretaker Notification in the Intensive Care Unit

OBJECTIVE

New approaches are needed to interpret large amounts of physiologic data continuously recorded in the ICU. We developed and prospectively validated a versatile platform (IRIS) for real-time ICU physiologic monitoring, clinical decision making, and caretaker notification.

METHODS

IRIS was implemented in the neurointensive care unit to stream multimodal time series data, including EEG, intracranial pressure (ICP), and brain tissue oxygenation (PbtO2), from ICU monitors to an analysis server. IRIS was applied for 364 patients undergoing continuous EEG, 26 patients undergoing burst suppression monitoring, and four patients undergoing intracranial pressure and brain tissue oxygen monitoring. Custom algorithms were used to identify periods of elevated ICP, compute burst suppression ratios (BSRs), and detect faulty or disconnected EEG electrodes. Hospital staff were notified of clinically relevant events using our secure API to route alerts through a password-protected smartphone application.

RESULTS

Sustained increases in ICP and concordant decreases in PbtO2 were reliably detected using user-defined thresholds and alert throttling. BSR trends computed by the platform correlated highly with manual neurologist markings (r2 0.633-0.781; p < 0.0001). The platform identified EEG electrodes with poor signal quality with 95% positive predictive value, and reduced latency of technician response by 93%.

CONCLUSION

This study validates a flexible real-time platform for monitoring and interpreting ICU data and notifying caretakers of actionable results, with potential to reduce the manual burden of continuous monitoring services on care providers.

SIGNIFICANCE

This work represents an important step toward facilitating translational medical data analytics to improve patient care and reduce health care costs.

Virtual resection predicts surgical outcome for drug-resistant epilepsy

Patients with drug-resistant epilepsy often require surgery to become seizure-free. While laser ablation and implantable stimulation devices have lowered the morbidity of these procedures, seizure-free rates have not dramatically improved, particularly for patients without focal lesions. This is in part because it is often unclear where to intervene in these cases. To address this clinical need, several research groups have published methods to map epileptic networks but applying them to improve patient care remains a challenge. In this study we advance clinical translation of these methods by: (i) presenting and sharing a robust pipeline to rigorously quantify the boundaries of the resection zone and determining which intracranial EEG electrodes lie within it; (ii) validating a brain network model on a retrospective cohort of 28 patients with drug-resistant epilepsy implanted with intracranial electrodes prior to surgical resection; and (iii) sharing all neuroimaging, annotated electrophysiology, and clinical metadata to facilitate future collaboration. Our network methods accurately forecast whether patients are likely to benefit from surgical intervention based on synchronizability of intracranial EEG (area under the receiver operating characteristic curve of 0.89) and provide novel information that traditional electrographic features do not. We further report that removing synchronizing brain regions is associated with improved clinical outcome, and postulate that sparing desynchronizing regions may further be beneficial. Our findings suggest that data-driven network-based methods can identify patients likely to benefit from resective or ablative therapy, and perhaps prevent invasive interventions in those unlikely to do so.

The sensitivity of network statistics to incomplete electrode sampling on intracranial EEG

Network neuroscience applied to epilepsy holds promise to map pathological networks, localize seizure generators, and inform targeted interventions to control seizures. However, incomplete sampling of the epileptic brain because of sparse placement of intracranial electrodes may affect model results. In this study, we evaluate the sensitivity of several published network measures to incomplete spatial sampling and propose an algorithm using network subsampling to determine confidence in model results. We retrospectively evaluated intracranial EEG data from 28 patients implanted with grid, strip, and depth electrodes during evaluation for epilepsy surgery. We recalculated global and local network metrics after randomly and systematically removing subsets of intracranial EEG electrode contacts. We found that sensitivity to incomplete sampling varied significantly across network metrics. This sensitivity was largely independent of whether seizure onset zone contacts were targeted or spared from removal. We present an algorithm using random subsampling to compute patient-specific confidence intervals for network localizations. Our findings highlight the difference in robustness between commonly used network metrics and provide tools to assess confidence in intracranial network localization. We present these techniques as an important step toward translating personalized network models of seizures into rigorous, quantitative approaches to invasive therapy.

High interictal connectivity within the resection zone is associated with favorable post-surgical outcomes in focal epilepsy patients

Patients with drug-resistant focal epilepsy are often candidates for invasive surgical therapies. In these patients, it is necessary to accurately localize seizure generators to ensure seizure freedom following intervention. While intracranial electroencephalography (iEEG) is the gold standard for mapping networks for surgery, this approach requires inducing and recording seizures, which may cause patient morbidity. The goal of this study is to evaluate the utility of mapping interictal (non-seizure) iEEG networks to identify targets for surgical treatment. We analyze interictal iEEG recordings and neuroimaging from 27 focal epilepsy patients treated via surgical resection. We generate interictal functional networks by calculating pairwise correlation of iEEG signals across different frequency bands. Using image coregistration and segmentation, we identify electrodes falling within surgically resected tissue (i.e. the resection zone), and compute node-level and edge-level synchrony in relation to the resection zone. We further associate these metrics with post-surgical outcomes. Greater overlap between resected electrodes and highly synchronous electrodes is associated with favorable post-surgical outcomes. Additionally, good-outcome patients have significantly higher connectivity localized within the resection zone compared to those with poorer postoperative seizure control. This finding persists following normalization by a spatially-constrained null model. This study suggests that spatially-informed interictal network synchrony measures can distinguish between good and poor post-surgical outcomes. By capturing clinically-relevant information during interictal periods, our method may ultimately reduce the need for prolonged invasive implants and provide insights into the pathophysiology of an epileptic brain. We discuss next steps for translating these findings into a prospectively useful clinical tool.

•

We analyze interictal iEEG recordings and neuroimaging from epilepsy patients

•

We find that high interictal strength selectivity is associated with better outcomes

•

This effect appears to be driven largely by connectivity within the resection zone

•

Interictal recordings can guide identification of seizure-generating networks

Field Distribution of Transcranial Static Magnetic Stimulation in Realistic Human Head Model

OBJECTIVE

The objective of this work was to characterize the magnetic field (B-field) that arises in a human brain model from the application of transcranial static magnetic field stimulation (tSMS).

MATERIALS AND METHODS

The spatial distribution of the B-field magnitude and gradient of a cylindrical, 5.08 cm × 2.54 cm NdFeB magnet were simulated in air and in a human head model using the finite element method and calibrated with measurements in air. The B-field was simulated for magnet placements over prefrontal, motor, sensory, and visual cortex targets. The impact of magnetic susceptibility of head tissues on the B-field was quantified.

RESULTS

Peak B-field magnitude and gradient respectively ranged from 179-245 mT and from 13.3-19.0 T/m across the cortical targets. B-field magnitude, focality, and gradient decreased with magnet-cortex distance. The variation in B-field strength and gradient across the anatomical targets largely arose from the magnet-cortex distance. Head magnetic susceptibilities had negligible impact on the B-field characteristics. The half-maximum focality of the tSMS B-field ranged from 7-12 cm³ .

SIGNIFICANCE

This is the first presentation and characterization of the three-dimensional (3D) spatial distribution of the B-field generated in a human brain model by tSMS. These data can provide quantitative dosing guidance for tSMS applications across various cortical targets and subjects. The finding that the B-field gradient is high near the magnet edges should be considered in studies where neural tissue is placed close to the magnet. The observation that susceptibility has negligible effects confirms assumptions in the literature.

Modeling transcranial electric stimulation in mouse: a high resolution finite element study

Mouse models are widely used in studies of various forms of transcranial electric stimulation (TES). However, there is limited knowledge of the electric field distribution induced by TES in mice, and computational models to estimate this distribution are lacking. This study examines the electric field and current density distribution in the mouse brain induced by TES. We created a high-resolution finite element mouse model incorporating ear clip electrodes commonly used in mouse TES to study, for example, electroconvulsive therapy (ECT). The electric field strength and current density induced by an ear clip electrode configuration were computed in the anatomically realistic, inhomogenous mouse model. The results show that the median electric field strength induced in the brain at 1 mA of stimulus current is 5.57 V/m, and the strongest field of 20.19 V/m was observed in the cerebellum. Therefore, to match the median electric field in human ECT at 800 mA current, the electrode current in mouse should be set to approximately 15 mA. However, the location of the strongest electric field in posterior brain regions in the mouse does not model well human ECT which targets more frontal regions. Therefore, the ear clip electrode configuration may not be a good model of human ECT. Using high-resolution realistic models for simulating TES in mice may guide the establishment of appropriate stimulation parameters for future in vivo studies.