Heterogeneity in human brain clearance adds resilience against tauopathy -- a computational model informed by glymphatic MRI

Neurotoxic protein fragments such as amyloid-beta and tau accumulate in characteristic staging patterns in Alzheimers disease (AD). The brain clears such metabolic substances via multiple different systems, including via the glymphatic (extracellular/extravascular) pathway. Here, we ask how the distinct features that characterize human glymphatic function would affect the prion-like cascade of protein invasion associated with AD. To address this question, we extract and analyze individual clearance rates from human glymphatic MRI (gMRI) data sets. These clearance rates define subject-specific maps of glymphatic clearance that vary both across cortical lobes and Braak staging regions. We apply these clearance maps as initial states in a computational network model linking misfolded proteins, tissue damage, and local clearance to simulate a series of individual proteinopathy trajectories. Our results show that the spatial heterogeneity in initial clearance induces characteristic propagation patterns, delaying and redirecting the disease progression. Moreover, reducing this spatial heterogeneity accelerates disease progression and induces staging patterns typically associated with AD. A comparison between well-rested subjects and subjects who underwent a single night of sleep deprivation did not reveal differences in initial clearance maps nor in simulated disease progression. These findings suggest that spatial heterogeneity in brain clearance may be a key factor for neurodegenerative resilience.