A continuous-time mathematical model and discrete approximations for the aggregation of β-Amyloid

Alzheimer's disease is a degenerative disorder characterized by the loss of synapses and neurons from the brain, as well as the accumulation of amyloid-based neuritic plaques. While it remains a matter of contention whether β-amyloid causes the neurodegeneration, β-amyloid aggregation is associated with the disease progression. Therefore, gaining a clearer understanding of this aggregation may help to better understand the disease. We develop a continuous-time model for β-amyloid aggregation using concepts from chemical kinetics and population dynamics. We show the model conserves mass and establish conditions for the existence and stability of equilibria. We also develop two discrete-time approximations to the model that are dynamically consistent. We show numerically that the continuous-time model produces sigmoidal growth, while the discrete-time approximations may exhibit oscillatory dynamics. Finally, sensitivity analysis reveals that aggregate concentration is most sensitive to parameters involved in monomer production and nucleation, suggesting the need for good estimates of such parameters.