Skip to content
DraconiticResearch
Search Sources API
201,683 studies31 sourcesUpdated 3 min ago
PreprintCC BY-NC-NDmedRxiv

A whole-brain modeling framework for tDCS montage optimization in drug-resistant epilepsy

  1. 1Neuroelectrics, Barcelona, Spain

Abstract

ObjectiveTranscranial direct current stimulation (tDCS) is a promising non-invasive therapy for drug-resistant epilepsy, but the montage optimization pipelines used to personalize it often rely on biophysical head models that capture only the induced electric field and ignore the network dynamics through which stimulation influences epileptogenic activity. We aimed to develop and evaluate a montage optimization framework that incorporates individualized network dynamics.

ApproachWe introduce a framework based on a digital twin, or neurotwin, a personalized replica of a patients brain that couples a biophysical head model of the stimulation-induced field with a personalized whole-brain model (WBM) of seizure dynamics constrained by structural connectivity and intracranial recordings. We applied it to produce optimal montages in 12 patients with drug-resistant epilepsy, compared them with conventional biophysics-based optimization, and assessed the framework retrospectively against the clinical outcomes of six treated patients.

Main resultsBecause the neurotwin optimization suppresses simulated spread by construction, the informative result is the divergence between pipelines: the network-informed and field-based optimizations produced substantially different montages despite delivering comparable inhibitory fields at the epileptogenic zone, and montages with nearly identical field distributions could yield markedly different network outcomes within the model. A minimum-replacement-set analysis traced these differences to field changes at a few propagation-zone and off-target parcels rather than to the overall field. The parcels whose inhibition most reduced spread tended to occupy highly connected (hub) positions in the structural connectome, a robust but partial association indicating that network structure shapes, but does not fully determine, where stimulation is most effective. In the retrospective analysis, model-predicted spread reduction tracked clinical seizure-frequency change in the expected direction, though not significantly (n = 6).

SignificanceWe hypothesize that accounting for individualized network dynamics may improve seizure control, a hypothesis now being tested prospectively (NCT06334952).

HighlightsO_LIA neurotwin couples a biophysical head model with a personalized whole-brain model (WBM) for tDCS optimization. C_LIO_LIneurotwin and field-based optimization yield divergent montages despite comparable EZ fields. C_LIO_LIField differences at a few target and off-target parcels, not overall field magnitude, drive the simulated network outcomes. C_LIO_LIParcels controlling simulated seizure spread tend to be structural network hubs (high nodal strength and centrality), an association that persists after adjusting for local dynamics. C_LIO_LIWBM-predicted spread reduction was directionally consistent with clinical seizure-frequency change (non-significant, n = 6). C_LI

Related

  1. Epilepsy Surgery vs Medical Management for Pediatric Drug-Resistant Focal Epilepsy
    T. Abel, E. Harford, D. A. Silliman +3
    medRxivPreprintAll rights reservedFull textNeurologymedRxiv
  2. Factors associated with drug-resistant epilepsy in children with epilepsy: a retrospective cohort study
    Qunyan Su, Yandan Yin, Haiying Cheng +8
    Translational pediatricsPeer-reviewedCC BY-NC-NDFull textGeneral medicineEurope PMC
  3. Diffusion MRI abnormalities localise the epileptogenic zone in paediatric drug-resistant epilepsy
    D. Veljanoski, A. Chari, K. K. Seunarine +7
    medRxivPreprintAll rights reservedFull textSurgery & anaesthesiamedRxiv
  4. Serum metabolomics identifies ganglioside dysregulation and pyruvate dehydrogenase complex–related metabolic disturbances in pediatric drug-resistant epilepsy
    Eman Y. Abu-Rish, Reem H. AlMalki, Amira Masri +4
    Frontiers in pediatricsPeer-reviewedCC BYFull textGeneral medicineEurope PMC
  5. Awake resection of temporal encephalocele associated with drug-resistant epilepsy: illustrative case
    Paulína Ondrušová, Gabriela Timárová, Eva Pakosová +9
    Journal of neurosurgery. Case lessonsPeer-reviewedCC BY-NC-NDFull textGeneral medicineEurope PMC
BibTeX
@misc{mercadal2026wholebrain,
  title = {{A whole-brain modeling framework for tDCS montage optimization in drug-resistant epilepsy}},
  author = {B. Mercadal and E. Lopez-Sola and J. Makhalova and F. Pizzo and R. Salvador and F. Bartolomei and G. Ruffini},
  howpublished = {medRxiv},
  year = {2026},
  month = sep,
  doi = {10.64898/2026.09.22.26363221},
  url = {https://www.medrxiv.org/content/10.64898/2026.09.22.26363221v1},
  note = {CC BY-NC-ND},
}
RIS
TY  - UNPB
TI  - A whole-brain modeling framework for tDCS montage optimization in drug-resistant epilepsy
AU  - B. Mercadal
AU  - E. Lopez-Sola
AU  - J. Makhalova
AU  - F. Pizzo
AU  - R. Salvador
AU  - F. Bartolomei
AU  - G. Ruffini
PY  - 2026
DA  - 2026/09/22
PB  - medRxiv
DO  - 10.64898/2026.09.22.26363221
UR  - https://www.medrxiv.org/content/10.64898/2026.09.22.26363221v1
AB  - ObjectiveTranscranial direct current stimulation (tDCS) is a promising non-invasive therapy for drug-resistant epilepsy, but the montage optimization pipelines used to personalize it often rely on biophysical head models that capture only the induced electric field and ignore the network dynamics through which stimulation influences epileptogenic activity. We aimed to develop and evaluate a montage optimization framework that incorporates individualized network dynamics.  ApproachWe introduce a framework based on a digital twin, or neurotwin, a personalized replica of a patients brain that couples a biophysical head model of the stimulation-induced field with a personalized whole-brain model (WBM) of seizure dynamics constrained by structural connectivity and intracranial recordings. We applied it to produce optimal montages in 12 patients with drug-resistant epilepsy, compared them with conventional biophysics-based optimization, and assessed the framework retrospectively against the clinical outcomes of six treated patients.  Main resultsBecause the neurotwin optimization suppresses simulated spread by construction, the informative result is the divergence between pipelines: the network-informed and field-based optimizations produced substantially different montages despite delivering comparable inhibitory fields at the epileptogenic zone, and montages with nearly identical field distributions could yield markedly different network outcomes within the model. A minimum-replacement-set analysis traced these differences to field changes at a few propagation-zone and off-target parcels rather than to the overall field. The parcels whose inhibition most reduced spread tended to occupy highly connected (hub) positions in the structural connectome, a robust but partial association indicating that network structure shapes, but does not fully determine, where stimulation is most effective. In the retrospective analysis, model-predicted spread reduction tracked clinical seizure-frequency change in the expected direction, though not significantly (n = 6).  SignificanceWe hypothesize that accounting for individualized network dynamics may improve seizure control, a hypothesis now being tested prospectively (NCT06334952).  HighlightsO_LIA neurotwin couples a biophysical head model with a personalized whole-brain model (WBM) for tDCS optimization. C_LIO_LIneurotwin and field-based optimization yield divergent montages despite comparable EZ fields. C_LIO_LIField differences at a few target and off-target parcels, not overall field magnitude, drive the simulated network outcomes. C_LIO_LIParcels controlling simulated seizure spread tend to be structural network hubs (high nodal strength and centrality), an association that persists after adjusting for local dynamics. C_LIO_LIWBM-predicted spread reduction was directionally consistent with clinical seizure-frequency change (non-significant, n = 6). C_LI
ER  -