Computational Causal Modeling of the Dynamic Biomarker Cascade in Alzheimer's Disease

Background. Alzheimer's disease (AD) is a major public health concern, and there is an urgent need to better understand its complex biology and develop effective therapies. AD progression can be tracked in patients through validated imaging and spinal fluid biomarkers of pathology and neuronal loss. We still, however, lack a coherent quantitative model that explains how these biomarkers interact and evolve over time. Such a model could potentially help identify the major drivers of disease in individual patients and simulate response to therapy prior to entry in clinical trials. A current theory of AD biomarker progression, known as the dynamic biomarker cascade model, hypothesizes AD biomarkers evolve in a sequential but temporally overlapping manner. A computational model incorporating assumptions about the underlying biology of this theory and its variations would be useful to test and refine its accuracy with longitudinal biomarker data from clinical trials. Methods. We implemented a causal model to simulate time-dependent biomarker data under the descriptive assumptions of the dynamic biomarker cascade theory. We modeled pathologic biomarkers (beta-amyloid and tau), neuronal loss biomarkers, and cognitive impairment as nonlinear first-order ordinary differential equations (ODEs) to include amyloid-dependent and nondependent neurodegenerative cascades. We tested the feasibility of the model by adjusting its parameters to simulate three specific natural history scenarios in early-onset autosomal dominant AD and late-onset AD and determine whether computed biomarker trajectories agreed with current assumptions of AD biomarker progression. We also simulated the effects of antiamyloid therapy in late-onset AD. Results. The computational model of early-onset AD demonstrated the initial appearance of amyloid, followed by biomarkers of tau and neurodegeneration and the onset of cognitive decline based on cognitive reserve, as predicted by the prior literature. Similarly, the late-onset AD computational models demonstrated the first appearance of amyloid or nonamyloid-related tauopathy, depending on the magnitude of comorbid pathology, and also closely matched the biomarker cascades predicted by the prior literature. Forward simulation of antiamyloid therapy in symptomatic late-onset AD failed to demonstrate any slowing in progression of cognitive decline, consistent with prior failed clinical trials in symptomatic patients. Conclusions. We have developed and computationally implemented a mathematical causal model of the dynamic biomarker cascade theory in AD. We demonstrate the feasibility of this model by simulating biomarker evolution and cognitive decline in early- and late-onset natural history scenarios, as well as in a treatment scenario targeted at core AD pathology. Models resulting from this causal approach can be further developed and refined using patient data from longitudinal biomarker studies and may in the future play a key role in personalizing approaches to treatment.