Computational disease models provide quantitative frameworks for understanding the complex dynamics of neurodegeneration in Alzheimer's disease (AD) and Parkinson's disease (PD). These models integrate molecular, cellular, and network-level mechanisms to simulate disease progression, test therapeutic interventions, and generate testable predictions. This mechanism page covers three major computational approaches: (1) neural network degeneration simulators, (2) protein aggregation kinetics models, and (3) network propagation models[@chen2021][@jakel2022].
Neural network degeneration simulators model how pathological proteins and cellular dysfunction spread through connected brain networks, disrupting neural activity and leading to cognitive and motor decline. These models build on the hypothesis that pathological proteins propagate along synaptic connections, causing network-level dysfunction that manifests as clinical symptoms[@stam_2020].
Computational disease models provide quantitative frameworks for understanding the complex dynamics of neurodegeneration in Alzheimer's disease (AD) and Parkinson's disease (PD). These models integrate molecular, cellular, and network-level mechanisms to simulate disease progression, test therapeutic interventions, and generate testable predictions. This mechanism page covers three major computational approaches: (1) neural network degeneration simulators, (2) protein aggregation kinetics models, and (3) network propagation models[@chen2021][@jakel2022].
Neural network degeneration simulators model how pathological proteins and cellular dysfunction spread through connected brain networks, disrupting neural activity and leading to cognitive and motor decline. These models build on the hypothesis that pathological proteins propagate along synaptic connections, causing network-level dysfunction that manifests as clinical symptoms[@stam_2020].
| Component | Description | Modeling Approach |
|-----------|-------------|------------------|
| Network topology | Structural and functional brain connectivity | Diffusion MRI, resting-state fMRI |
| Node dynamics | Neuronal activity in each brain region | Neural mass models, mean-field approximations |
| Pathology spread | Inter-node transmission of pathological proteins | Linear diffusion, prion-like propagation |
| Network failure | Breakdown of network integration/segregation | Graph theory metrics |
Network dynamics are typically modeled using coupled differential equations:
dX_i/dt = -X_i + f(Σ_j C_ij * X_j - μ_i) + P_i(t)
Where:
AD network models focus on:
PD network models incorporate:
Network models are validated against:
Protein aggregation kinetics models describe how misfolded proteins (amyloid-beta, tau, alpha-synuclein) nucleate, grow, and spread in the brain. These models use classical aggregation theory, nucleation-growth models, and prion-like propagation frameworks[@reach_2018].
Protein aggregation follows a nucleation-dependent mechanism:
| Stage | Description | Mathematical Model |
|-------|-------------|-------------------|
| Nucleation | Formation of stable oligomeric seeds | Homogeneous/heterogeneous nucleation rate k_n |
| Elongation | Addition of monomers to existing seeds | Elongation rate k_e |
| Secondary nucleation | New seeds from existing fibril surfaces | Rate k_2 |
| Fragmentation | Breakage creating new fibril ends | Rate k_f |
| Clearance | Autophagy, protease degradation | Rate k_c |
The master equation approach:
dM/dt = -k_n M^n - k_e M·N + k_f F - k_c M
dO/dt = k_n M^n - k_e M·O - k_c O
dF/dt = k_e M·N + k_e M·O - k_f F - k_c F
Where M = monomer, O = oligomer, F = fibril concentrations.
Aβ42 exhibits faster aggregation kinetics than Aβ40 due to:
Tau exhibits prion-like propagation:
α-Synuclein aggregation is central to PD pathogenesis:
| Feature | WT α-Syn | A53T α-Syn |
|---------|----------|------------|
| Nucleation rate | Baseline | ~5x faster |
| Oligomer toxicity | Moderate | High |
| Fibril formation | Slower | Faster |
| PD risk | Normal | Enhanced |
Key model parameters:
The reaction-diffusion model incorporates spatial spread:
∂C/∂t = D∇²C + R(C) - λC
Where:
Network propagation models combine connectivity data with pathology spread dynamics to predict spatiotemporal patterns of neurodegeneration. These models have been particularly successful in explaining Braak staging in PD and the spreading patterns of tau in AD[@vogel2023][@zhou2023].
The network diffusion equation models pathology spread across brain regions:
τ(t+1) = τ(t) + αC(τ(t) - τ(t-1)) - βτ(t)
Where:
Propagation follows connectivity patterns:
The Braak hypothesis proposes staging based on brainstem-to-cortex progression:
| Stage | Affected Regions | Years |
|-------|----------------|-------|
| 1-2 | Brainstem (dorsal motor nucleus, coeruleus) | 0-2 |
| 3-4 | Substantia nigra, basal forebrain | 2-5 |
| 5-6 | Cortex (motor, premotor, association) | 5-12 |
Network propagation models reproduce this pattern when:
Tau propagation models incorporate:
| Parameter | AD (Tau) | PD (α-Syn) |
|-----------|-----------|------------|
| Propagation rate (α) | 0.08-0.15 | 0.10-0.20 |
| Clearance rate (β) | 0.03-0.08 | 0.02-0.05 |
| Origin tau burden | 0.3-0.5 SUVr | Baseline |
| Time to cortex | 8-15 years | 10-20 years |
Network propagation models are validated against:
Deterministic models fail to capture individual variability. Stochastic models incorporate:
dτ_i = αΣ_j C_ij τ_j dt - βτ_i dt + σdW_i
Where σdW_i = Wiener process for noise.
Individual predictions require Monte Carlo approaches:
Comprehensive computational models integrate across scales:
| Scale | Modeling Method | Time Scale |
|-------|--------------|----------|
| Molecular | MD simulations, Markov models | ns - ms |
| Cellular | ODE/PDE models | hours - days |
| Network | Neural mass models | ms - seconds |
| System | Clinical progression models | years - decades |
Computational models enable:
Individualized models incorporate:
| Intervention | Model Prediction |
|-------------|--------------|
| Anti-Aβ antibodies | Reduce Aβ burden 30-50% |
| Anti-tau immunotherapy | Slow propagation 40-60% |
| α-Synuclein aggregation inhibitor | Reduce oligomer toxicity |
| Deep brain stimulation | Normalize network dynamics |
Computational models inform:
Modern approaches incorporate:
Individual digital twins integrate:
Integration with: