- 1Research Academy of Medicine Combining Sports, Ningbo No. 2 Hospital, Ningbo, China
- 2Faculty of Sports Science, Ningbo University, Ningbo, China
- 3Key Laboratory of Exercise and Health Sciences of Ministry of Education, Shanghai University of Sport, Shanghai, China
- 4Department of Neurology, University Medical Center Göttingen, Georg-August University, Göttingen, Germany
- 5Department of Rehabilitation Medicine, Ningbo No. 2 Hospital, Ningbo, China
Background: Multifocal transcranial direct current stimulation (m-tDCS) may modulate distributed motor networks in a polarity-dependent and task-state-dependent manner to support performance near exhaustion. To this end, this study aims to test whether m-tDCS targeting lower-limb-specific cortical areas could optimize late-stage performance and phase-specific muscle coordination during cycling.
Methods: Two independent trials were conducted: (i) a tolerability assessment (Trial 1) and (ii) a randomized, double-blind, sham-controlled parallel study (Trial 2). In Trial 1, participants completed the tolerability test and recorded pain and side effects during a 21-min stimulation period. In Trial 2, healthy adults completed an incremental cycling test; late-stage performance was operationalized a priori as the 85%–100% peak power output (PPO) phase, during which the time-to-exhaustion (TTE), work (W), mean power (P), revolutions per minute (RPM), heart rate (HR), blood lactate level (ΔL), ratings of perceived exertion (RPE), and EMG-derived muscle contribution ratio (MCR) and knee co-activation index (CAI) were analyzed across propulsion and pull.
Results: (1) m-tDCS was well-tolerated; pain ratings declined progressively across the stimulation, with typical transient sensations. (2) At the 85%–100% PPO phase, m-tDCS increased W and RPM relative to sham without altering the mean power or ΔL; HR decreased after m-tDCS, and RPE rose only after sham. (3) At the coordination level, m-tDCS preserved quadriceps MCR during propulsion and reduced antagonistic activation during pull to prevent the CAI increase observed in the sham.
Conclusion: m-tDCS did not augment peak mechanical output but preserved late-stage endurance via phase-specific coordination tuning, which is consistent with improved neural efficiency near exhaustion. These findings refine mechanistic interpretations of the effects of tDCS on endurance and support m-tDCS as a safe, coordination-centric adjunct for high-intensity cycling.
1 Introduction
Endurance performance in cycling is influenced by the interplay between central drive, peripheral fatigue, and neuromuscular coordination, especially during the high-intensity terminal phases of exercise. A growing body of research indicates that non-invasive brain stimulation (NIBS) can exert promising ergogenic effects on motor capacity across diverse populations; however, these effects are modality- and protocol-dependent, with different modalities of NIBS producing distinct physiological responses.
Transcranial direct current stimulation (tDCS), one of the most widely used and safe NIBS techniques in sports science and rehabilitation, modulates cortical excitability and network dynamics in a polarity-specific and state-dependent manner. Its ergogenic effects are confirmed to be gated by the ongoing “task state,” such that the same current can yield different behavioral consequences at rest compared with during an engaged motor task (Nitsche and Paulus, 2000; Stagg and Nitsche, 2011; Silvanto et al., 2008; Polanía et al., 2018). In contrast, transcranial magnetic stimulation (TMS), another NIBS technique, delivers brief, focal magnetic pulses that directly evoke action potentials; it is also commonly used to probe corticospinal function and central fatigue rather than as an on-task ergogenic intervention (Hallett, 2007; Lefaucheur et al., 2020). Thus, TMS should be considered a tool to provide complementary mechanistic benchmarks, whereas tDCS may offer a practical means to optimize motor control during fatiguing activity (Nitsche and Paulus, 2000; Nitsche and Paulus, 2001; Stagg and Nitsche, 2011; Polanía et al., 2018).
However, evidence supporting the beneficial effects of tDCS on exercise performance remains mixed: some studies report improvements in endurance performance and perceived exertion (e.g., longer time to exhaustion or reduced perceived effort) (Vitor-Costa et al., 2015; Okano et al., 2015; Angius et al., 2016), whereas meta-analyses indicate generally small, protocol-sensitive effects with substantial heterogeneity across tasks (Machado et al., 2019; Holgado et al., 2019; Kaushalya et al., 2022; Holgado et al., 2024; Flor et al., 2025). One approach to reduce variability and enhance functional specificity is to use multifocal electrode montages that target distributed motor networks rather than a single scalp site; both modeling work and empirical work suggest that multi-electrode configurations can improve focality and effective intensity at the target while engaging interconnected cortical regions (Dmochowski et al., 2011; Huang and Parra, 2019).
Cycling performance requires phase-specific coordination of muscle synergies, with task-relevant contributions during the propulsion phase and potential antagonistic activation during the pull phase. Electromyographic studies have demonstrated robust muscle synergy structures and rate-dependent coordination (Neptune et al., 1997; Hug and Dorel, 2009; Hug et al., 2011). From this perspective, indices such as the muscle contribution ratio (MCR) and co-activation index (CAI) provide concise and interpretable summaries of how neuromuscular effort is distributed across different phases, capturing both task-relevant activation and the inhibition of non-task-related co-activation.
Moreover, growing evidence suggests that prefrontal cortical mechanisms play a significant role in regulating effort and terminating tasks during endurance exercise (Qi et al., 2024). This supports the notion that stimulation strategies aimed at optimizing network efficiency may influence tolerance at high intensities (Robertson and Marino, 2016; De Wachter et al., 2021). Our previous research (Figure 1) has also found a significant correlation between the activation of the prefrontal and motor cortices and lower-limb performance during incremental loading to exhaustion (Liang, 2024). Building on these findings, multifocal tDCS (m-tDCS) targeting these key regions could potentially produce network-level modulation effects.
Therefore, this study tested whether m-tDCS targeting lower-limb-specific cortical areas could optimize late-stage performance and phase-specific muscle coordination during cycling via focusing on terminal-phase behavioral outcomes and on MCR and CAI across propulsion and pull phases as mechanistic endpoints. Two independent trials were conducted. Trial 1 tested the tolerability of the m-tDCS protocol, whereas Trial 2 evaluated its effects on lower-limb cycling endurance motor performance during an incremental loading cycling exercise using multimodal performance measures such as electromyography, physiological responses, and subjective effort ratings.
2 Methods
2.1 Participants
Sixty healthy young adults who met the eligibility criteria were recruited to participate in Trial 1 and Trial 2. The eligibility criteria were established based on Liang (2024).
Trial 1. A total of 15 healthy young adults voluntarily participated in the tolerance test for tailored multifocal stimulation (Table 1); to avoid participants guessing the experimental purpose and optimize blind effects, participants in Trial 1 were not allowed to take part in Trial 2.
Trial 2. A total of 45 healthy young adults with the same characteristics participated in the tDCS and endurance exercise tests. Participants were randomly allocated to two groups to either receive one session of m-tDCS or one session of sham stimulation. Before trial testing, all participants had to pass a medical screening, were informed of the procedures, and provided informed written consent. The study was approved by the Ethics Committee of the Shanghai University of Sport under the registration number 102772023RT031, and its protocol adhered to the 1964 Declaration of Helsinki and its revisions.
2.2 Experiment protocol
Two trials were used to evaluate the m-tDCS effects on endurance capacity. Trial 1 was designed to verify the physical tolerance on tailored tDCS, whereas Trial 2, a randomized, sham-controlled, and double-blinded study, was designed to investigate the effects of m-tDCS on lower-limb cycling endurance motor performance.
During Trial 1, participants received 21 min of multifocal tDCS while sitting; the physical tolerance and side effects at different time points throughout the stimulation were monitored concurrently using the pain and side effect scales. During Trial 2, participants were asked to visit the laboratory on two consecutive days. The first visit was to investigate the baseline for maximum endurance capacity on the cycling ergometer. After 24 h, the participants visited the laboratory for the second visit. On the second visit, subjects first received either a 21-min m-tDCS or sham stimulation and then completed the incremental loading cycling immediately after stimulation, during which the physiological performance, psychological performance, and dominant muscle activity at termination phases were also measured to evaluate cycling endurance performance. To achieve the best motor performance, participants were verbally encouraged throughout the test.
2.3 Multifocal transcranial direct current stimulation
The Neuroelectrics Starstim stimulator (Neuroelectrics, Barcelona, Spain) was used to deliver m-tDCS and sham stimulation through a set of seven circular Ag/AgCl electrodes with a surface area of 1.13 cm2. Based upon our previous findings (Liang et al., 2024), seven electrodes were placed in the C1, C2, FC1, FC2, AF7, P9, and P07 positions according to the 10/10 international EEG system, during which the participants received an anodal current over bilateral motor cortices (C1 and C2) and bilateral superior cortices (FC1 and FC2) and a cathodal current over the left lateral occipital cortex (P07 and P9) and left frontal pole (AF7) (Figure 2). The current intensity at each electrode was set to 1.053 mA from C1, 0.77 mA from C2, 1.221 mA from FC1, 0.954 mA from FC2, −1.016 mA from AF7, −1.430 mA from P9, and −1.552 mA from P07; the corresponding current densities were 0.93, 0.68, 1.08, 0.84, 0.90, 1.27, and 1.37 mA·cm−2, respectively. To mimic cutaneous sensations while maintaining blinding, the stimulator in the sham condition delivered a 30-s ramp-up immediately followed by a 30-s ramp-down to 0 mA, after which the current was set to zero for the remainder of the session. During stimulation, both the investigator and the participants were blinded to the type of stimulation the participants received.
2.4 Simulation of electric field characteristics for the stimulation protocol
SimNIBS 4.0 with the algorithm and finite element model was used to calculate the electric field of m-tDCS protocol based on Ruffini et al. (2014). MRI data, the Colin 27, were segmented into the scalp, skull, cerebrospinal fluid, gray matter, and white matter. Conductivities of the anatomical tissues were set to σ-skin = 0.330 S/m; σ-bone = 0.008 S/m; σ-cerebrospinal fluid = 1.790 S/m; σ-gray matter = 0.400 S/m; and σ-white matter = 0.150 S/m. Electrodes were set up as cylinders with a radius of 10 mm and a thickness of 3 mm. Every tissue was regarded as isotropic and homogeneous. In order to best fit the spatial topography of the targeted brain regions, the target cortical solution was represented by the error relative to no intervention, a metric serving as a measure of quality of fit. The current intensity at each electrode was set up following the above m-tDCS protocol (Figure 3).
2.5 Tolerability and safety of the stimulation protocol
The pain perceived and side effects during stimulation were monitored to evaluate the tolerability and safety of m-tDCS. The participants were asked to report the subjective pain level and any of the adverse effects under the electrodes at the stimulation onset (30-s up); at 30 s, 1 min, 5 min, 10 min, 15 min, and 20 min; and at the end of stimulation (30-s down), respectively. The pain intensity was rated on an 11-point numerical rating scale (NRS) consisting of integers ranging from 0 to 10 (0: no pain; 3–7: moderate pain; and 7–10: worst pain possible). The tingling side effect under the electrodes was monitored, including headache, itching, burning sensation under electrodes, sleepiness, trouble concentrating, and acute mood change.
2.6 Incremental loading exercise protocol
The incremental loading exercise was adopted as the motor task to determine the dominant lower-limb maximum capacity of each participant. Participants had to cycle for 3 min without any load on a cycle ergometer (282E, Monark, Sweden), and subsequently completed the incremental cycling. During cycling, the incremental loading was increased by 20 W per minute for female participants and by 25 W per minute for male participants until voluntary exhaustion. During cycling, the criterion to stop the exercise was in keeping with Liang (2024). The cycling time and the peak power output (PPO) were recorded simultaneously. Based on a previous study, we divided the whole incremental loading exercise into four phases: 0%–50% PPO, 50%–75% PPO, 75%–85% PPO, and 85%–100% PPO (O’Bryan et al., 2014). The performance at 85%–100% PPO, corresponding to the maximum motor capacity, was utilized to evaluate effects of m-tDCS on lower-limb motor capacity via assessing time-to-exhaustion (TTE), work (W), average power (P), and revolutions per minute (RPM) throughout this phase.
2.7 Evaluation of physiological and psychological performances
The heart rate (HR), blood lactate level (ΔL), and rating of perceived effort (RPE) were monitored to assess physiological and psychological performances at the termination phase. HR was measured using a wireless heart-rate belt (Polar H10, Polar Electro Oy, Finland), with the sensor worn at the sternal position and the elastic band adjusted to ensure the constant positioning of the receiver during cycling. Furthermore, the blood lactate level (ΔL) was measured using the portable Lactate Scout 4 Solo analyzer (EKF Diagnostics, Germany) before and after the incremental loading exercise intervention. During testing, a second blood sample was taken at the right earlobe to avoid unnecessary distractions (e.g., sweat).
To assess the psychological performance related to subjective feelings of fatigue and loading (Borg, 1982), the 20-scale Borg rating of perceived effort was implemented to quantify exercise fatigue at each incremental loading; the collection time was consistent with the participants’ HR at that time, which was the last 10 s of each loading.
2.8 Evaluation of lower-limb muscle activity
Surface electromyography (sEMG) of the dominant lower limb as the key parameter was used to evaluate the effect of m-tDCS on lower-limb motor capacity, data for which were collected from eight active lower-limb muscles during cycling using an Ultium EMG system (Noraxon, United States) with a sampling frequency of 2,000 Hz. The monitored muscles included the rectus femoris (RF), vastus lateralis (VL), vastus medialis (VM), tibialis anterior (TA), long head of the biceps femoris (BF), semitendinosus (ST), medial gastrocnemius (GM), and lateral gastrocnemius (GL). For each loading phase of cycling, sEMG data were recorded starting from the 10th second of the loading phase and continued for a duration of 20 s, and the cycling period was divided into two phases: the propulsion and pull phases based on the recorded ankle joint motion. The sEMG of the dominant lower-limb muscle was analyzed for root mean square (RMS) amplitude, which was used to calculate the muscle contribution ratio (MCR) and co-activation index (CAI) of the knee joint at 85%–100% PPO, to further qualify the m-tDCS effects on lower-limb motor capacity.
2.9 Data processing
2.9.1 Peak power output processing
The calculation of PPO adopted the following formulas that have been previously published in Jeukendrup et al. (1996):
Here, Wout is the last loading value in watts (w) and t is the cycling lasting time during the last loading in seconds (s).
2.9.2 Electromyography data processing
The sEMG data were processed using the Visual 3D™ software (C-Motion, United States). The specific processing steps for the sEMG signals are as follows: DC offset removal, band-pass filtering (10 Hz–400 Hz, second-order Butterworth filter), full-wave rectification, and moving average filtering (window width: 250 ms) to obtain the linear envelope of the EMG signal with a cutoff frequency of 6 Hz. RMS values for the propulsion and pull phases were calculated based on the motion trajectory of the ankle joint along the Z-axis, averaged over 10 consecutive cycling cycles. The RMS values were then normalized using the maximum value from the maximal voluntary contraction (MVC). Finally, the data from the 10 cycling cycles were averaged to provide the final RMS value for the statistical analysis of the muscle activity level, the muscle contribution ratio, and CAI of the knee joint. The formulas for MCR and CAI of the knee joint are as follows:
The MCR was calculated to assess the relative contribution of each of the eight monitored muscles to the overall muscle activity during cycling. The MCR was defined as
where RMSi refers to the RMS value of the EMG signal for the ith muscle, representing the level of muscle activation during cycling.
The denominator represents the sum of the RMS values of all eight monitored muscles, which gives the total muscle activity across all muscles.
The CAI was calculated to assess the level of co-activation between the agonist and antagonist muscles during knee flexion and extension. The CAI was defined as
where “RMS of Agonist” refers to the RMS value of the EMG signal of the agonist muscle (e.g., quadriceps for knee extension) and “RMS of Antagonist” refers to the RMS value of the EMG signal of the antagonist muscle (e.g., hamstrings for knee extension).
The CAI ranges from 0 to 1, where a value close to 0 indicates minimal co-activation of the antagonist muscle, suggesting that the agonist muscle primarily controls the movement. A value approaching 1 indicates significant co-activation of both muscles, which may reflect the need for joint stability or reduced movement amplitude during the specific phase of motion.
2.10 Statistical analysis
All collected data were expressed as mean ± standard error (M ± SEM).
Trial 1. Mauchly’s test was used to assess the sphericity of the repeated-measures factor. If the sphericity assumption was violated, the Greenhouse–Geisser correction was applied to adjust the degrees of freedom. A one-way repeated-measures analysis of variance (ANOVA) was conducted to analyze the NRS scores at different time points during m-tDCS. Post-hoc pairwise comparisons were performed using the Bonferroni correction. The effect size was calculated using partial eta-squared
Trial 2. Student’s unpaired t-test was used to compare the baseline characteristics (age, body mass index, lower-limb motor capacity, HR, RPE, and ΔL) between the two intervention groups (m-tDCS and sham) to ensure comparability prior to stimulation. To evaluate the effect of m-tDCS on lower-limb motor capacity, a mixed-model ANOVA was conducted, with time (pre-stimulation and post-stimulation) as the within-subjects factor and group (m-tDCS and sham) as the between-subjects factor. The following outcome measures collected at 85%–100% PPO were analyzed: time to TTE, W, P, RPM, HR, RPE, ΔL, MRC, and CAI. When significant interaction effects were detected, Bonferroni-corrected post-hoc pairwise comparisons were performed. For main effects of time or protocol, post-hoc analysis was conducted using paired Student’s t-test (for time) or unpaired Student’s t-test (for group). If significant baseline differences were observed between groups, further analysis was performed using unpaired Student’s t-test to compare the change in variables from pre- to post-stimulation between the two groups. Effect sizes were calculated as partial eta-squared
3 Results
3.1 Testing the performance of the stimulation protocol’s tolerability and side effects
Nevertheless, the pain levels decreased progressively as the stimulation continued (see Table 2). A significant change in pain levels was observed throughout the stimulation period, with a main effect of time. No participant reported extreme pain ratings, and there were no dropouts in the experiment. A total of 73.3% of the participants experienced a tingling sensation under the electrodes during m-tDCS, with 53.3% reporting itching, 40% reporting a burning sensation, and 6.67% reporting drowsiness.
3.2 Effects of m-tDCS on motor performance
3.2.1 Comparison of the participant-based characteristics and motor capacity
All participants successfully completed the incremental exercise test, which enabled a valid comparison of motor capacity. No significant differences in participant demographics and lower-limb maximum motor capacity were found between the m-tDCS and sham groups (p > 0.05; see Table 3).
3.2.2 Cycling endurance performance at 85%–100% PPO
Participants completed the test to their maximum motor capacity. No participants reported any sports injury during cycling.
3.2.2.1 Motor performance of the lower limb
A significant stimulation protocol × time interaction was found for W (p < 0.05) (Table 4) (Figure 4). A between-group comparison at baseline showed no significant difference (F = 0.92, p = 0.342), whereas post-stimulation, W of the m-tDCS group was significantly greater than that of the sham group (F = 4.81, p = 0.034).
Figure 4. Effects of m-tDCS on lower-limb motor performance at 85%–100% PPO. TTE indicates cycling time; W indicates the work; P indicates the power output, RPM indicates revolutions per minute. # indicates the significant differences between groups: p < 0.05.
A significant main effect of the stimulation protocol was observed for TTE (p < 0.05) (Table 4) (Figure 4). TTE was significantly higher in the m-tDCS group than in the sham group at both pre- (t = −2.26, p = 0.029) and post-stimulation (t = −2.82, p = 0.007). However, there was no significant change in TTE between groups (t = −1.92, p = 0.062).
A significant main effect of the stimulation protocol was also found for RPM (p < 0.05) (Table 4) (Figure 4). At baseline, no difference in RPM was recorded between m-tDCS and sham groups (t = −1.93, p = 0.061); however, after stimulation, RPM was significantly higher in the m-tDCS group than in the sham group (t = −2.22, p = 0.032).
No significant main effects of stimulation protocol or stimulation time, nor their interaction, were observed for P (p > 0.05).
3.2.2.2 Physiological performance
A significant main effect of stimulation time was found for HR (p < 0.05) (Table 3) (Figure 5). Post-hoc analysis indicated that HR was significantly reduced after m-tDCS (t = 3.65, p = 0.001), whereas no significant change was observed in the sham group (t = 1.41, p = 0.172). No significant main effects of stimulation protocol or stimulation time, nor their interaction, were detected for ΔL (p > 0.05).
Figure 5. Effects of m-tDCS on the physiological performance at 85%–100% PPO. HR indicates the heart rate; ΔL indicates the blood lactate level. * indicates the significant differences between pre- and post-intervention: p < 0.05.
3.2.2.3 Psychological performance
Analysis of RPE showed a significant main effect of stimulation time (p < 0.05) (Table 3) (Figure 6). Post-hoc analysis revealed that RPE significantly increased after sham stimulation (t = −2.09, p = 0.049), whereas no significant change was observed following tDCS (t = −1.39, p = 0.180).
Figure 6. Effects of m-tDCS on the psychological performance at 85%–100% PPO. RPE indicates the rate of perception effort. * indicates the significant differences between pre- and post-intervention: p < 0.05.
3.2.3 Contribution of the dominant lower-limb muscle before and after stimulation
3.2.3.1 Performance of the lower-limb muscle in the propulsion phase
Figures 7, 8 show the overall view of MCR performance in the propulsion phase. The interaction between stimulation time and stimulation protocol significantly affected RF-MCR (p < 0.05). Within-group analysis revealed that MCR significantly decreased after sham stimulation (t = 4.38, p < 0.001), whereas no significant change was observed following m-tDCS (t = 0.67, p = 0.513). Between-group analysis showed no significant difference in MCR between m-tDCS and sham groups at baseline (t = 0.64, p = 0.528); however, after stimulation, MCR of the m-tDCS group was significantly higher than that of the sham group (t = −2.59, p = 0.013). The interaction between stimulation time and protocol was also significant for ST-MCR (p < 0.05). Within-group analysis showed a significant increase in MCR following sham stimulation (t = −1.95, p = 0.066), whereas no significant change occurred after m-tDCS (t = 0.67, p = 0.513). Between-group analysis revealed no significant difference between m-tDCS and sham groups at baseline (t = −0.84, p = 0.404), but MCR of the m-tDCS group was significantly lower than that in the sham group after stimulation (t = 5.77, p < 0.001). For TA, the interaction between the stimulation time and protocol was also significant (p < 0.05). Within-group analysis revealed a significant decrease in MCR after sham stimulation (t = 3.83, p = 0.001), whereas no significant change occurred after m-tDCS (t = 0.18, p = 0.861). Between-group analysis showed no significant difference in MCR between m-tDCS and sham groups at baseline (t = 0.38, p = 0.704), but MCR of the m-tDCS group was significantly higher than that of the sham group after stimulation (t = −2.15, p = 0.038).
Figure 7. Overall view of MCR changes at different interventions during the propulsion phase at 85%–100% PPO. MCR indicates muscle contribution ratio; RF indicates rectus femoris; VL indicates vastus lateralis; VM indicates vastus medialis; BF indicates biceps femoris, long head; ST indicates semitendinosus; GM indicates medial gastrocnemius; GL indicates lateral gastrocnemius. (a) MCR changes in the sham group. (b) MCR changes in the m-tDCS group.
Figure 8. Performance of the muscle contribution ratio of the lower-limb muscle in the propulsion phase. MCR indicates the muscle contribution ratio; RF indicates rectus femoris; VL indicates vastus lateralis; VM indicates vastus medialis; BF indicates biceps femoris, long head; ST indicates semitendinosus; GM indicates medial gastrocnemius; GL indicates lateral gastrocnemius. The asterisk indicates the significant difference between pre- and post-intervention: * represents p < 0.05, ** represents p < 0.01, and *** represents p < 0.001; the well number indicates the significant difference between groups: # represents p < 0.05, and ### represents p < 0.001.
There was a significant main effect of the stimulation protocol on BF-MCR (p < 0.05). No significant difference was found in MCR between m-tDCS and sham groups at baseline (t = 1.17, p = 0.247), but MCR of the m-tDCS group was significantly lower than that of the sham group after stimulation (t = 2.14, p = 0.039). A significant main effect of the stimulation protocol was also observed for VM-MCR (p < 0.05). The MCR of the m-tDCS group was significantly higher than that of the sham group at both baseline (t = −2.44, p = 0.019) and post-stimulation (t = −3.40, p = 0.02). Further analysis of the change in MCR showed no significant difference between the two intervention groups (t = −0.80, p = 0.428).
No significant interaction between stimulation time and protocol was found for GL-, GM-, or VL-MCR (p > 0.05), and no significant main effects of time or protocol were observed for these muscle groups (p > 0.05).
3.2.3.2 Performance of the lower-limb muscle in the pull phase
Figures 9, 10 show the overall view of MCR performance in the pull phase. The interaction between stimulation protocol and stimulation time had a significant effect on RF-MCR (p < 0.05) (Figure 9). Within-group analysis revealed that RF-MCR was significantly increased after sham stimulation (t = −2.09, p = 0.040), whereas it significantly decreased following m-tDCS (t = 3.59, p = 0.01). Between-group analysis showed no significant difference in RF-MCR between m-tDCS and sham groups at baseline (t = 0.35, p = 0.726), but MCR was significantly lower in the m-tDCS group than in the sham group after stimulation (t = 3.90, p < 0.01). The interaction between stimulation protocol and stimulation time was also significant for the VM (p < 0.05). Within-group analysis showed that VM-MCR significantly decreased after sham stimulation (t = 2.78, p = 0.007), whereas no significant change was observed after m-tDCS (t = −0.71, p = 0.479). Between-group analysis indicated that VM-MCR was significantly higher in the m-tDCS group than in the sham group at both baseline (t = −4.18, p < 0.001) and post-stimulation (t = −8.67, p < 0.001). Further analysis of the change in MCR revealed that the reduction in MCR after sham stimulation was significantly greater than that after m-tDCS (t = −2.65, p = 0.009).
Figure 9. Overall view of MCR changes at different interventions during the pull phase at 85%–100% PPO. MCR indicates the muscle contribution ratio; RF indicates rectus femoris; VL indicates vastus lateralis; VM indicates vastus medialis; BF indicates biceps femoris, long head; ST indicates semitendinosus; GM indicates medial gastrocnemius; GL indicates lateral gastrocnemius. (a) MCR changes in the sham group. (b) MCR changes in the m-tDCS group.
Figure 10. Performance of the muscle contribution ratio of the lower-limb muscle in the pull phase. MCR indicates muscle contribution ratio; RF is rectus femoris; VL indicates vastus lateralis; VM indicates vastus medialis; BF indicates biceps femoris, long head; ST indicates semitendinosus; GM indicates medial gastrocnemius; GL indicates lateral gastrocnemius. The asterisk indicates the significant difference between pre- and post-intervention: * represents p < 0.05, ** represents p < 0.01, and *** represents p < 0.001; the well number indicates the significant difference between groups: # represents p < 0.05, and ### represents p < 0.001.
A significant main effect of the stimulation protocol was found for the BF (p < 0.05). BF-MCR was significantly higher in the m-tDCS group than in the sham group at both baseline (t = 3.87, p < 0.001) and post-stimulation (t = 5.63, p < 0.001). However, there was no significant difference in MCR between m-tDCS and sham groups (t = 1.78, p = 0.076). A significant main effect of stimulation protocol was also observed for the GL (p < 0.05). At baseline, there was no significant difference in MCR between m-tDCS and sham groups (t = 1.40, p = 0.164). However, after stimulation, MCR was significantly lower in the m-tDCS group than in the sham group (t = 2.83, p = 0.005).
Significant main effects of time were observed for VL-, ST-, TA-, GM-, and GL-MCR (p < 0.05). Specifically, within-group analysis showed that VL-MCR significantly decreased after sham stimulation (t = 2.67, p = 0.009), but no significant change was found after m-tDCS (t = 0.30, p = 0.762). The ST-MCR increased significantly after sham stimulation (t = −2.00, p = 0.049), whereas no significant change was found following m-tDCS (t = −1.15, p = 0.255). Both sham and m-tDCS groups showed significant reductions in TA-MCR (sham: t = 2.30, p = 0.024; m-tDCS: t = 2.20, p = 0.031). The GM-MCR increased significantly in both sham and m-tDCS groups (sham: t = −2.96, p = 0.004; m-tDCS: t = −4.24, p < 0.001). Finally, the GL-MCR increased significantly after sham stimulation (t = −3.44, p = 0.001), but no significant change was observed following m-tDCS (t = −1.78, p = 0.079).
3.2.4 Performance comparison of antagonist and agonist muscles in the knee joint before and after stimulation
Based on the characteristics of the movement cycle, the muscles monitored in this study were categorized as follows: during the propulsion phase, the RF, VL, and VM were classified as agonist muscles, whereas the BF, ST, GM, and GL acted as antagonist muscles. Analysis of the CAI for active muscles during the propulsion phase revealed no significant interaction effects between m-tDCS and sham groups (p > 0.05), and there were no significant main effects of stimulation protocol or stimulation time (p > 0.05) (Figure 11).
Figure 11. Co-activation index of agonist muscles around the knee at 85%–100% PPO. The asterisk indicates the significant difference between pre- and post-intervention: *** represents p < 0.001; the well number indicates the significant difference between groups: ## represents p < 0.01.
In contrast to the propulsion phase, the BF, ST, GM, and GL acted as agonist muscles, whereas the RF, VL, and VM were classified as antagonist muscles during the pull phase. Analysis of the CAI for agonist muscles during the pull phase revealed a significant interaction effect between stimulation time and protocol (p < 0.05). Within-group analysis showed a significant increase in CAI after sham stimulation compared to baseline (t = −5.49, p < 0.001), whereas no significant difference was observed before and after m-tDCS (t = −1.45, p = 0.152). Between-group analysis revealed that the CAI in the m-tDCS group was significantly higher than that in the sham group at baseline (t = −2.97, p = 0.003), but no significant difference was found between the two groups after stimulation (t = 0.32, p = 0.748). Further analysis of the changes in CAI between sham stimulation and m-tDCS revealed that the increase in CAI in the sham group was significantly greater than that in the m-tDCS group (t = 3.61, p < 0.001).
4 Discussion
This study aimed to evaluate the effects of m-tDCS targeting lower-limb-specific cortical areas on cycling endurance motor performance in healthy adults. The findings provide compelling evidence that m-tDCS is a safe and non-invasive brain stimulation technique that can significantly enhance lower-limb endurance capacity. This improvement is accompanied by changes in lower-limb muscle activity and coordination, suggesting that m-tDCS optimizes endurance performance by modulating cortical excitability. Additionally, significant differences in the RPM were observed during the high-intensity phase, indicating enhanced motor control and coordination post-stimulation.
The multifocal montage is designed to engage a distributed motor network, potentially increasing focality at cortical targets and improving current steering relative to single-site approaches (Dmochowski et al., 2011). Within this framework, the present propulsion-specific maintenance of quadriceps MCR and the pull-phase restraint of antagonistic co-activation (stable CAI) align with enhanced neural efficiency—i.e., a more economical recruitment of task-relevant synergies and suppression of non-task-relevant activity (Falconer and Winter, 1985; Li et al., 2020). Because P and ΔL were unchanged, the performance pattern (higher W at post-stimulation vs. sham and higher RPM) is better explained by coordination-driven endurance support than by increased maximal power output or altered metabolic accumulation. Such synergy-level modulation aligns with prior reports that tDCS can facilitate endurance-type outcomes without necessarily increasing peak strength (Vitor-Costa et al., 2015; Angius et al., 2016; Angius et al., 2018; Alix-Fages et al., 2019; Chinzara et al., 2022).
A substantial body of work indicates that anodal or extracephalic tDCS may extend time-to-task failure and attenuate perceived exertion under certain protocols (Vitor-Costa et al., 2015; Angius et al., 2016; Angius et al., 2018). Meta-analytic summaries suggest small-to-moderate, parameter-dependent effects on endurance with mixed evidence for strength and power (Alix-Fages et al., 2019; Chinzara et al., 2022). The current findings refine this picture in two ways: (i) by using a multifocal montage that likely acts on a broader motor network (Dmochowski et al., 2011) and (ii) by linking behavioral outcomes to phase-dependent adjustments in MCR/CAI. Rather than a generalized boost in the maximal output, our results highlight preservation of terminal-phase work and cadence via coordination tuning—a distinction with practical relevance for a prolonged or fatiguing task. The significant changes in MCR and CAI observed in this study provide objective, physiological evidence of how m-tDCS modulates muscle activation patterns to enhance endurance performance.
Compared with previous studies, this study also contributes to the growing body of literature on tDCS by providing new insights into its effects on lower-limb motor performance, particularly endurance (Vitor-Costa et al., 2015; Angius et al., 2016; Angius et al., 2018; Alix-Fages et al., 2019; Chinzara et al., 2022). The changes in muscle coordination, as reflected by MCR and CAI, offer a novel mechanism through which m-tDCS enhances motor performance (Falconer and Winter, 1985; Li et al., 2020; Dmochowski et al., 2011). These findings support the idea that tDCS, especially multifocal stimulation, can improve neural efficiency, muscle activation, and coordination, leading to enhanced endurance (Dmochowski et al., 2011; Angius et al., 2018; Chinzara et al., 2022; Pilloni et al., 2020). From a theoretical standpoint, these results deepen our understanding of how brain stimulation affects the neural networks responsible for motor control and endurance (Dmochowski et al., 2011; Chinzara et al., 2022). The modulation of muscle coordination and the efficiency of muscle activation by m-tDCS could represent a key mechanism underlying the observed endurance improvements (Dmochowski et al., 2011; Falconer and Winter, 1985; Li et al., 2020; Alix-Fages et al., 2019).
Practically, these findings have significant implications for both sports performance and rehabilitation. The ability of m-tDCS to improve endurance performance and muscle coordination makes it a valuable tool for athletes in endurance sports. By enhancing neuromuscular coordination and efficiency, m-tDCS could help athletes optimize performance during prolonged physical exertion. In rehabilitation, m-tDCS could be used to assist individuals with motor impairments in improving muscle coordination and endurance (Pilloni et al., 2020; Li et al., 2020). Given the observed changes in MCR and CAI, m-tDCS may help enhance neuromuscular function in clinical populations, particularly those recovering from neurological disorders such as stroke or multiple sclerosis (Pilloni et al., 2020; Li et al., 2020). Therefore, based on our findings and from an applied standpoint, we suggest that m-tDCS is most likely to be effective when paired with task engagement. Specifically, in athletes, a single pre-task m-tDCS session is feasible to deliver immediately before high-intensity cycling as a practical use case, which may help preserve late-stage endurance by supporting phase-specific neuromuscular coordination. In rehabilitation, brief pre-therapy m-tDCS paired with task-specific practice may improve tolerance and neuromuscular efficiency (Pilloni et al., 2020; Li et al., 2020). Importantly, tDCS after-effects begin to develop during stimulation and usually persist for approximately 30 min–60 min—sometimes up to 90 min—showing a short onset and gradual decline. Consequently, more durable improvements likely require multi-session protocols.
Despite the valuable insights provided by this study, several limitations should be considered. The sample size was relatively small, and the study focused only on healthy young adults. Future research should involve larger and more diverse populations, including older adults and individuals with neurological conditions, to assess the generalizability of these findings. Moreover, although this study focused on endurance performance, future research should explore the effects of m-tDCS on other motor abilities, such as strength and coordination, to gain a broader understanding of its effects on motor function (Alix-Fages et al., 2019; Chinzara et al., 2022). Additionally, this study did not include long-term follow-up to assess the durability of the observed effects. Future studies should incorporate follow-up assessments to determine whether the performance improvements persist over time (Pilloni et al., 2020). Furthermore, the neural mechanisms underlying the effects of m-tDCS on muscle coordination remain to be elucidated. Although the results suggest changes in muscle coordination via increased MCR and reduced CAI, further investigation using neuroimaging techniques such as functional MRI or electroencephalography could provide more detailed insights into the brain activity changes that accompany these improvements in motor performance (Dmochowski et al., 2011). Finally, it is theorized that m-tDCS offers improved focality and intensity compared to the classical single-site tDCS approach; however, we did not conduct a direct comparison between these two techniques in our study. This comparison could further clarify the distinct advantages and potential differences in their effects on performance. We acknowledge the importance of this issue, and thus future research should explore and compare the effects of both techniques to provide a more comprehensive understanding of their relative efficacy.
5 Conclusion
m-tDCS, delivered with a multifocal montage, did not increase peak mechanical capacity but helped preserve late-stage work and cadence during high-intensity cycling, accompanied by phase-specific optimization of muscle coordination—maintenance of propulsion-relevant contributions and pull-phase restraint of antagonistic co-activation. These findings refine the view of the effects of tDCS on endurance: coordination-centric, network-level efficiency rather than a uniform boost in the output of motor capacity.
Data availability statement
The original contributions presented in the study are included in the article/supplementary material; further inquiries can be directed to the corresponding authors.
Ethics statement
The studies involving humans were approved by the Ethics Committee of the Shanghai University of Sport. The studies were conducted in accordance with the local legislation and institutional requirements. Written informed consent for participation in this study was provided by the participants’ legal guardians/next of kin.
Author contributions
ZL: Conceptualization, Funding acquisition, Methodology, Project administration, Writing – original draft, Writing – review and editing. PR: Writing – review and editing. XG: Writing – review and editing. YF: Funding acquisition, Project administration, Writing – review and editing. AA: Supervision, Writing – review and editing. YL: Funding acquisition, Project administration, Supervision, Writing – original draft, Writing – review and editing.
Funding
The authors declare that financial support was received for the research and/or publication of this article. This study was supported by the China Postdoctoral Science Foundation (2024M761525), Ningbo Clinical Research Center for Orthopedics and Exercise Rehabilitation (2024L004), Research Academy of Medicine Combining Sports (2023001), Zhejiang Province Rehabilitation Medical Society Research Special Fund (ZKKY2024012), and No. 2 Hospital of Ningbo City HwaMei Research Fund (2024HMKYA53).
Acknowledgements
The authors thank all participants for their commitment and the staff of the Human Performance Laboratory for assistance with participant screening, data acquisition, and stimulation setup.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Generative AI statement
The authors declare that Generative AI was used in the creation of this manuscript. Generative AI was used solely for grammar and language editing of the manuscript (ChatGPT and OpenAI GPT-4). No text, data analysis, images, references, or scientific content were generated by AI. All authors reviewed the edits and take full responsibility.
Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
References
Alix-Fages C., Romero-Arenas S., Castro-Alonso M., Colomer-Poveda D., Río-Rodriguez D., Jerez-Martínez A., et al. (2019). Short-term effects of anodal transcranial direct current stimulation on endurance and maximal force production: a systematic review and meta-analysis. J. Clin. Med. 8, 536. doi:10.3390/jcm8040536
Angius L., Pageaux B., Hopker J. G., Marcora S. M., Mauger A. R. (2016). Transcranial direct current stimulation improves isometric time to exhaustion of the knee extensors. Neuroscience 339, 363–375. doi:10.1016/j.neuroscience.2016.10.028
Angius L., Mauger A. R., Hopker J., Pascual-Leone A., Santarnecchi E., Marcora S. M. (2018). Bilateral extracephalic transcranial direct current stimulation improves endurance performance in healthy individuals. Brain Stimul. 11 (1), 108–117. doi:10.1016/j.brs.2017.09.017
Chinzara T. T., Buckingham G., Harris D. J. (2022). Transcranial direct current stimulation and sporting performance: a systematic review and meta-analysis of transcranial direct current stimulation effects on physical endurance, muscular strength and visuomotor skills. Eur. J. Neurosci. 55, 468–486. doi:10.1111/ejn.15540
De Wachter J., Van Cutsem J., Pattyn N., De Pauw K., Roelands B. (2021). Prefrontal cortex oxygenation during endurance performance: a systematic review. Front. Physiology 12, 761232. doi:10.3389/fphys.2021.761232
Dmochowski J. P., Datta A., Bikson M., Su Y., Parra L. C. (2011). Optimized multi-electrode stimulation increases focality and intensity at target. J. Neural Eng. 8 (4), 046011. doi:10.1088/1741-2560/8/4/046011
Falconer K., Winter D. A. (1985). Quantitative assessment of co-contraction at the ankle joint in walking. Electromyogr. Clin. Neurophysiology 25 (2–3), 135–149.
Flor G. Z., Monteiro W., da Silva R. C. F., de Oliveira B. R. R., Marquez G., Lattari E. (2025). Acute effects of anodal transcranial direct current stimulation on endurance and maximal voluntary contraction in lower limbs: a systematic review and meta-analysis. Exp. Brain Res. 243 (3), 57. doi:10.1007/s00221-025-07008-9
Hallett M. (2007). Transcranial magnetic stimulation: a primer. Neuron 55 (2), 187–199. doi:10.1016/j.neuron.2007.06.026
Holgado D., Vadillo M. A., Sanabria D. (2019). The effects of transcranial direct current stimulation on objective and subjective indexes of exercise performance: a systematic review and meta-analysis. Brain Stimul. 12 (2), 242–250. doi:10.1016/j.brs.2018.12.002
Holgado D., Sanabria D., Vadillo M. A., Román-Caballero R. (2024). Zapping the brain to enhance sport performance? An umbrella review of the effect of transcranial direct current stimulation on physical performance. Neurosci. and Biobehav. Rev. 164, 105821. doi:10.1016/j.neubiorev.2024.105821
Huang Y., Parra L. C. (2019). Can transcranial electric stimulation with multiple electrodes reach deep targets? Brain Stimul. 12, 30–40. doi:10.1016/j.brs.2018.09.010
Hug F., Dorel S. (2009). Electromyographic analysis of pedaling: a review. J. Electromyogr. Kinesiol. 19 (2), 182–198. doi:10.1016/j.jelekin.2007.10.010
Hug F., Turpin N. A., Couturier A., Dorel S. (2011). Consistency of muscle synergies during pedaling across different mechanical constraints. J. Neurophysiology 106 (1), 91–103. doi:10.1152/jn.01096.2010
Jeukendrup A., Saris W. H., Brouns F., Kester A. D. (1996) A new validated endurance performance test. Med. Sci. Sports Exerc., 28(2): 266–270. doi:10.1097/00005768-199602000-00017
Kaushalya F. S., Romero-Arenas S., García-Ramos A., Colomer-Poveda D., Marquez G. (2022). Acute effects of transcranial direct current stimulation on cycling and running performance: a systematic review and meta-analysis. Eur. J. Sport Sci. 22 (2), 113–125. doi:10.1080/17461391.2020.1856933
Lefaucheur J.-P., Aleman A., Baeken C., Benninger D. H., Brunelin J., Di Lazzaro V., et al. (2020). Evidence-based guidelines on the therapeutic use of repetitive transcranial magnetic stimulation (rTMS): an update (2014–2018). Clin. Neurophysiol. 131 (2), 474–528. doi:10.1016/j.clinph.2019.11.002
Li G., Shourijeh M. S., Ao D., Patten C., Fregly B. J. (2020). How well do commonly used co-contraction indices approximate lower-limb joint stiffness trends during gait? Front. Bioeng. Biotechnol. 8, 588908. doi:10.3389/fbioe.2020.588908
Liang Z. (2024). Functional near-infrared spectroscopy study of hemodynamic in the prefrontal and motor cortices and its implications for endurance capacity. Mol. Cell. Biomechanics 21 (3), 350. doi:10.62617/mcb350
Liang Z., Xue G., Ramasawmy P. (2024). Multifocal tDCS with anodal configuration over motor cortex induce higher excitability than cathodal configuration. Med. Sci. Sports Exerc. 56 (10s), 725–726. doi:10.1249/01.mss.0001058872.18386.a1
Machado D. G., Unal G., Andrade S. M., Moreira A., Altimari L. R., Brunoni A. R., et al. (2019). Effect of transcranial direct current stimulation on exercise performance: a systematic review and meta-analysis. Brain Stimul. 12 (3), 593–605. doi:10.1016/j.brs.2018.12.227
Neptune R. R., Kautz S. A., Hull M. L. (1997). The effect of pedaling rate on coordination in cycling. J. Biomechanics 30 (10), 1051–1058. doi:10.1016/S0021-9290(97)00071-7
Nitsche M. A., Paulus W. (2000). Excitability changes induced in the human motor cortex by weak transcranial direct current stimulation. J. Physiology 527 (3), 633–639. doi:10.1111/j.1469-7793.2000.t01-1-00633.x
Nitsche M. A., Paulus W. (2001). Sustained excitability elevations induced by transcranial DC motor cortex stimulation in humans. Neurology 57 (10), 1899–1901. doi:10.1212/WNL.57.10.1899
Okano A. H., Fontes E. B., Montenegro R. A., Farinatti P. d. T. V., Cyrino E. S., Li L. M., et al. (2015). Brain stimulation modulates the autonomic nervous system, rating of perceived exertion and performance during maximal exercise. Br. J. Sports Med. 49 (18), 1213–1218. doi:10.1136/bjsports-2012-091658
O’Bryan S. J., Brown N. A., Billaut F., Rouffet D. M. (2014). Changes in muscle coordination and power output during sprint cycling. Neurosci. Lett. 576, 11–16. doi:10.1016/j.neulet.2014.05.023
Pilloni G., Choi C., Coghe G., Krupp L., Moffat M., et al. (2020). Walking in multiple sclerosis improves with tDCS: a randomized, double-blind, sham-controlled study. Ann. Clin. Transl. Neurology 7 (11), 2310–2319. doi:10.1002/acn3.51224
Polanía R., Nitsche M. A., Ruff C. C. (2018). Studying and modifying brain function with non-invasive brain stimulation. Nat. Neurosci. 21 (2), 174–187. doi:10.1038/s41593-017-0054-4
Qi S., Yu J., Li L., Dong C., Ji Z., Cao L., et al. (2024). Advances in non-invasive brain stimulation: enhancing sports performance function and insights into exercise science. Front. Hum. Neurosci. 18, 1477111. doi:10.3389/fnhum.2024.1477111
Ruffini G., Fox M. D., Ripolles O., Miranda P. C., Pascual-Leone A. (2014). Optimization of multifocal transcranial current stimulation for weighted cortical pattern targeting from realistic modeling of electric fields. Neuroimage 89, 216–225. doi:10.1016/j.neuroimage.2013.12.002
Silvanto J., Muggleton N., Walsh V. (2008). State-dependency in brain stimulation studies of perception and cognition. Trends Cognitive Sci. 12 (12), 447–454. doi:10.1016/j.tics.2008.09.004
Keywords: multifocal transcranial direct current stimulation, lower-limb cortical networks, incremental cycling test, endurance performance, muscle coordination
Citation: Liang Z, Ramasawmy P, Guo X, Fang Y, Antal A and Liu Y (2025) Multifocal tDCS targeting lower-limb cortical areas preserves late-stage endurance and tunes phase-specific coordination during incremental cycling. Front. Physiol. 16:1706478. doi: 10.3389/fphys.2025.1706478
Received: 30 September 2025; Accepted: 07 November 2025;
Published: 02 December 2025.
Edited by:
Daniel Gomes Da Silva Machado, Federal University of Rio Grande do Norte, BrazilReviewed by:
Susannah L. Williamson, I Corps at U.S. Army, United StatesFederica Marchiotto, University of Verona, Italy
Copyright © 2025 Liang, Ramasawmy, Guo, Fang, Antal and Liu. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Yufei Fang, ZnlmMzY5QDEyNi5jb20=; Andrea Antal, YWFudGFsQGd3ZGcuZGU=; Yu Liu, eXVsaXVAc3VzLmVkdS5jbg==