Macroscopic resting state model predicts theta burst stimulation response: A randomized trial

Repetitive transcranial magnetic stimulation (rTMS) is a promising alternative therapy for treatment-resistant depression, although its limited remission rate indicates room for improvement. As depression is a phenomenological construction, the biological heterogeneity within this syndrome needs to be considered to improve the existing therapies. Whole-brain modeling provides an integrative multi-modal framework for capturing disease heterogeneity in a holistic manner. Computational modelling combined with probabilistic nonparametric fitting was applied to the resting-state fMRI data from 42 patients (21 women), to parametrize baseline brain dynamics in depression. All patients were randomly assigned to two treatment groups, namely active (i.e., rTMS, n = 22) or sham (n = 20). The active treatment group received rTMS treatment with an accelerated intermittent theta burst protocol over the dorsomedial prefrontal cortex. The sham treatment group underwent the identical procedure but with the magnetically shielded side of the coil. We stratified the depression sample into distinct covert subtypes based on their baseline attractor dynamics captured by different model parameters. Notably, the two detected depression subtypes exhibited different phenotypic behaviors at baseline. Our stratification could predict the diverse response to the active treatment that could not be explained by the sham treatment. Critically, we further found that one group exhibited more distinct improvement in certain affective and negative symptoms. The subgroup of patients with higher responsiveness to treatment exhibited blunted frequency dynamics for intrinsic activity at baseline, as indexed by lower global metastability and synchrony. Our findings suggested that whole-brain modeling of intrinsic dynamics may constitute a determinant for stratifying patients into treatment groups and bringing us closer towards precision medicine.