Abstract
Transcranial direct current stimulation (tDCS) is a variant of noninvasive neuromodulation, which promises treatment for brain diseases like major depressive disorder. In experiments, long-lasting aftereffects were observed, suggesting that persistent plastic changes are induced. The mechanism underlying the emergence of lasting aftereffects, however, remains elusive. Here we propose a model, which assumes that tDCS triggers a homeostatic response of the network involving growth and decay of synapses. The cortical tissue exposed to tDCS is conceived as a recurrent network of excitatory and inhibitory neurons, with synapses subject to homeostatically regulated structural plasticity. We systematically tested various aspects of stimulation, including electrode size and montage, as well as stimulation intensity and duration. Our results suggest that transcranial stimulation perturbs the homeostatic equilibrium and leads to a pronounced growth response of the network. The stimulated population eventually eliminates excitatory synapses with the unstimulated population, and new synapses among stimulated neurons are grown to form a cell assembly. Strong focal stimulation tends to enhance the connectivity within new cell assemblies, and repetitive stimulation with well-chosen duty cycles can increase the impact of stimulation even further. One long-term goal of our work is to help in optimizing the use of tDCS in clinical applications.
Author Summary
Noninvasive brain stimulation techniques like tDCS have the potential to directly interfere with neural activity, but may also trigger activity-dependent plasticity. We propose a model to study the mechanism of tDCS and persistent aftereffects that may be induced as a consequence of homeostatic structural plasticity. Based on the idea that tDCS perturbs the ongoing activity of neurons, our model predicts that the stimulation also triggers a rearrangement of synapses among stimulated and unstimulated neurons, eventually leading to network remodeling and cell assembly formation. Focal and strong stimulation leads to stronger cell assemblies, and so does repetitive stimulation with optimized stimulation protocols. This is the first original work studying possible long-lasting aftereffects of transcranial stimulation at the mesoscopic neuronal network level using a computational model.
INTRODUCTION
Transcranial direct current stimulation (tDCS) is a noninvasive brain stimulation technique, where a weak constant current (1−2 mA) is applied to the brain via large electrodes attached to the scalp (Edwards et al., 2013). It induces weak electric fields that are typically not sufficient to trigger action potentials directly, but can polarize the membrane of neurons by fractions of millivolts (Joucla & Yvert, 2009), depending on the orientation of the electric field vector relative to the somato-dendritic axis of the neuron (Wiethoff, Hamada, & Rothwell, 2014; Gluckman et al., 1996; Radman, Ramos, Brumberg, & Bikson, 2009). This membrane potential deflection can influence spike timing and firing rates of neurons that are part of an active network (Bikson, Radman, & Datta, 2006; Vöröslakos et al., 2018). Similar to other methods of neuromodulation, tDCS is claimed to have a potential for alleviating symptoms of certain brain diseases, such as major depressive disorder (Nitsche, Boggio, Fregni, & Pascual-Leone, 2009; Loo et al., 2012) or chronic pain (Garcia-Larrea, 2016; Ngernyam et al., 2015).
Although there is a record of promising applications of tDCS, both positive and negative outcomes have been reported in the literature (Horvath, Forte, & Carter, 2015). Typical issues are due to insufficient sensitivity of measurements, or large intersubject and intrasubject variability (Wiethoff et al., 2014). Positive evidence includes immediate changes of neural activity caused by tDCS, observed both in humans and in rodents. Positron emission tomography (PET) in humans revealed that tDCS can influence the activity of neurons in different brain regions (Lang et al., 2005), but the most affected region varies with electrode montage (Kuo et al., 2013), skull thickness (Opitz, Paulus, Will, Antunes, & Thielscher, 2015), individual geometry of cortex (Opitz et al., 2015), preexisting lesions (Minjoli et al., 2017), and other aspects. Systematic transcutaneous current stimulation experiments in rats (Vöröslakos et al., 2018) could establish quantitative relations between the externally applied current, the induced electric field, the associated membrane potential deflection, and the resulting firing rate change.
In addition to the instant impact on activity during stimulation, a sustained modulation of neural activity was also observed in humans after stimulation was turned off. Lasting aftereffects of tDCS, measured as motor evoked potentials (MEP) triggered by transcranial magnetic stimulation (TMS), were first reported by Nitsche & Paulus (2000), and later confirmed in motor cortex (Nitsche & Paulus, 2001) and somatosensory cortex as sensory evoked potentials (SEP) (Matsunaga, Nitsche, Tsuji, & Rothwell, 2004). Animal studies suggested that the elevated activity and excitability is not due to reverberating networks (Gartside, 1968a). Rather, changes in synaptic protein synthesis (Gartside, 1968b) point towards increased synaptic plasticity. In turn, blocking either brain-derived neurotrophic factor (BDNF; Fritsch et al., 2010), NMDA receptors (Nitsche et al., 2003), or calcium channels (Monte-Silva et al., 2013) leads to a reduction of the stimulation-induced increments of the field potential in mice, or MEP in humans. Recent evidence suggests that multiple forms of plasticity are in fact contributing to tDCS aftereffects. Monte-Silva et al. (2013) observed that fast facilitation, or early-LTP (e-LTP), was induced after a single tDCS session (13 min) and lasted for at least 2 hr after stimulation. In contrast, 26-min stimulation resulted in a reduced MEP amplitude. More interestingly, repetitive tDCS with 20-min pauses interspersed (13 − 20 − 13 min) resulted in late facilitation, or late-LTP (l-LTP). An elevated MEP was observed one day after the second stimulation, but not immediately after it. Functional LTP-like plastic changes of existing synapses were observed in DCS (Ranieri et al., 2012). Given the timescales of l-LTP, structural plasticity involving network remodeling also seems to play a role for the aftereffects. Structural changes at a slower timescale, however, can easily be underestimated because of difficulties measuring synapse turnover and changes in neuronal morphology in vivo. In summary, it is likely that both Hebbian and homeostatic as well as functional and structural forms of plasticity underlie tDCS aftereffects.
Quantitative models of network remodeling have previously been described in the literature. Butz, Steenbuck, & van Ooyen (2014) first introduced the term homeostatic structural plasticity with reference to previously published versions of the theory (Butz, van Ooyen, & Wörgötter, 2009; Butz & van Ooyen, 2013; van Ooyen, 2011), which was based on ample experimental evidence that structural plasticity (Trachtenberg et al., 2002; Oray, Majewska, & Sur, 2004; Holtmaat & Svoboda, 2009; Pfeiffer et al., 2018) as well as homeostatic regulation of activity (Turrigiano & Nelson, 2004; Lee et al., 2013; Keck et al., 2013) are constantly taking place in many brain areas. This homeostatic structural plasticity model was able to provide explanations for cortical reorganization after stroke (Butz et al., 2009) and lesion (Butz-Ostendorf & van Ooyen, 2017), and for the formation of certain global network features during development (Butz et al., 2014; Gallinaro & Rotter, 2018). In this model, changing the number of synaptic contacts between two neurons leads to an apparent facilitation or depression of this specific connection, and the model may therefore also account for some cases of functional plasticity. Based on these previous insights it seemed natural to explore the contribution of homeostatic structural plasticity to the long-lasting aftereffects of transcranial brain stimulation.
In the present work, we hypothesize that employing proper stimulation protocols and adequate current strengths, tDCS is potent enough to polarize single neurons in a network (Vöröslakos et al., 2018). Based on this assumption, we assess the effect of such membrane potential deflections on neuronal firing rates. In a neural network model with homeostatic structural plasticity, we then systematically explore the influence of various stimulation parameters known from tDCS practice, such as electrode size and montage, stimulation strength, and repetitive stimulation protocols. Our results suggest that tDCS can indeed induce substantial network remodeling and cell assembly formation, and focused strong and/or repetitive stimulation with well-chosen duty cycles can effectively boost the connectivity of the cell assemblies formed. The enhanced cell assembly might contribute to the empirical finding of profound plastic responses and enhanced therapeutic effects observed in current tDCS applications with a high-definition montage (Kuo et al., 2013) and repetitive stimulation (Monte-Silva et al., 2013). Our analysis also provides explanations for some of the negative results in tDCS practice.
Methods
Neuron Model
Model of Transcranial DC Stimulation
The electric field (EF) induced by tDCS can directly affect the membrane potential of neurons. Following Vöröslakos et al. (2018), we assumed that a strong enough EF will cause a small but significant membrane potential depolarization or hyperpolarization on some neurons in the network. The effective membrane potential deflection is determined by the orientation of the electric field vector relative to the somato-dendritic axis of the neuron (Wiethoff et al., 2014; Gluckman et al., 1996; Radman et al., 2009). When the electric field is properly aligned with the axis (apical dendrite closer to anode than soma), the somatic membrane potential is depolarized and the neuronal firing rate is increased. In contrast, if the electric field is perpendicular to the axis, it cannot influence the activity of this particular neuron. As a consequence, cells with extended and nonisotropic morphology, such as pyramidal neurons, should generally be more influenced by tDCS than the more compact inhibitory interneurons, which is also confirmed by Vöröslakos et al. (2018). Therefore, we assume only excitatory neurons to be sensitive to tDCS because of their spatial extent and nonisotropic morphology. We then asked whether such a polarization could also cause significant changes in the firing rate and, as a consequence, trigger structural plasticity and network remodeling. As our model neurons are actually point neurons with no spatial extent, we simply imposed an equivalent membrane potential bias ΔV on the soma of the neuron (Kayyali & Durand, 1991; Gluckman et al., 1996); see Figure 1A. This membrane potential bias also reflects the angle θ between the EF vector and somato-dendritic axis of the neuron with a factor ; see Figure 1B. The smallest magnitude of a membrane potential deflection reported in tDCS experiments to trigger physiological effects was in the range of 0.1 mV (Jackson et al., 2016; Vöröslakos et al., 2018).
Modeling the effect of tDCS on cortical networks. (A) It is assumed that transcranial stimulation leads to a weak polarization of the neuron’s membrane potential (left). For a point neuron, this is achieved by injecting a current of suitable strength into its soma (right). (B) Firing rate modulation with the angle θ for three different values of ΔVm (dotted lines on C). (C) Firing rate of a neuron, the ongoing activity of which is modulated by tDCS, for different values of θ and membrane polarization ΔVm. The contour lines correspond to 7 Hz, 8 Hz, and 9 Hz in white, orange, and maroon. (D) Electrode montages used in tDCS. (E) The region of interest subject to tDCS is modeled as a recurrent network of excitatory and inhibitory neurons. (F) Excitatory-to-excitatory synapses require the combination of a bouton (empty triangle) and a spine (red dot). The growth rate of both types of synaptic elements depends linearly on firing rate. (G) The network is grown from scratch before each tDCS stimulation experiment.
Modeling the effect of tDCS on cortical networks. (A) It is assumed that transcranial stimulation leads to a weak polarization of the neuron’s membrane potential (left). For a point neuron, this is achieved by injecting a current of suitable strength into its soma (right). (B) Firing rate modulation with the angle θ for three different values of ΔVm (dotted lines on C). (C) Firing rate of a neuron, the ongoing activity of which is modulated by tDCS, for different values of θ and membrane polarization ΔVm. The contour lines correspond to 7 Hz, 8 Hz, and 9 Hz in white, orange, and maroon. (D) Electrode montages used in tDCS. (E) The region of interest subject to tDCS is modeled as a recurrent network of excitatory and inhibitory neurons. (F) Excitatory-to-excitatory synapses require the combination of a bouton (empty triangle) and a spine (red dot). The growth rate of both types of synaptic elements depends linearly on firing rate. (G) The network is grown from scratch before each tDCS stimulation experiment.
Relative Strength of Background Activity and tDCS
The effect of tDCS on a neuron with ongoing activity was assessed with single neuron simulations. The background input impinging onto the neuron was approximated by a spike train with Poisson statistics and rate νext = 18.1 kHz, coupled to the neuron with synapses of strength Jext = 0.1 mV. Given the parameters of our neuron model, this ongoing background activity leads to a fluctuating subthreshold membrane potential with a mean value μ = νextτmJext = 18.1 mV (Brunel, 2000). Different values of membrane polarization caused by tDCS (from 0.1 mV to 1.2 mV) were considered in our study, as described above. The firing rate of each condition was estimated from simulations of 100-s duration.
Network Model
Although there are a variety of EF distributions induced by different tDCS montages, we assume simple uniform EF distributions in our model of the most affected area (Jackson et al., 2016). This most affected area is modeled as an inhibition-dominated recurrent network (Brunel, 2000), comprising 10,000 excitatory and 2,500 inhibitory neurons. All connections involving inhibitory neurons were taken to be static. Excitatory and inhibitory synapses had fixed synaptic weights of JE = 0.1 mV and JI = −0.8 mV, respectively. All these connections were randomly established, with 10% connection probability. In contrast, excitatory-to-excitatory (E-E) connections were subject to a growth rule called homeostatic structural plasticity (Gallinaro & Rotter, 2018; Butz & van Ooyen, 2013; Diaz-Pier, Naveau, Ostendorf, & Morrison, 2016). The network had initially no E-E connections whatsoever, and they were grown according to the specified rule during a growth period of 750 s for all simulations in the paper. Each neuron in the network received Poissonian external input at a rate of rext = 30 kHz. For the parameters chosen here, the network automatically entered an asynchronous-irregular state (Brunel, 2000). In all figures and simulations, transcranial DC stimulation was only applied after the end of the growth period. All network parameters are once more listed in Table 2.
Homeostatic Structural Plasticity
Protocols of Transcranial DC Stimulation
As suggested by current tDCS practice, many factors are essential to the outcome of a stimulation. For example, the traditional montage of two large sponge electrodes of size 5 cm × 7 cm induces a diffusive and weak EF. In contrast, high-definition montage using a small anodal electrode surrounded by several small cathodal electrodes induces a focal and relatively strong EF for the same stimulation current (Edwards et al., 2013). High-definition montage induces higher current densities, affects smaller populations, and possibly opposite field polarity at the edge of the cathodes. This method also exhibits better performance in tDCS practice, as compared with conventional montage (Kuo et al., 2013). To test these factors in our model we employed three different scenarios and systematically changed the size of the stimulated focus and the intensity of the stimulation in all of them. A summary of the parameters used in the different stimulation protocols described in this section can be found in Table 4.
Figure . | Protocol . | fG1 . | ΔV1[mV] . | fG2 . | ΔV2[mV] . | fG3 . | ΔV3[mV] . | Growth [s] . | Repetition . | tstim or t1[s] . | trelax or t2[s] . |
---|---|---|---|---|---|---|---|---|---|---|---|
2B | uni-group | 10% | 0.1 | 90% | 0 | − | − | 750 | no | 150 | 300 |
2D | uni-group | 10% | 0.1 | 90% | 0 | − | − | 750 | no | 150 | 300 |
3A | tri-group | 30% | 0.1 | 30% | − 0.1 | 40% | 0 | 750 | no | 150 | 300 |
3B | bi-group | 30% | 0.1 | 70% | − 0.1 | − | − | 750 | no | 150 | 300 |
4A | bi-group | 10%, 30%, 50%, 70% | −1.2, 1.2 1 | 1 − fG1 | − ΔV1 | − | − | 750 | no | 150 | 5,850 |
4B | uni-group | 10%, 30%, 50%, 70% | −1.2, 1.2 | 1 − fG1 | 0 | − | − | 750 | no | 150 | 5,850 |
4C | tri-group | 10%, 20%, 30%, 40% | −1.2, 1.2 | fG1 | − ΔV1 | 1 − fG1 − fG2 | 0 | 750 | no | 150 | 5,850 |
4I | all | 50% | −1.2, 1.2 | 50% | −1.2, 1.2 | − | − | 750 | no | 150 | 5,850 |
5D | repetitive | 10% | 0.1 | 90% | 0 | − | − | 750 | yes2 | multiple3 | |
5E | repetitive | 10% | multiple4 | 90% | 0 | − | − | 750 | 80 | 75 | 150 |
6-on-off | repetitive | 10% | 0.1 | 90% | 0 | − | − | 750 | 3 | 150 | 150 |
6-alternating1 | repetitive | 10% | ± 0.05 | 90% | 0 | − | − | 750 | 3 | 150 | 150 |
6-alternating2 | repetitive | 10% | ±0.1 | 90% | 0 | − | − | 750 | 3 | 150 | 150 |
Figure . | Protocol . | fG1 . | ΔV1[mV] . | fG2 . | ΔV2[mV] . | fG3 . | ΔV3[mV] . | Growth [s] . | Repetition . | tstim or t1[s] . | trelax or t2[s] . |
---|---|---|---|---|---|---|---|---|---|---|---|
2B | uni-group | 10% | 0.1 | 90% | 0 | − | − | 750 | no | 150 | 300 |
2D | uni-group | 10% | 0.1 | 90% | 0 | − | − | 750 | no | 150 | 300 |
3A | tri-group | 30% | 0.1 | 30% | − 0.1 | 40% | 0 | 750 | no | 150 | 300 |
3B | bi-group | 30% | 0.1 | 70% | − 0.1 | − | − | 750 | no | 150 | 300 |
4A | bi-group | 10%, 30%, 50%, 70% | −1.2, 1.2 1 | 1 − fG1 | − ΔV1 | − | − | 750 | no | 150 | 5,850 |
4B | uni-group | 10%, 30%, 50%, 70% | −1.2, 1.2 | 1 − fG1 | 0 | − | − | 750 | no | 150 | 5,850 |
4C | tri-group | 10%, 20%, 30%, 40% | −1.2, 1.2 | fG1 | − ΔV1 | 1 − fG1 − fG2 | 0 | 750 | no | 150 | 5,850 |
4I | all | 50% | −1.2, 1.2 | 50% | −1.2, 1.2 | − | − | 750 | no | 150 | 5,850 |
5D | repetitive | 10% | 0.1 | 90% | 0 | − | − | 750 | yes2 | multiple3 | |
5E | repetitive | 10% | multiple4 | 90% | 0 | − | − | 750 | 80 | 75 | 150 |
6-on-off | repetitive | 10% | 0.1 | 90% | 0 | − | − | 750 | 3 | 150 | 150 |
6-alternating1 | repetitive | 10% | ± 0.05 | 90% | 0 | − | − | 750 | 3 | 150 | 150 |
6-alternating2 | repetitive | 10% | ±0.1 | 90% | 0 | − | − | 750 | 3 | 150 | 150 |
The stimulation intensities are − 1.2, − 0.8, − 0.4, 0.4, 0.8, 1.2 mV
The repetition round (n) were matched with n×t1 = 6,000s
The combinations used are (75,75), (75,150), (150,75), (150,150), (150,300), (300,150) s
The stimulation intensities are 0.02, 0.04, 0.06, 0.08, 0.1, 0.2, 0.3, 0.4, 0.5 mV
Uni-group.
The first protocol we considered was a simplified scenario, in which only a subgroup of excitatory neurons in a large network was polarized by tDCS according to the above described protocol, while the remaining neurons were not affected and received only baseline external input. The focality of the stimulation is quantified by the percentage of excitatory neurons stimulated fG1. The more focused a stimulation is, the smaller is the subgroup of neurons affected by tDCS. The intensity of stimulation, on the other hand, is quantified by the amount of polarization. A stronger EF would lead to stronger membrane polarization of the soma of the model neurons. After a certain stimulation time tstim, tDCS is turned off and the network is allowed to relax for a period of trelax. Table 4 shows the values of fG1 and ΔV, as well as tstim and trelax, used for the different figures.
Bi-group.
Neurons in biological brains may not be uniformly polarized by stimulation. This is reflected in the bi-group scenario, in which a subgroup of neurons containing a fraction fG1 of all excitatory neurons is polarized by ΔV1 (similarly to the uni-group scenario), while the remaining excitatory neurons fG2 are stimulated with the same magnitude, but opposite polarity ΔV2. Similarly to the uni-group scenario, after a certain stimulation time tstim, tDCS is turned off and the network is simulated for a relaxation period trelax. The effect of stimulation on connectivity IG was calculated as described below.
Tri-group.
We designed yet another protocol, the tri-group scenario, to study the interaction of two actively stimulated subgroups with an unstimulated background. Two subgroups of excitatory neurons of the same size fG1 and fG2 are stimulated with the same magnitude, but different polarity ΔV1 and ΔV2. The remaining excitatory neurons in the network fG3 remain unstimulated. The resulting effect of stimulation on connectivity is measured as described below.
Repetitive patterns.
To examine the effects of repetitive on-off stimulation, a certain fraction fG1 of the excitatory neurons was stimulated in multiple cycles with the uni-group protocol. Each cycle corresponds to a stimulation period of length t1 followed by a pause of length t2. The number of cycles nc in each scenario was arranged to achieve a total DC stimulation time of nct1 = 6,000 s.
Repetitive alternating stimulation is similar to the repetitive on-off protocol based on the uni-group scenario. The difference is that, instead of pausing, neurons are stimulated with opposite polarity and same magnitude. In Table 4 we compiled a summary of all parameters for the stimulation protocols considered in our study.
Measurements and Calculations
Firing rate.
The firing rate of a neuron was calculated from its spike count, in a 5-s activity recording, unless stated otherwise. The mean firing rate of a population was taken to be the arithmetic mean of firing rates across neurons in the group.
Synaptic connectivity.
Let (Aij) be the n × n connectivity matrix of a network with n neurons. Its columns correspond to the axons, and its rows correspond to the dendrites of the neurons involved. The specific entry Aij of this matrix represents the total number of synapses from the presynaptic neuron j to the postsynaptic neuron i. The mean connectivity of this network is then given by , where t is the observing time point.
Time integral of the connectivity.
The parameter Ak is the amplitude of a component that decays with time constant τk. We then computed the total integral of the connectivity by integrating the sum of exponentials, amounting to (see Figure S2 in the Supporting Information). This way we can also account for connectivity transients that persist for a very long time, extrapolating beyond the duration of our simulations.
Results
Immediate Firing Rate Modulation by Transcranial DC Stimulation
We assume that the direct current applied to the brain during transcranial stimulation induces small deflections of the somatic membrane potential of neurons (Vöröslakos et al., 2018) and study the consequences of this deflection on neuronal firing rates. A polarization of the membrane ΔVm in the range between − 1.2 mV and 1.2 mV, which is not strong enough to elicit spikes in a neuron at rest, can nevertheless induce appreciable firing rate changes in a neuron with ongoing activity. Figures 1B and1C show how the firing rate of a model neuron driven by background input is modulated by both the strength of the depolarization and the angle θ between the electric field (EF) and the somato-dendritic axis. Even for a polarization as weak as ± 0.1 mV, which is about the weakest depolarization known to cause observable physiological effects in tDCS experiments (Vöröslakos et al., 2018), the firing rate change was found to be larger than ± 10% (Figure 1B, light gray curve). This suggests very clearly that tDCS can have an appreciable impact on neuronal activity, even if the stimulation intensity is apparently subthreshold. As neuronal spiking can affect synaptic connectivity via activity-dependent plasticity, this raises the question whether transcranial stimulation can trigger plastic effects as well. To find an answer to this question, we set up a plastic network representing the tissue most affected by tDCS (Figures 1E–G) and study the effect of stimulation.
Network Remodeling Triggered by Transcranial DC Stimulation
Different electrode montages are used in tDCS (Figure 1D), and they are thought to trigger different electric field distributions in the whole brain. We only modeled the most affected region stimulated by the peak current intensity. To explore the homeostatic response of the network, and the plastic processes associated with it, we first considered a simplified setting. In the uni-group scenario, only a subset of excitatory neurons in a larger network is stimulated (blue region in Figure 1E and Figure 2A). As shown in Figure 1F, tDCS disrupts the homeostatic equilibrium of the stimulated neurons by increasing their firing rate, initially leading to a deletion of synapses between stimulated neurons (see Methods for details of the structural plasticity model). When the stimulation has ceased, the firing rate of stimulated neurons drops because of a lack of recurrent input (Figure 2B), and the homeostatic process now triggers the formation of new synapses, predominantly among the stimulated neurons (Figure 2C). Figure 2F illustrates the process of cell assembly formation, similarly to what has been described previously by Gallinaro & Rotter (2018). Before and after the stimulation, assuming equilibrium in both cases, each neuron receives the same external input and fires at its target rate (here, 8 Hz). Thus, the total number of input synapses from excitatory neurons will not have changed through stimulation. What has changed, however, is the source of input synapses: Before stimulation, input comes from both groups of neurons—to be stimulated (blue) and background (empty)—without any bias. During stimulation, however, synapses are broken up, and when stimulation is turned off, the stimulated neurons have more free synaptic elements to offer. Background neurons, which are only indirectly affected by stimulation and deviate less from their target rate, can only offer a few synaptic elements to form new connections. Since the formation of new synapses is based on the availability of free elements, this leads to a bias for connections to be formed among stimulated neurons.
tDCS triggers the formation of cell assemblies. (A) A subgroup comprising 10% of all excitatory neurons in a larger network is stimulated by tDCS. (B) Average firing rate of directly stimulated (blue) and unstimulated (gray) excitatory neurons before, during, and after applying a depolarizing stimulus. (C) Average connectivity among stimulated neurons (blue), among unstimulated neurons (dark gray), and between neurons belonging to different groups (light gray) upon depolarizing stimulation. (D–E) Similar to (B–C), but for a hyperpolarizing stimulus. Shaded areas on (B–E) indicate the stimulation period. (F) Illustration explaining the process of structural plasticity that happened after a depolarizing tDCS. The stimulation triggers the removal of interpopulation synapses, and accelerates the growth of synapses among stimulated neurons, leading to the formation of cell assemblies.
tDCS triggers the formation of cell assemblies. (A) A subgroup comprising 10% of all excitatory neurons in a larger network is stimulated by tDCS. (B) Average firing rate of directly stimulated (blue) and unstimulated (gray) excitatory neurons before, during, and after applying a depolarizing stimulus. (C) Average connectivity among stimulated neurons (blue), among unstimulated neurons (dark gray), and between neurons belonging to different groups (light gray) upon depolarizing stimulation. (D–E) Similar to (B–C), but for a hyperpolarizing stimulus. Shaded areas on (B–E) indicate the stimulation period. (F) Illustration explaining the process of structural plasticity that happened after a depolarizing tDCS. The stimulation triggers the removal of interpopulation synapses, and accelerates the growth of synapses among stimulated neurons, leading to the formation of cell assemblies.
A similar process happens for hyperpolarizing DC (Figure 2D and2E). In this case, however, the connectivity among stimulated neurons increases during tDCS because of hyperpolarization and a resulting drop in firing rate. In summary, any perturbation to the equilibrium of the network firing rate dynamics, no matter whether it is depolarizing or hyperpolarizing, will trigger an increased synaptic turnover and network remodeling by deleting between-group synapses and forming new synapses within the stimulated group to form a cell assembly.
The Effect of Montage, Focality, and Intensity of Transcranial DC Stimulation
Stimulation is able to induce cell assembly formation in the uni-group scenario, as illustrated in Figure 2. However, neurons affected by tDCS might not be uniformly depolarized or hyperpolarized. Parameters like stimulation montage, focality, or intensity certainly influence the degree to which each neuron in the stimulated population is affected, and to what extent its membrane potential is depolarized or hyperpolarized. Therefore, we investigated two alternative stimulation scenarios that capture some of the complexities of neuron polarization in real tissue: tri-group stimulation (Figure 3A) and bi-group stimulation (Figure 3B), the details of which are described in the Methods section. Similarly to the simplest scenario illustrated in Figure 2, the stimulated neurons again form a cell assembly (Figures 3E and3F) also under more general conditions.
Interactions between subpopulations and cell assembly formation in more complex stimulation paradigms. (A) Tri-group scenario: 30% of all neurons in a network (G1) are depolarized by 0.1 mV, another 30% (G2) are hyperpolarized by − 0.1 mV, and the rest of 40% receives no stimulus. (B) Bi-group scenario: 30% (G1) are hyperpolarized by − 0.1 mV, and the remaining 70% (G2) are depolarized by 0.1 mV. (C,E) Group averages of firing rates in G1 (blue) and in G2 (yellow) before, during, and after stimulation. (D,F) Group averages of the connectivity within G1 (blue), within G2 (yellow), and between G1 and G2 (gray).
Interactions between subpopulations and cell assembly formation in more complex stimulation paradigms. (A) Tri-group scenario: 30% of all neurons in a network (G1) are depolarized by 0.1 mV, another 30% (G2) are hyperpolarized by − 0.1 mV, and the rest of 40% receives no stimulus. (B) Bi-group scenario: 30% (G1) are hyperpolarized by − 0.1 mV, and the remaining 70% (G2) are depolarized by 0.1 mV. (C,E) Group averages of firing rates in G1 (blue) and in G2 (yellow) before, during, and after stimulation. (D,F) Group averages of the connectivity within G1 (blue), within G2 (yellow), and between G1 and G2 (gray).
We performed a systematic study covering different degrees of stimulation focality and intensity and compared the effects in all three scenarios: the bi-group (Figure 4A), uni-group (Figure 4B), and tri-group (Figure 4C). Higher stimulus intensity is implemented as a stronger membrane polarization, which results from a higher tDCS current density. Focality, quantified as the percentage of neurons in the network affected by membrane polarization, describes how focused stimulation is. More focused stimulation should have a polarizing effect on a smaller percentage of neurons. In each scenario, the connectivity in a newly formed cell assembly increases with absolute stimulation intensity and decreases with the size of the stimulated population (Figures 4D–F). We conclude that strong and focused stimulation (like high-definition stimulation) leads to stronger effects on the connectivity of the cell assembly. We further compared the effects of bi-group, uni-group, and tri-group scenarios and found that the montage can greatly influence the outcome. When the polarization is very strong (above 0.8 mV) and focused, the effect is much stronger in the uni-group scenario as compared with the bi-group (Figure 4G) and tri-group (Figure 4H) scenario. But if the stimulus is weak, its effect in the bi-group scenario is larger than in the uni-group scenario. Therefore, using opposite polarities for stimulation could slightly boost cell assembly formation, provided the stimulus is weak. However, for strong and/or focused stimulation, uni-group stimulation leads to more pronounced cell assemblies.
Comparison of tDCS effects with different electrode montage, as well as stimulus focality and intensity. (A) Bi-group stimulation scenario. (B) Uni-group stimulation scenario. (C) Tri-group stimulation scenario. (D–F) Integrated G1 cell assembly connectivity () at different focality and intensity levels for scenarios (A–C). (G,H) Difference between D and E, as well as F and E, respectively. (I) Integrated G1 cell assembly connectivity integrals () for different stimulation intensity levels for a specialized bi-group scenario, where G1 and G2 comprise half of the excitatory population, respectively. The white squares correspond to situations where the difference between stimulation intensities of both groups amounts to 0.8 mV.
Comparison of tDCS effects with different electrode montage, as well as stimulus focality and intensity. (A) Bi-group stimulation scenario. (B) Uni-group stimulation scenario. (C) Tri-group stimulation scenario. (D–F) Integrated G1 cell assembly connectivity () at different focality and intensity levels for scenarios (A–C). (G,H) Difference between D and E, as well as F and E, respectively. (I) Integrated G1 cell assembly connectivity integrals () for different stimulation intensity levels for a specialized bi-group scenario, where G1 and G2 comprise half of the excitatory population, respectively. The white squares correspond to situations where the difference between stimulation intensities of both groups amounts to 0.8 mV.
The application of hyperpolarizing DC to all neurons in the background population can either amplify or attenuate the effect of the actual depolarizing stimulus. Two aspects might contribute to this phenomenon. Stimulating the background with reversed polarity increases the discrepancy of the stimulated group compared with the background (from ΔVm to 2ΔVm), but it may reduce the firing because of an activation of inhibitory neurons in the network. To disentangle the situation, we fixed the sizes of both the stimulated and the unstimulated group at 50% and then systematically changed the stimulus strength for both G1 and G2 in the range between − 1.2 mV and 1.2 mV. The effect on G1 connectivity for different polarizations of G1 and G2 is displayed in Figure 4I. The values along the diagonal are very small, as there is no cell assembly formation when both groups experience the same stimulation. When the difference in stimulation of the two populations is large irrespective of its sign, the impact on G1 connectivity is also large (upper left and bottom right corners). We then checked whether the relative difference between the polarization of G1 and G2 is sufficient to predict the stimulation outcome. The white squares in Figure 4I indicate simulations in which the difference between G1 and G2 polarization is the same (0.8 mV), but the actual connectivities for individual groups are different. The strongest effect was achieved when the polarization of one of the two groups is 0 mV, which corresponds to the uni-group scenario. This supports the idea that network effects might influence the interaction between two groups, and that uni-group stimulation can achieve better outcomes than alterantive scenarios, provided stimulation is very strong.
The Effect of Repetitive Transcranial DC Stimulation
Repetitive stimulation was simulated in our model by repeating stimulation of duration t1 in the uni-group scenario (Figure 2) multiple times, with a pause of duration t2 between successive stimulation periods (Figures 5A and5B). The connectivity of the stimulated subpopulation generally increased upon repetition (Figure 5C). Figure 5D summarizes the outcome of different combinations of t1 and t2. Compared with long uninterrupted DC stimulus (single stimulation cycle with t1 = 6,000 s), repetitive stimulation (total stimulation time of 6,000 s distributed over multiple cycles of shorter duration t1) led to higher final connectivity. We found that repetitive stimulation generally potentiated the effect of tDCS on cell assembly connectivity. Figure 5E demonstrates that after multiple repetitions, the connectivity appears to saturate at a level that essentially depends on the imposed polarization. As a consequence, a single stimulation with weak intensity for a very long time does not necessarily lead to high connectivity, while repetitive stimulation at high intensity may lead to (much) higher connectivity (see Figure S1 in the Supporting Information). In our model we also tried very strong stimulation, repeated for several rounds. This led to a very high assembly connectivity and eventually also to a very high firing rate of the excitatory population. High firing rates, in turn, induced a strong homeostatic response of the network and fast deletion of synapses, putting the network in an unfavorable and somewhat pathological state (data not shown).
Repetitive stimulation boosts network remodeling. (A) A subnetwork of excitatory neurons (10%) is stimulated with a train of DC stimuli. Stimulation time is t1, followed by a pause of duration t2. (B, C) Average firing rate and connectivity during a train of stimuli. (D) For the same total stimulation time (6,000 s), the boosting depends on the exact repetition protocol. (E) The peak connectivity reached depends on the stimulation intensity; an asymmetric repetitive protocol (t1 = 75 s, t2 = 150 s) was used for all simulations here.
Repetitive stimulation boosts network remodeling. (A) A subnetwork of excitatory neurons (10%) is stimulated with a train of DC stimuli. Stimulation time is t1, followed by a pause of duration t2. (B, C) Average firing rate and connectivity during a train of stimuli. (D) For the same total stimulation time (6,000 s), the boosting depends on the exact repetition protocol. (E) The peak connectivity reached depends on the stimulation intensity; an asymmetric repetitive protocol (t1 = 75 s, t2 = 150 s) was used for all simulations here.
Repetitive stimulation can also be performed in cycles of alternating polarities, instead of a simple on-off protocol. Figure 6 shows the connectivity changes for two stimulation patterns: on-off, in which periods of depolarizing stimulation are followed by periods of no stimulation, and alternating, in which periods of depolarization are followed by periods of hyperpolarization. Simply substituting the off period by stimulation with different polarity seems to boost cell assembly connectivity (compare light green and dark brown traces in Figure 6B). However, if the alternating pattern has the same overall amplitude as the on-off stimulation (compare light brown and green traces in Figure 6B), the effect on cell assembly connectivity is the same as on the on-off pattern. Figure 6C depicts the final connectivity after three repetitions in 30 independent trials (mean and standard deviation are indicated in the inset).
Comparison of three different scenarios for repetitive DC stimulation. (A) 10% of excitatory neurons were stimulated, using the same temporal protocol (t1 = 150 s, t2 = 150 s) in each case, but different amplitudes and polarities were employed, as indicated by the three different curves. (B) Evolution of average connectivity for the different stimulation scenarios; color code matches the stimulus curves in panel A. (C) Histograms of the connectivity reached after three cycles in the different scenarios extracted from 30 independent depolarizing simulations; mean values and standard deviations are shown in the inset.
Comparison of three different scenarios for repetitive DC stimulation. (A) 10% of excitatory neurons were stimulated, using the same temporal protocol (t1 = 150 s, t2 = 150 s) in each case, but different amplitudes and polarities were employed, as indicated by the three different curves. (B) Evolution of average connectivity for the different stimulation scenarios; color code matches the stimulus curves in panel A. (C) Histograms of the connectivity reached after three cycles in the different scenarios extracted from 30 independent depolarizing simulations; mean values and standard deviations are shown in the inset.
Discussion
We explored the plastic changes in network structure that can be induced by transcranial direct current stimulation (tDCS), exploiting the homeostatic response of synaptic growth and decay. We demonstrated that weak subthreshold DC stimulation induces changes of neuronal firing rates and, thus, triggers network remodeling and cell assembly formation. Depolarized neurons first reduce the number of excitatory input synapses during stimulation, but then create new excitatory synapses predominantly with other stimulated neurons after stimulation is off. Interestingly, hyperpolarization also causes new synapses being formed preferentially among stimulated neurons. Stimulation triggers a profound and sustainable reorganization of network connectivity and leads to the formation of cell assemblies. With the help of our model, we explored different parameters of tDCS stimulation and found that strong and focused stimulation generally enhances the newly formed cell assemblies. We also observed that repetitive stimulation with well-chosen duty cycles can boost the induction of structural changes, and repetitive stimulation with alternating polarization may induce even higher connectivity changes.
We used network connectivity as a direct readout of stimulation effects, which is possible in model simulations, but cannot easily be done in experiments. However, the factors that we found to amplify the overall impact of stimulation are not unheard of in tDCS practice. Strong and focused stimulation, for example, which results from a high-definition electrode montage, does indeed lead to a stronger readout (MEP) and potentiates the therapeutic effects as compared with a conventional montage (Kuo et al., 2013). While applying the same total current, a high-definition montage induces stronger electric fields in smaller brain volumes (Edwards et al., 2013). Moreover, a high-definition montage narrows down the most affected brain region. We also found in our model that both factors indeed contribute to the induction of higher connectivity. Moreover, repetitive stimulation can boost connectivity, provided the duty cycles are chosen right. In fact, it has been demonstrated in experiments (Monte-Silva et al., 2013) that two 13-min stimulations interrupted by a 20-min pause yields stronger MEP aftereffects than a single, uninterrupted 26-min stimulation, while a repetition with a 24-h pause in between could not accumulate the aftereffects at all. In our model, we likewise found that multiple stimulation episodes with properly chosen pauses can achieve better effects than a single, uninterrupted stimulation.
Other computational approaches have been employed previously to analyze the neuron-scale mechanisms underlying tDCS or DCS. Most notably, Bikson et al. (2006) has explored several aspects of this: extracellular potassium concentration, polarization of the axonal terminal, action potential timing, and inhibitory neurons. Joucla & Yvert (2009) provided an estimate of membrane potential changes for large axons exposed to an electric field, and Aspart, Ladenbauer, & Obermayer (2016) conceived the influence of the electric field on neuronal dendrites as external input to the soma. Another computational approach based on modern neural imaging methods has shed light on the question of how strong the stimulation effects actually are. Spherical head models were first used to estimate the 3-D current flow for any given electrode montage (Miranda, Lomarev, & Hallett, 2006). Later, fMRI-based modeling was employed to devise individualized treatment of stroke or depressive patients (Datta et al., 2009; Ho et al., 2014; Huang et al., 2017). Our present work adopted insight and parameters from both approaches. In addition, we developed a new and original computational model to explore the impact of structural plasticity at the level of networks. This provides a bridge between the level of single neurons and the level of large-scale networks. Although our model contributes new explanations for some core observations in tDCS practice, there are still important issues left that cannot be appropriately addressed with our highly simplified model lacking relevant features of brain geometry. Also, the exact rules of growth and the timescales involved in homeostatic structural plasticity remain to be elucidated in experiments. To treat the influence of tDCS on network dynamics and structural plasticity of multiple brain regions would require a “network of networks” approach, which is, however, beyond the scope of our current study.
What are the actual effects of tDCS on network activity and function? Although robust and sustained effects of tDCS using relatively weak stimulation currents (1–2 mA) have been demonstrated (Nitsche et al., 2009; Nitsche & Paulus, 2000), Horvath et al. (2015) pointed to the difficulty reproducing positive results. Recently, Vöröslakos et al. (2018) have shown that the amount of membrane polarization due to tDCS depends on the strength of the applied current, and that there should be indeed no effect expected for very low intensities. Our simulation results suggest, however, that repetition could boost the impact on connectivity. The peak connectivity reached after sufficiently many repetitions, however, depends on stimulus intensity. Very weak stimulation cannot achieve high connectivity changes, even if repeated ad infinitum. Strong stimulation within a safe range could achieve higher connectivity, but too strong stimulation may lead to unfavorable network dynamics. Our model predicts very clearly that the accumulated effect achieved by stimulation depends not only on the exact repetition pattern, but also on stimulation intensity. On the other hand, a quantitative assessment of the aftereffects is difficult. In our work, the effect of tDCS on the network is quantified by measuring anatomical connectivity among stimulated and nonstimulated neurons. Such measurement is currently not possible in experiments, neither in vivo nor in vitro. Transcranial stimulation perturbs neuronal firing rates transiently and leads to the formation of cell assemblies, which persist after tDCS has been switched off and neuronal activity is back to baseline. Therefore, considering the homeostatic nature of structural plasticity, it is actually impossible to measure the effect of tDCS using simple neuronal activity measures. The question is, what are the effects of altered connectivity on the activity and the function of neuronal networks, and how can these effects be measured? This is a very interesting question, and the answer is complicated. Even if newly formed cell assemblies do not affect spontaneous activity as the firing rate of the neurons may be homeostatically regulated, they might still influence the evoked responses of neurons. Interestingly, Horvath et al. (2015) reviewed many tDCS studies and found that stimulation has a reliable effect only on the MEP amplitude, out of many potential biomarkers that were tested. The debate about the effects of tDCS on network function should, therefore, include the measures to quantify the outcome of a stimulation.
Another important issue raised by our work is that the total effect of stimulation might be too weak for detection. The connectivity changes triggered by a single cycle of polarization at ΔVm = 0.1mV can only be detected if the full connectome is available for quantification. While possible in simulations, such a scenario is unrealistic in an experimental setting. Our simulation results suggest, however, that the outcome should increase upon repetitive stimulation and, therefore, possibly becomes easier to measure. The measurement time window of tDCS effects adds another puzzle to this question. The connectivity of the stimulated plastic network undergoes constant changes. During and after stimulation, for instance, total connectivity decreases and increases fast, constituting the homeostatic response. In contrast, the newly formed cell assembly persists for much longer periods and decays only with a slower time constant. It is not yet clear, however, which parameters influence this time constant, and it might be that different current intensity and electrode size have an impact on it. In fact, Jamil et al. (2017) recently observed in experiments that the current intensity might interact with the duration of stimulation needed for the homeostatic reversal of plasticity. If the exact stimulation protocol indeed influences the timescale of the aftereffect, naively comparing tDCS effects under different stimulation conditions “before” and “after” does not provide sufficient information regarding its outcome. In view of this, using a measure that takes the dynamics of the changes triggered by stimulation into account, such as the IG measure introduced in this work, could quantify the effects of stimulation much more reliably.
In general, one needs to interpret the results and predictions of our work on network remodeling induced by tDCS with due caution. Our current work, however, could be a first step toward the goal of understanding and optimizing tDCS performance. More experiments addressing the impact of tDCS in human and in animal brains are definitely needed, and the results of our simulation study might indicate some new directions.
SUPPORTING INFORMATION
Supporting information for this article is available at http://doi.org/10.1162/netn_a_00097.
AUTHOR CONTRIBUTIONS
Han Lu: Conceptualization; Data curation; Formal analysis; Investigation; Methodology; Project administration; Resources; Visualization; Writing - Original Draft. Jlia V. Gallinaro: Conceptualization; Formal analysis; Investigation; Methodology; Project administration; Resources; Software; Writing - Review & Editing. Stefan Rotter: Conceptualization; Funding acquisition; Investigation; Methodology; Project administration; Resources; Supervision; Validation; Writing - Review & Editing.
FUNDING INFORMATION
Han Lu, Universitätsklinikum Freiburg. Han Lu, Neurex. Júlia V. Gallinaro, Erasmus Mundus/ EuroSPIN. Stefan Rotter, Carl Zeiss Stiftung. Stefan Rotter, Deutsche Forschungsgemeinschaft (http://dx.doi.org/10.13039/501100001659), Award ID: INST 39/963-1 FUGG. Stefan Rotter, Deutsche Forschungsgemeinschaft, Award ID: EXC 1086. Stefan Rotter, University of Freiburg, Award ID: open access publication fund.
ACKNOWLEDGMENTS
The authors thank Claus Normann, Lukas Frase, Andre Russowsky Brunoni, Sandra Diaz-Pier, and Benjamin Merkt for useful discussions. We also thank Uwe Grauer from the Bernstein Center Freiburg as well as Bernd Wiebelt and Michael Janczyk from the Freiburg University Computing Center for their assistance with HPC applications.
TECHNICAL TERMS
- Somato-dendritic axis:
Imaginary axis linking the soma and the distral dendrites of a pyramidal neuron.
- PET:
Positron-emission tomography, based on a radioactive tracer that is introduced into the body before the induction of metabolic processes.
- MEP/SEP:
The motor evoked potential or sensory evoked potential is a signal recorded from muscles or from the surface of the brain, respectively, as a result of a stimulus.
- BDNF:
Growth factor that acts on certain neurons. It supports survival of neurons, growth, and differentiation of new neurons and synapses.
- Cell assembly:
A group of neurons with enhanced mutual synaptic connectivity; arises after repetitive activation.
- High-definition montage:
Instead of using large sponge electrodes, high-definition montage employs smaller gel-based electrodes to allow for more focal stimulation.
- Point neuron:
A simplified neuron model that only represents the somatic membrane potential, disregarding the properties of spatially extended dendrites and axons.
- Poisson process:
Simple statistical model to account for a barrage of spikes, commonly used to describe the synaptic inputs driving a neuron.
- Set-point hypothesis:
In the model of homeostatic structural plasticity, neural firing rate is actively maintained at a set-point by adjusting synaptic input and output.
- Spherical head model:
Standard three-shell head model, which omits individual features, and includes the brain, intermediate skull, and outer-layer scalp.
REFERENCES
Author notes
Competing Interests: The authors have declared that no competing interests exist.
Handling Editor: Sandra Diaz