The effectiveness of transcranial magnetic stimulation (TMS) depends highly on the coil orientation relative to the subject's head. This implies that the direction of the induced electric field has a ...large effect on the efficiency of TMS. To improve future protocols, knowledge about the relationship between the coil orientation and the direction of the induced electric field on the one hand, and the head and brain anatomy on the other hand, seems crucial. Therefore, the induced electric field in the cortex as a function of the coil orientation has been examined in this study.
The effect of changing the coil orientation on the induced electric field was evaluated for fourteen cortical targets. We used a finite element model to calculate the induced electric fields for thirty-six coil orientations (10 degrees resolution) per target location. The effects on the electric field due to coil rotation, in combination with target site anatomy, have been quantified.
The results confirm that the electric field perpendicular to the anterior sulcal wall of the central sulcus is highly susceptible to coil orientation changes and has to be maximized for an optimal stimulation effect of the motor cortex. In order to obtain maximum stimulation effect in areas other than the motor cortex, the electric field perpendicular to the cortical surface in those areas has to be maximized as well. Small orientation changes (10 degrees) do not alter the induced electric field drastically.
The results suggest that for all cortical targets, maximizing the strength of the electric field perpendicular to the targeted cortical surface area (and inward directed) optimizes the effect of TMS. Orienting the TMS coil based on anatomical information (anatomical magnetic resonance imaging data) about the targeted brain area can improve future results. The standard coil orientations, used in cognitive and clinical neuroscience, induce (near) optimal electric fields in the subject-specific head model in most cases.
The inconsistent response to transcranial electric stimulation in the stroke population is attributed to, among other factors, unknown effects of stroke lesion conductivity on stimulation strength at ...the targeted brain areas. Volume conduction models are promising tools to determine optimal stimulation settings. However, stroke lesion conductivity is often not considered in these models as a source of inter-subject variability. The goal of this study is to propose a method that combines MRI, EEG, and transcranial stimulation to estimate the conductivity of cortical stroke lesions experimentally. In this simulation study, lesion conductivity was estimated from scalp potentials during transcranial electric stimulation in 12 chronic stroke patients. To do so, first, we determined the stimulation configuration where scalp potentials are maximally affected by the lesion. Then, we calculated scalp potentials in a model with a fixed lesion conductivity and a model with a randomly assigned conductivity. To estimate the lesion conductivity, we minimized the error between the two models by varying the conductivity in the second model. Finally, to reflect realistic experimental conditions, we test the effect rotation of measurement electrode orientation and the effect of the number of electrodes used. We found that the algorithm converged to the correct lesion conductivity value when noise on the electrode positions was absent for all lesions. Conductivity estimation error was below 5% with realistic electrode coregistration errors of 0.1° for lesions larger than 50 ml. Higher lesion conductivities and lesion volumes were associated with smaller estimation errors. In conclusion, this method can experimentally estimate stroke lesion conductivity, improving the accuracy of volume conductor models of stroke patients and potentially leading to more effective transcranial electric stimulation configurations for this population.
Many human cortical regions are targeted with transcranial magnetic stimulation (TMS). The stimulus intensity used for a certain region is generally based on the motor threshold stimulation intensity ...determined over the motor cortex (M1). However, it is well known that differences exist in coil-target distance and target site anatomy between cortical regions. These differences may well make the stimulation intensity derived from M1 sub-optimal for other regions. Our goal was to determine in what way the induced electric fields differ between cortical target regions. We used finite element method modeling to calculate the induced electric field for multiple target sites in a realistic head model. The effects on the electric field due to coil-target distance and target site anatomy have been quantified. The results show that a correction based on the distance alone does not correctly adjust the induced electric field for regions other than M1. In addition, a correction based solely on the TMS-induced electric field (primary field) does not suffice. A precise adjustment should include coil-target distance, the secondary field caused by charge accumulation at conductivity discontinuities and the direction of the field relative to the local cerebrospinal fluid–grey matter boundary.
Transcranial direct current stimulation (tDCS) is a noninvasive brain stimulation technique able to induce long-lasting changes in cortical excitability that can benefit cognitive functioning and ...clinical treatment. In order to both better understand the mechanisms behind tDCS and possibly improve the technique, finite element models are used to simulate tDCS of the human brain. With the detailed anisotropic head model presented in this study, we provide accurate predictions of tDCS in the human brain for six of the practically most-used setups in clinical and cognitive research, targeting the primary motor cortex, dorsolateral prefrontal cortex, inferior frontal gyrus, occipital cortex, and cerebellum. We present the resulting electric field strengths in the complete brain and introduce new methods to evaluate the effectivity in the target area specifically, where we have analyzed both the strength and direction of the field. For all cerebral targets studied, the currently accepted configurations produced sub-optimal field strengths. The configuration for cerebellum stimulation produced relatively high field strengths in its target area, but it needs higher input currents than cerebral stimulation does. This study suggests that improvements in the effects of transcranial direct current stimulation are achievable.
Electrical activity of the myocardium is recorded with the 12-lead ECG. ECG simulations can improve our understanding of the relation between abnormal ventricular activation in diseased myocardium ...and body surface potentials (BSP). However, in equivalent dipole layer (EDL)-based ECG simulations, the presence of diseased myocardium breaks the equivalence of the dipole layer. To simulate diseased myocardium, patches with altered electrophysiological characteristics were incorporated within the model. The relation between diseased myocardium and corresponding BSP was investigated in a simulation study.
Activation sequences in normal and diseased myocardium were simulated and corresponding 64-lead BSP were computed in four models with distinct patch locations. QRS-complexes were compared using correlation coefficient (CC). The effect of different types of patch activation was assessed. Of one patient, simulated electrograms were compared to electrograms recorded during invasive electro-anatomical mapping.
Hundred-fifty-three abnormal activation sequences were simulated. Median QRS-CC of delayed versus dyssynchronous were significantly different (1.00 vs. 0.97,
< 0.001). Depending on the location of the patch, BSP leads were affected differently. Within diseased regions, fragmentation, low bipolar voltages and late potentials were observed in both recorded and simulated electrograms.
A novel method to simulate cardiomyopathy in EDL-based ECG simulations was established and evaluated. The new patch-based approach created a realistic relation between ECG waveforms and underlying activation sequences. Findings in the simulated cases were in agreement with clinical observations. With this method, our understanding of disease progression in cardiomyopathies may be further improved and used in advanced inverse ECG procedures.
Background:
The detection and localization of electrophysiological substrates currently involve invasive cardiac mapping. Electrocardiographic imaging (ECGI) using the equivalent dipole layer (EDL) ...method allows the noninvasive estimation of endocardial and epicardial activation and repolarization times (AT and RT), but the RT validation is limited to
in silico
studies. We aimed to assess the temporal and spatial accuracy of the EDL method in reconstructing the RTs from the surface ECG under physiological circumstances and situations with artificially induced increased repolarization heterogeneity.
Methods:
In four Langendorff-perfused pig hearts, we simultaneously recorded unipolar electrograms from plunge needles and pseudo-ECGs from a volume-conducting container equipped with 61 electrodes. The RTs were computed from the ECGs during atrial and ventricular pacing and compared with those measured from the local unipolar electrograms. Regional RT prolongation (cooling) or shortening (pinacidil) was achieved by selective perfusion of the left anterior descending artery (LAD) region.
Results:
The differences between the computed and measured RTs were 19.0 ± 17.8 and 18.6 ± 13.7 ms for atrial and ventricular paced beats, respectively. The region of artificially delayed or shortened repolarization was correctly identified, with minimum/maximum RT roughly in the center of the region in three hearts. In one heart, the reconstructed region was shifted by ~2.5 cm. The total absolute difference between the measured and calculated RTs for all analyzed patterns in selectively perfused hearts (
n
= 5) was 39.6 ± 27.1 ms.
Conclusion:
The noninvasive ECG repolarization imaging using the EDL method of atrial and ventricular paced beats allows adequate quantitative reconstruction of regions of altered repolarization.
Promising results have been reported in noninvasive estimation of cardiac activation times (AT) using the equivalent dipole layer (EDL) source model in combination with the boundary element method ...(BEM). However, the assumption of equal anisotropy ratios in the heart that underlies the EDL model does not reflect reality. In the present study, we quantify the errors of the nonlinear AT imaging based on the EDL approximation. Nine different excitation patterns (sinus rhythm and eight ectopic beats) were simulated with the monodomain model. Based on the bidomain theory, the body surface potential maps (BSPMs) were calculated for a realistic finite element volume conductor with an anisotropic heart model. For the forward calculations, three cases of bidomain conductivity tensors in the heart were considered: isotropic, equal, and unequal anisotropy ratios in the intra- and extracellular spaces. In all inverse reconstructions, the EDL model with BEM was employed: AT were estimated by solving the nonlinear optimization problem with the initial guess provided by the fastest route algorithm. Expectedly, the case of unequal anisotropy ratios resulted in larger localization errors for almost all considered activation patterns. For the sinus rhythm, all sites of early activation were correctly estimated with an optimal regularization parameter being used. For the ectopic beats, all but one foci were correctly classified to have either endo- or epicardial origin with an average localization error of 20.4 mm for unequal anisotropy ratio. The obtained results confirm validation studies and suggest that cardiac anisotropy might be neglected in clinical applications of the considered EDL-based inverse procedure.
This study presents a novel non-invasive equivalent dipole layer (EDL) based inverse electrocardiography (
i
ECG) technique which estimates both endocardial and epicardial ventricular activation ...sequences. We aimed to quantitatively compare our
i
ECG approach with invasive electro-anatomical mapping (EAM) during sinus rhythm with the objective of enabling functional substrate imaging and sudden cardiac death risk stratification in patients with cardiomyopathy. Thirteen patients (77% males, 48 ± 20 years old) referred for endocardial and epicardial EAM underwent 67-electrode body surface potential mapping and CT imaging. The EDL-based
i
ECG approach was improved by mimicking the effects of the His-Purkinje system on ventricular activation. EAM local activation timing (LAT) maps were compared with
i
ECG-LAT maps using absolute differences and Pearson’s correlation coefficient, reported as mean ± standard deviation 95% confidence interval. The correlation coefficient between
i
ECG-LAT maps and EAM was 0.54 ± 0.19 0.49–0.59 for epicardial activation, 0.50 ± 0.27 0.41–0.58 for right ventricular endocardial activation and 0.44 ± 0.29 0.32–0.56 for left ventricular endocardial activation. The absolute difference in timing between
i
ECG maps and EAM was 17.4 ± 7.2 ms for epicardial maps, 19.5 ± 7.7 ms for right ventricular endocardial maps, 27.9 ± 8.7 ms for left ventricular endocardial maps. The absolute distance between right ventricular endocardial breakthrough sites was 30 ± 16 mm and 31 ± 17 mm for the left ventricle. The absolute distance for latest epicardial activation was median 12.8 IQR: 2.9–29.3 mm. This first in-human quantitative comparison of
i
ECG and invasive LAT-maps on both the endocardial and epicardial surface during sinus rhythm showed improved agreement, although with considerable absolute difference and moderate correlation coefficient. Non-invasive
i
ECG requires further refinements to facilitate clinical implementation and risk stratification.
•Conventional anodal M1 tDCS leads to variable electric fields in stroke patients.•Optimizing electrode locations increases stimulation strength in stroke patients.•Individual brain anatomy and ...function require consideration in tDCS for stroke patients.•Lack of individualization may partially explain mixed tDCS effects in clinical RCTs.
Transcranial direct current stimulation (tDCS) is a promising tool to improve and speed up motor rehabilitation after stroke, but inconsistent clinical effects refrain tDCS from clinical implementation. Therefore, this study aimed to assess the need for individualized tDCS configurations in stroke, considering interindividual variability in brain anatomy and motor function representation.
We simulated tDCS in individualized MRI-based finite element head models of 21 chronic stroke subjects and 10 healthy age-matched controls. An anatomy-based stimulation target, i.e. the motor hand knob, was identified with MRI, whereas a motor function-based stimulation target was identified with EEG. For each subject, we simulated conventional anodal tDCS electrode configurations and optimized electrode configurations to maximize stimulation strength within the anatomical and functional target. The normal component of the electric field was extracted and compared between subjects with stroke and healthy, age-matched controls, for both targets, during conventional and optimized tDCS.
Electrical field strength was significantly lower, more variable and more frequently in opposite polarity for subjects with stroke compared to healthy age-matched subjects, both for the anatomical and functional target with conventional, i.e. non-individualized, electrode configurations. Optimized, i.e. individualized, electrode configurations increased the electrical field strength in the anatomical and functional target for subjects with stroke but did not reach the same levels as in healthy subjects.
Considering individual brain structure and motor function is crucial for applying tDCS in subjects with stroke. Lack of individualized tDCS configurations in subjects with stroke results in lower electric fields in stimulation targets, which may partially explain the inconsistent clinical effects of tDCS in stroke trials.