The prion hypothesis states that misfolded proteins can act as infectious agents that template the misfolding and aggregation of healthy proteins to transmit a disease. Increasing evidence suggests that pathological proteins in neurodegenerative diseases adopt prion-like mechanisms and spread across the brain along anatomically connected networks. Local kinetic models of protein misfolding and global network models of protein spreading provide valuable insight into several aspects of prion-like diseases. Yet, to date, these models have not been combined to simulate how pathological proteins multiply and spread across the human brain. Here, we create an efficient and robust tool to simulate the spreading of misfolded protein using three classes of kinetic models, the Fisher–Kolmogorov model, the Heterodimer model and the Smoluchowski model. We discretize their governing equations using a human brain network model, which we represent as a weighted Laplacian graph generated from 418 brains from the Human Connectome Project. Its nodes represent the anatomic regions of interest and its edges are weighted by the mean fibre number divided by the mean fibre length between any two regions. We demonstrate that our brain network model can predict the histopathological patterns of Alzheimer’s disease and capture the key characteristic features of finite-element brain models at a fraction of their computational cost: simulating the spatio-temporal evolution of aggregate size distributions across the human brain throughout a period of 40 years takes less than 7 s on a standard laptop computer. Our model has the potential to predict biomarker curves, aggregate size distributions, infection times, and the effects of therapeutic strategies including reduced production and increased clearance of misfolded protein.

No abstract

The prion hypothesis states that misfolded proteins can act as infectious agents that trigger the misfolding and aggregation of healthy proteins to transmit a variety of neurodegenerative diseases. Increasing evidence suggests that pathogenic proteins in Alzheimer's disease adapt prion-like mechanisms and spread across the brain along an anatomically connected network. Local kinetics models of protein misfolding and global network models of protein diffusion provide valuable insight into the dynamics of prion-like diseases. Yet, to date, these models have not been combined to simulate how pathological proteins multiply and spread across the human brain. Here we model the prion-like spreading of Alzheimer's disease by combining misfolding kinetics and network diffusion through a connectivity-weighted Laplacian graph created from 418 brains of the Human Connectome Project. The nodes of the graph represent anatomic regions of interest and the edges represent their connectivity, weighted by the mean fiber number divided by the mean fiber length. We show that our brain network model correctly predicts the neuropathological pattern of Alzheimer's disease and captures the key characteristic features of whole brain models at a fraction of their computational cost. To illustrate the potential of brain network modeling in neurodegeneration, we simulate biomarker curves, infection times, and two promising therapeutic strategies to delay the onset of neurodegeneration: reduced production and increased clearance of misfolded protein.

The prion-like hypothesis of neurodegenerative diseases states that the accumulation of misfolded proteins in the form of aggregates is responsible for tissue death and its associated neurodegenerative pathology and cognitive decline. Some disease-specific misfolded proteins can interact with healthy proteins to form long chains that are transported through the brain along axonal pathways. Since aggregates of different sizes have different transport properties and toxicity, it is important to follow independently their evolution in space and time. Here, we model the spreading and propagation of aggregates of misfolded proteins in the brain using the general Smoluchowski theory of nucleation, aggregation, and fragmentation. The transport processes considered here are either anisotropic diffusion along axonal bundles or discrete Laplacian transport along a network. In particular, we model the spreading and aggregation of both amyloid-β and τ molecules in the brain connectome. We show that these two models lead to different size distributions and different propagation along the network. A detailed analysis of these two models reveals the existence of four different stages with different dynamics and invasive properties.