Abstract
Repetitive transcranial magnetic stimulation (rTMS) is a promising intervention tool for the noninvasive modulation of brain activity and behavior in neuroscience research and clinical settings. However, the resting-state dynamic evolution of large-scale functional brain networks following rTMS has rarely been investigated. Here, using resting-state fMRI images collected from 23 healthy individuals before (baseline) and after 1 Hz rTMS of the left frontal (FRO) and occipital (OCC) lobes, we examined the different effects of rTMS on brain dynamics across the human cortex. By fitting a pairwise maximum entropy model (pMEM), we constructed an energy landscape for the baseline and poststimulus conditions by fitting a pMEM. We defined dominant brain states (local minima) in the energy landscape with synergistic activation and deactivation patterns of large-scale functional networks. We calculated state dynamics including appearance probability, transitions and duration. The results showed that 1 Hz rTMS induced increased and decreased state probability, transitions and duration when delivered to the FRO and OCC targets, respectively. Most importantly, the shortest path and minimum cost between dominant brain states were altered after stimulation. The absolute sum of the costs from the source states to the destinations was lower after OCC stimulation than after FRO stimulation. In conclusion, our study characterized the dynamic trajectory of state transitions in the energy landscape and suggested that local rTMS can induce significant dynamic perturbation involving stimulated and distant functional networks, which aligns with the modern view of the dynamic and complex brain. Our results suggest low-dimensional mapping of rTMS-induced brain adaption, which will contribute to a broader and more effective application of rTMS in clinical settings.