Computational modeling of transcranial direct current stimulation (tDCS) in obesity: Impact of head fat and dose guidelines.

Recent studies show that acute neuromodulation of the prefrontal cortex with transcranial direct current stimulation (tDCS) can decrease food craving, attentional bias to food, and actual food intake. These data suggest potential clinical applications for tDCS in the field of obesity. However, optimal stimulation parameters in obese individuals are uncertain. One fundamental concern is whether a thick, low-conductivity layer of subcutaneous fat around the head can affect current density distribution and require dose adjustments during tDCS administration. The aim of this study was to investigate the role of head fat on the distribution of current during tDCS and evaluate whether dosing standards for tDCS developed for adult individuals in general are adequate for the obese population. We used MRI-derived high-resolution computational models that delineated fat layers in five human heads from subjects with body mass index (BMI) ranging from "normal-lean" to "super-obese" (20.9 to 53.5 kg/m(2)). Data derived from these simulations suggest that head fat influences tDCS current density across the brain, but its relative contribution is small when other components of head anatomy are added. Current density variability between subjects does not appear to have a direct and/or simple link to BMI. These results indicate that guidelines for the use of tDCS can be extrapolated to obese subjects without sacrificing efficacy and/or treatment safety; the recommended standard parameters can lead to the delivery of adequate current flow to induce neuromodulation of brain activity in the obese population.

Introduction

Obesity is a major public health concern worldwide. In the United States alone, 78 million adults and approximately 12.5 million children and adolescents were obese between 2009 and 2010 (Ogden et al., 2012). Research indicates that these numbers will continue to rise. The largest increase will be in severe obesity, with its accompanying surge in comorbid conditions and related healthcare costs (Finkelstein et al., 2012; Wang et al., 2011). The medical, social, and economic consequences of obesity have focused global attention on the condition and spawned numerous public health initiatives. Still, therapeutic options remain limited. New treatment strategies are required to halt the rise in obesity and limit future economic and societal costs.

A growing body of evidence, mostly from human neuroimaging studies, suggests that dysregulation in brain regions that process cognitive and reward aspects of food and eating behavior may be a key component of obesity (Alonso-Alonso and Pascual-Leone, 2007; Appelhans, 2009; Carnell et al., 2012; Dagher, 2012; Volkow et al., 2013; Zheng et al., 2009). Thus, modulating brain activity with neurotechnologies may open new therapeutic avenues. Compared to other neuromodulatory techniques, transcranial direct current stimulation (tDCS) offers significant advantages due to its relative safety, noninvasiveness, low-cost, and portability (Nitsche et al., 2008).

By delivering a weak direct current to the scalp via two electrodes – anode and cathode – tDCS can modulate the transmembrane potential of neurons, modify excitability, and induce plasticity changes. Over time, these can translate into clinical effects in diverse patient populations (Brunoni et al., 2012; Nitsche and Paulus, 2011; Nitsche et al., 2008). Preliminary, single-session data support a potential role for tDCS in the modulation of appetite and eating behavior in humans. In a randomized, sham-controlled, crossover study conducted in 23 subjects, Fregni et al. (2008) reported an acute decrease in food craving, as well as a reduction in snack consumption and eye gaze fixation to food following 20 min of tDCS applied over the prefrontal cortex. Similarly, a study in 19 subjects by Goldman et al. (2011) found that prefrontal tDCS caused a transient increase in the self-reported ability to resist food and a reduction in food cravings, particularly for sweet foods and carbohydrates; however there was no effect on ad libitum food intake. More recently, Montenegro et al. (2012) replicated the reduction of the desire to eat in 9 overweight subjects following a single session of prefrontal tDCS, and they also found an enhancement of the effect when tDCS was combined with a bout of aerobic exercise. Altogether, these three small studies provide initial rationale for the use of tDCS in clinical trials in the field of obesity.

To optimize stimulation parameters in obese subjects requires knowing the potential influence of head fat on current density distribution. It is well-established that head anatomy and variations in tissue layers, including fat (Shahid et al., 2011; Truong et al., 2012), affect how current density is distributed in the brain (Bikson et al., 2012a,b; Sadleir et al., 2010; Wagner et al., 2007). The identical tDCS montage applied to subjects with different head anatomy can produce varied intensity and pattern of current flow, which in turn may influence efficacy or safety. Even across anatomically typical adults, variation in peak cortical current density can vary > two-fold (Datta et al., 2012).

Therefore, the presence of a thickened layer of fat around the head in obese individuals could influence brain current flow and result in neuromodulation during tDCS administration. Investigating if and how to alter tDCS dose to accommodate variations in BMI is timely. Interest in the use of this technology in obese subjects is growing, for both the modulation of craving-related processes, and more broadly, for neuropsychiatric treatment of patients who often have obesity as a comorbidity. The purpose of this study was to systematically examine the role of head fat on the distribution of current during tDCS using MRI-derived high-resolution computational models, and to evaluate whether tDCS dosing standards developed for adults in general are adequate for the obese population.

Methods

Subjects

To determine the effect of head fat on current density distribution during tDCS, we created models from MRI images of five human subjects categorized according to BMI, from normal (18.5–24.9 kg/m2) to super obese (> 50 kg/m2). Subjects were a 35-year-old female with a BMI of 53.5 kg/m2 (S1), a 47-year-old female with a BMI of 43.4 kg/m2 (S2), a 22-year-old female with a BMI of 38.3 kg/m2 (S3), and a 25-year-old female with a BMI of 20.9 kg/m2 (S4). We also included a 36-year-old male subject with a BMI of 25.1 kg/m2 (S#) who participated in prior tDCS computational modeling studies (Bikson et al., 2010; Datta et al., 2012). Subjects S1–S4 underwent an MRI as part of research studies related to eating behavior and obesity. Study procedures were approved by the Institutional Review Board of Beth Israel Deaconess Medical Center.

MRI data collection and segmentation

We performed high-resolution T1-weighted magnetization prepared rapid gradient echo (MPRAGE) MRI scans at the Center for Biomedical Imaging, Boston University School of Medicine, using a 3-T Philips Achieva scanner (Philips Medical Systems, Best, The Netherlands) equipped with a Synergy-L Sensitivity Encoding (SENSE) head coil. Acquisition parameters were: TE = 3.2 ms; TR = 6.92 ms; flip angle = 8°; FOV = 256 mm; resolution = 256 × 256; slice thickness = 1.2 mm; no gap; and voxel size of 1 × 1 × 1.2 mm. The scans were segmented into 7 tissues: air, skin, fat, skull, cerebral spinal fluid (CSF), gray matter, and white matter.

Automated segmentation algorithms from Statistical Parametric Mapping (SPM8, Wellcome Trust Centre for Neuroimaging, London, UK) were used in conjunction with updated tissue probability maps (Rorden et al., 2012) to generate an initial segmentation of air, skin, skull, CSF, gray matter, and white matter. Additional post-processing algorithms smoothed artifacts and corrected for discontinuities (Huang et al., 2012). We added fat segmentation through a threshold flood fill of skin (fat has high signal intensity on T1-weighted MRI images) and manually corrected lingering errors in continuity and detail in all tissues with ScanIP 4.2 (Simpleware Ltd, Exeter, UK). Tissue continuity was verified after the results of the automatic segmentation, visualizing the data extensively. Further manual adjustments were performed to guarantee continuity and improve the accuracy of the segmentation so that the tissue masks matched closely the real anatomy of the individual. All these procedures were carried out by two team members using the visualization options and tools provided in the ScanIP 4.2 software.

Two models (S# and S4) were artificially “fattened” by dilating the segmentation of fat. Fat was merged with the outer surface of skin, and then dilated isometrically up to 10 mm. This dilation caused fat to overtake the original skin surface by expanding outward. To recover the skin we made a duplicate of this merged fat and skin segmentation mask and dilated it an additional 3 mm to form the new skin surface. As a result of this transformation skin and fat were still distinct segmentation masks; fat became thicker, and the thickness of skin was fixed at 3 mm. No tissues other than skin and fat were altered in these models.

We measured the thicknesses of skin, fat, bone, and CSF for each model from the segmentation data. Measurements were performed over both motor strips (C3 and C4), 5 times each, and averaged. Tissue thickness was measured in three-dimensions by sampling a patch of tissue with a bounding box. The volume of the tissue within the bounding box was determined by summating the segmentation voxels (1 mm3). This volume was then projected from the area of the diagonal plane within the bounding box so that thickness (length ‘L’) equals tissue volume divided by plane area (L = V / A). The measurements were recorded from 10:20 positions C3 and C4 where the plane was tangential to the scalp in a coronal slice.

Modeling of tDCS

Stimulation electrodes, sponge pads, and gels were modeled in SolidWorks (Dassault Systèmes Corp., Waltham, MA) and imported into ScanIP for meshing. Three montages were modeled: 5 × 7 cm pads with anode over the motor strip (C3) and cathode over the contralateral supra-orbital (M1-SO); 5 × 7 cm pads with anode over the inferior frontal gyrus (F8) and cathode over the contralateral supra-orbital (IFG-SO); and a high-definition (HD) electrode ring configuration designed for anodal stimulation over the motor strip (4 × 1 over C3; 5 cm radius from center electrode to outer electrodes). An adaptive tetrahedral meshing algorithm was used in ScanIP to generate meshes between 6 × 106 and 14 × 106 quadratic elements.

Finite element method (FEM) models were created in COMSOL multiphysics 3.5a (COMSOL, Inc., Burlington, MA) using the aforementioned meshes. Models were created using electrostatic volume conductor physics with material conductivities defined as follows: (in S/m): air, 1 × 10− 15; skin, 0.465; fat, 0.025; skull, 0.01; CSF, 1.65; gray matter, 0.276; white matter, 0.126; electrode, 5.99 × 107; saline-soaked sponge, 1.4; and conductive gel, 0.3. These conductivity values used were the same as previously published modeling work drawing on data from a combination of in vivo and in vitro measurements (Datta et al., 2011; Gabriel et al., 1996). We applied boundary conditions to simulate direct current stimulation. Internal boundaries between tissues were assigned the continuity condition (n ∗ (J1 − J2) = 0), and the Laplace equation (∇ ∗ (σ∇V) = 0) was solved. The resulting cortical electric field was interpreted as a correlate for modulation (Bikson et al., 2004; Tranchina and Nicholson, 1986). The surfaces of the cathodes were grounded (V = 0), while the surfaces of the anodes had a current density of 1 A/m2. All other exterior surfaces were electrically insulated. Peak electric field data are provided in the text as absolute values or as mean (μ) and standard deviation (σ).

To examine the role of head fat and other tissue thickness in current distribution, we multiplied conductivity values by the average tissue thickness measured under C3 and C4 for each subject model. Also, a “weighted thickness” measure was calculated by summing each tissue thickness weighted by the associated conductivity. Linear associations between peak current intensity, tissue thickness and BMI were evaluated using Spearman's rank correlation coefficient as the data did not follow a normal distribution.

Results

Optimized segmentation, including fat delineation

Fig. 1 shows the segmentation of all head tissue layers for each of the five subjects. We were able to delineate essential tissues, including the fat layer, with great detail in all cases. In all the tissue layers, including fat, we observed a high degree of interindividual variability in head anatomy.

Fig. 1

Segmentation of five subjects with varying BMI (S1, S2, S3, S4, S#), six tissue compartment models [skin, fat, skull, cerebral spinal fluid (CSF), gray matter and white matter]. High-resolution MRI scans were segmented using a combination of automated and manual techniques. Specific anatomical considerations, such as continuity of CSF, were verified or corrected. Images are shown on the same scale.

Current distribution in three tDCS montages

We tested two standard (5 × 7 cm pads) tDCS montages and a HD-tDCS montage with the ring (4 × 1) electrode configuration, using a total of 1 mA of current injection in each case. Models predicted substantial interindividual variability in current peak values and distribution (Fig. 2). 4 × 1-ring HD-tDCS resulted in more focal and consistent distribution of electric currents. From lowest to highest peak amplitude (sensitivity), the individuals ranked S3, S2, S#, S1, and S4; this ranking was the same in each of the montages.

Fig. 2

Resulting peak electric field magnitude simulated in three montages (M1–SO, 4 × 1 HD-tDCS over M1, IFG-SO) across subjects. Variations in intensity occur across individuals, but these individual variations are consistent in ranking across montages (S3 < S2 < S# < S1 < S4). For each of the models, scale is adjusted to the current density peak, which is also provided as value for reference.

Our results corroborate with previous modeling studies using single heads which shows diffuse and clustered brain current flow with conventional tDCS pad montages, while 4 × 1 HD-tDCS produces focal current flow with moderately lower peak electric field for the same total applied current (Datta et al., 2009; Dmochowski et al., 2011). Comparing across subjects, we observed differences in both peak brain electric field intensity and individual variation between montages. 4 × 1 HD-tDCS was the least intense on average (μ = 0.190 V/m) with the most variation in peak electric field (σ = 0.094 V/m). However, spatial targeting was controlled across subjects to/within the ring. Conventional tDCS M1-SO and IFG-SO montages were comparable in intensity with relatively less variation in peak electric field across subjects (μ = 0.317 and 0.330 V/m; σ = 0.041 and 0.039 V/m); the location of the peak and indeed overall distribution varied across subjects.

Given the differences across montages, it is notable that ranking of sensitivity, at least for this sample set, remained fixed across the montages evaluated. Simultaneously, differences across subjects were as significant as across montages. Implications of these findings for rational tDCS design are considered in the discussion.

Role of head fat and other tissue thickness in current distribution

Table 1 shows individual data for specific tissue thicknesses as well as total thickness weighed by conductivity. We observed a positive trend between BMI and head fat thickness (Spearman's rho = 0.8; p = 0.107), but did not predict a consistent monotonic (linear) association between current intensity and subject BMI, or thickness of either skin, fat, skull, or CSF. For example, though the highest predicted peak was in the individual with the lowest BMI (S4; 20.9 kg/m2), the subject with the intermediate BMI in the sample (S3; BMI: 38.3 kg/m2) corresponded to the lowest overall predicted intensity. Ranking order of sensitivity to tDCS was thus not consistent with the increasing/decreasing ranking order for either BMI, fat or other tissue thickness, or total tissue thickness (distance of the electrode to brain).

Table 1

Quantifying individual differences. BMI and thickness of tissues surrounding the brain were measured at EEG 10–20 positions C3 and C4. Total thickness and total thickness weighted by conductivity are also listed. Images are shown on the same scale.

	S1	S2	S3	S4	S#
BMI	53.5	43.4	38.3	20.9	25.1
Skin (mm)	3.79	5.38	5.23	3.10	5.79
Fat (mm)	4.98	5.18	2.30	1.83	2.43
Bone (mm)	3.17	4.90	5.68	4.79	4.75
CSF (mm)	2.71	2.45	2.71	2.26	2.30
Total (mm)	14.65	17.91	15.92	11.97	15.26
Weighted conductance (mS)	6.39	6.72	7.02	5.26	6.60

In post-hoc analysis, in this sample, we observed that the rank of sensitivity to tDCS was consistent with the decreasing “weighted thickness”—which was calculated by summing each tissue thickness weighted by the associated conductivity (Spearman's rho = − 1; p = 0.044). This empirical “weighted thickness” metric is more influenced by the relatively conductive layers (skin and CSF) that may reflect the role of shunting through these tissues in the decreasing peak cortical electric field.

Effect of intrasubject fat layer dilation

In the absence of an evident relationship between BMI and brain current flow across subjects (where other tissues were also different), we evaluated the influence of increasing fat thickness within an individual (with all other factors being equal). In the fat dilated head models, we observed a drop in peak electric field in an extreme, but not physiologically typical, scenario (fat thickness > 10–15 mm) (Hayman et al., 2003) (Fig. 3). Increasing the layer of fat at a physiologically observed range (a few millimeters (Hayman et al., 2003)) did not have a significant effect or even a consistent effect on the direction of change. In the discussion we address that insensitivity to increasing fat thickness does not equate to a negligible role for baseline fat.

Fig. 3

Influence of fat thickness in isolation. Fat was dilated isometrically with 3 mm of skin cover; other tissues were unchanged. A moderate increase in the thickness of fat caused little change in peak electric field. There was a slight increase (5.7%) in S# and a slight decrease (8.9%) in S4. Increasing the thickness of fat beyond that physiologically observed led to noticeable decreases in intensity in both S# (15.6%) and S4 (25.7%). Varying the thickness of tissues surrounding the brain not only changes the overall conductance, but also changes the orientation of the electrodes with respect to the brain.

Skin current density across subjects and montages

As a proxy for skin tolerability, we analyzed the current density of the skin surface (Fig. 4). Variations in current density magnitude were minimal between subjects. The spatial distribution of current density resembled that of previous studies resulting in hot spots along the edges of the contact boundary (Minhas et al., 2010) although given the complexity of electrode design and skin interface (Kronberg and Bikson, 2012), we propose that these models are useful primarily for qualitative comparisons across heads.

Fig. 4

Skin current density across subjects and montages. The largest change in current density magnitude was between the montages utilizing HD electrodes versus conventional pads; interindividual differences are relatively minor. Peak current density was calculated per mA of stimulation.

Discussion

Variations in brain current flow with increasing BMI

In this study, we used high-resolution MRI-derive computational models to predict the effect of BMI and head fat on current density distribution. We used computational models to systematically address this problem in three selected tDCS montages that were simulated in five human subjects with different BMIs. We found that current density variability between subjects does not appear to have a direct and/or simple link to BMI. For example, we observed that peak amplitudes in an extreme case of obesity (BMI > 50 kg/m2) were comparable to those found in non-obese cases. Further, simulated dilation of the fat layer revealed a within-subject significant effect only at supraphysiological values of fat thickness.

When combined with previous modeling studies, our results suggest that head fat contributes to current density distribution in conjunction with other anatomical differences. Ultimately, the variation among individuals is likely the result of a multitude of factors, not just BMI. According to our data, differences in head fat thickness contribute an extra ~ 10% variability in peak cortical current density in addition to the previously reported 1.5-to-3-fold variability that exists across normal (non-obese) individuals (Datta et al., 2012).

Effect of fat on tDCS current distribution and implications for low BMI

Our findings show that BMI does not, in itself, significantly predict brain current flow intensity, nor do physiological increases in individual BMI profoundly influence current flow. Yet these findings do not diminish the validity of studies indicating that fat influences current flow and that the omission of fat in computational models (e.g., representation as skin) reduces precision (Shahid et al., 2011). In the first case, the hypothetical removal of fat will influence current flow. In the second case, failure to implement fat in computational models may change predicted brain current by up to 60% (Truong et al., 2012). Our results reinforce the utility of individual, MRI-derived computational models, and their value in guiding and supporting the development of new clinical applications of tDCS.

The largest peak amplitudes in all montages were found in the lowest BMI (BMI 20.9 kg/m2), which also corresponded to the smallest head size. This observation, in combination with the impact of model omission, suggests that current distribution may alter significantly in individuals with lower-than-normal BMI. We did not sample the underweight BMI spectrum (< 18.5 kg/m2). However, future studies are needed to address and clarify this potential issue and its clinical implications for the future use of tDCS in low BMI individuals, including those with anorexia nervosa or cachexia (Hecht, 2010).

Clinical efficacy and safety considerations

In contrast to noninvasive approaches such as transcranial magnetic stimulation (TMS) and electroconvulsive therapy (ECT), tDCS dose is currently not typically adjusted across subjects and the latest clinical guidelines do not recommend individual titration (Brunoni et al., 2012; Nitsche and Paulus, 2011; Nitsche et al., 2008). Nonetheless, just as other neuromodulation modalities benefit from individual dose adjustment, it is increasingly recognized that variation in head anatomy may alter current delivery to the brain during tDCS (Datta et al., 2012) and that this may account for variability in cognitive and behavioral outcomes. Here it was considered if at one extreme, the increased head fat (a relatively resistive tissue) associated with high-BMI/obesity, will alter brain current flow in a meaningful manner. Our modeling data suggest that compared to variations seen in healthy lean subjects (Datta et al., 2012), head fat has a minor influence on current density distribution, meaning its relative contribution is small when other sources of variability related to head anatomy are added. Therefore, no special considerations regarding tDCS dose and safety may be needed for use in clinical trials involving overweight or obese individuals.

Whether individualized dose optimization is warranted on the population as a whole, including obese individuals remains an open question that is beyond the scope of this report. Generally, efficacy is normalized by considering electric field intensity in the target and/or distribution (focality) across the heads, computational models allow for these consideration and even automatic optimization (Dmochowski et al., 2011).

For safety and tolerability, effects on the brain and at the skin can be considered independently (Bikson et al., 2009). In regards to brain electric fields, there is currently no evidence that conventionally used protocols approach hazardous levels (Bikson et al., 2009; Liebetanz et al., 2009) that may allow a safety factor for normal variation across subjects. Nonetheless, as with any new investigation, caution should be used when applying this technique in clinical applications, including in the field of obesity and eating disorders. For example, these models, which predict only current flow, do not consider potential differences in neurophysiological changes and/or sensitivity for the same brain electric field.

Skin related adverse effects (tolerability, and pain/sensation) are presumed linked to current density at the skin. Though skin current density predictions are limited by model sensitivity to idiosyncratic details of the electrode/skin surface (Kronberg and Bikson, 2012), our results (consistent with (Truong et al., 2012)) suggest that gross changes in the anatomy in obesity do not influence current density at the skin. However, as the prevalence of skin conditions tends to be higher in obese than in normal-BMI individuals (Scheinfeld, 2004), we recommend careful interviewing of subjects and skin inspection as suggested in current tDCS guidelines (Nitsche and Paulus, 2011; Nitsche et al., 2008).

These results are in line with recent reports of an acute decrease in self-reported measures of food craving following one session of tDCS over the dorsolateral prefrontal cortex in overweight/obese subjects (BMI 25.2 to 43.5 kg/m2; Montenegro et al. (2012); BMI > 30 kg/m2 in 31.6% of the sample; Goldman et al. (2011)). Neither of these investigations mentioned adverse effects related to tDCS administration using standard electrode sponges at 2 mA intensity for 20 min. The combination of modeling and experimental evidence suggests that the current guidelines are both safe and sufficient for neuromodulation of brain activity across the normal-to-obese BMI spectrum.

Computational models in tDCS design and special considerations at extremes of BMI

This study makes predictions based on computational models with precision limited by the accuracy of segmentation (Fig. 1) and tissue conductivity assignments. Other permutations and refinements, such as the addition of further tissue masks and anisotropy, only have value in informing clinical guidelines if 1) extra precision is added rationally rather than for complexity; and 2) relative changes in current flow across models have a significant effect on clinical dosing decisions (Bikson and Datta, 2012; Bikson et al., 2012b). However, if as supported by prior and the present study, individual adjustment is warranted across the general population, then added complexity must be balanced against costs associated with higher volume modeling.

Conclusions

In sum, evidence indicates that current guidelines for the administration of tDCS in the general population can be extended to those who have obesity. High-resolution computational models that include head fat provide individualized prediction of tDCS current density, and can accurately guide and support tDCS protocols in emerging clinical applications.