TY - JOUR
T1 - Personalized whole-brain Ising models with heterogeneous nodes capture differences among brain regions
AU - Craig, Adam
AU - Chen, Sida
AU - Tang, Qianyuan
AU - Zhou, Changsong
N1 - This work was partially supported by STI 2030-Major Projects (No. 2022ZD0208500), the Hong Kong Research Grant Council Senior Research Fellow Scheme (SRFS2324-2S05) and General Competitive Fund (GRF12202124), Guangdong and Hong Kong Universities “1+1+1” Joint Research Collaboration Scheme (2025A0505000011), and Hong Kong Baptist University (HKBU) Seed Funding for Collaborative Research Grants (RC-SFCRG/23–24/SCI/06). Data were provided by the Human Connectome Project, WU-Minn Consortium (Principal Investigators: David Van Essen and Kamil Ugurbil; 1U54MH091657) funded by the 16 NIH Institutes and Centers that support the NIH Blueprint for Neuroscience Research; and by the McDonnell Center for Systems Neuroscience at Washington University. The authors thank Zhao Chang of the HKBU Department of Physics for his help in using the AMICO pipeline and interpreting its output.
Copyright © 2026 The Authors. Published by Elsevier Inc. All rights reserved.
PY - 2026/5/15
Y1 - 2026/5/15
N2 - Multiple lines of research have studied how complex brain dynamics emerge from underlying connectivity by using Ising models as simplified neural mass models. However, limitations on parameter estimation have prevented their use with individual, high-resolution human neuroimaging data. Furthermore, most studies focus only on connectivity, ignoring node heterogeneity, even though real brain regions have different structural and dynamical properties. Here we present an improved approach to fitting Ising models to 360-region functional MRI data: derivation of an initial guess model from group data, optimization of simulation temperature, and two stages of Boltzmann learning, first with group data, then with individual data. Our implementation uses GPU acceleration to mitigate the high computational cost of this approach. We then analyze how data binarization threshold affects goodness-of-fit, the role of the external field in model behavior, consistency among models fitted to different scans of the same individual, and correlations between model parameters and features from structural MRI, including measures of myelination and cortical folding. We find that binarizing fMRI data at higher thresholds decreases correlation between model and data functional connectivity but increases the heterogeneity of node external fields and their correlations with structural features. A choice of threshold that achieves both goodness-of-fit and intrinsic heterogeneity of regions results in a model that better reflects the reality of the brain as a network of intrinsically heterogeneous nodes. By enabling personalized, biophysically interpretable modeling of structure-function mapping across the whole brain, this approach can aid understanding of individual differences in brain network organization and bridge the gap between the network-focused methodology of connectomics and the region-focused paradigm typical of translational research.
AB - Multiple lines of research have studied how complex brain dynamics emerge from underlying connectivity by using Ising models as simplified neural mass models. However, limitations on parameter estimation have prevented their use with individual, high-resolution human neuroimaging data. Furthermore, most studies focus only on connectivity, ignoring node heterogeneity, even though real brain regions have different structural and dynamical properties. Here we present an improved approach to fitting Ising models to 360-region functional MRI data: derivation of an initial guess model from group data, optimization of simulation temperature, and two stages of Boltzmann learning, first with group data, then with individual data. Our implementation uses GPU acceleration to mitigate the high computational cost of this approach. We then analyze how data binarization threshold affects goodness-of-fit, the role of the external field in model behavior, consistency among models fitted to different scans of the same individual, and correlations between model parameters and features from structural MRI, including measures of myelination and cortical folding. We find that binarizing fMRI data at higher thresholds decreases correlation between model and data functional connectivity but increases the heterogeneity of node external fields and their correlations with structural features. A choice of threshold that achieves both goodness-of-fit and intrinsic heterogeneity of regions results in a model that better reflects the reality of the brain as a network of intrinsically heterogeneous nodes. By enabling personalized, biophysically interpretable modeling of structure-function mapping across the whole brain, this approach can aid understanding of individual differences in brain network organization and bridge the gap between the network-focused methodology of connectomics and the region-focused paradigm typical of translational research.
KW - Diffusion tensor imaging
KW - Functional MRI
KW - Ising model
KW - Structural MRI
KW - Structure-function mapping
KW - Whole-brain model
UR - https://www.scopus.com/pages/publications/105035454768
U2 - 10.1016/j.neuroimage.2026.121886
DO - 10.1016/j.neuroimage.2026.121886
M3 - Journal article
C2 - 41895545
AN - SCOPUS:105035454768
SN - 1053-8119
VL - 332
JO - NeuroImage
JF - NeuroImage
M1 - 121886
ER -