Studying the Progression of Amyloid Pathology and Its Therapy Using Translational Longitudinal Model of Accumulation and Distribution of Amyloid Beta.

Long-term effects of amyloid targeted therapy can be studied using a mechanistic translational model of amyloid beta (Aβ) distribution and aggregation calibrated on published data in mouse and human species. Alzheimer disease (AD) pathology is modeled utilizing age-dependent pathological evolution for rate constants and several variants of explicit functions for Aβ toxicity influencing cognitive outcomes (Adas-cog). Preventive Aβ targeted therapies were simulated to minimize the Aβ difference from healthy physiological levels. Therapeutic targeted simulations provided similar predictions for mouse and human studies. Our model predicts that: (1) at least 1 year (2 years for preclinical AD) of treatment is needed to observe cognitive effects; (2) under the hypothesis with functional importance of Aβ, a 15% decrease in Aβ (using an imaging biomarker) is related to 15-20% cognition improvement by immunotherapy. Despite negative outcomes in clinical trials, Aβ continues to remain a prospective target demanding careful assessment of mechanistic effect and duration of trial design.

WHAT IS THE CURRENT KNOWLEDGE ON THE TOPIC?

☑ Multiple clinical trials targeting amyloid, the hallmark of AD pathology, have failed. The key efficacy markers are CSF amyloid and brain PET, whereas neither clear correlations between them and cognitive scores have been proven nor quantification of amyloid toxic effect proposed.

WHAT QUESTION DOES THIS STUDY ADDRESS?

☑ The questions this study address are (1) whether a mechanistic translational model can allow for prediction of long‐term clinical trials and (2) what efficacy on cognition would be demonstrated under different amyloid toxicity hypotheses.

WHAT THIS STUDY ADDS TO OUR KNOWLEDGE

☑ The mechanistic model allows for prediction as to whether long‐term therapy leads to normalization of amyloid levels in different compartments. Quantitative formulation of toxicity hypotheses allows for using the trial results for hypothesis testing.

HOW MIGHT THIS CHANGE DRUG DISCOVERY, DEVELOPMENT, AND/OR THERAPEUTICS?

☑ Our model, validated on multiple pathology progression data, allows for more comprehensive design of clinical trials. Variants of scenarios can be analyzed depending on the disease stage and type of therapy and optimized for higher efficacy and safety.

Alzheimer disease (AD) leading to the decrement of cognitive functions is hypothesized to be caused by abnormal accumulation of amyloid beta (Aβ) peptide,1 which is found in soluble and insoluble forms in the postmortem brains of subjects with AD. Two peptide variants, Aβ40 and Aβ42, produced from amyloid precursor protein (APP), form oligomers and polymer chains,2 which differ by physicochemical properties, including solubility. The brain levels of soluble and insoluble forms of Aβ are distinguished in subjects with AD with respect to healthy individuals.3 Main techniques to observe the progression of amyloid pathology are the measurement of Aβ changes in cerebrospinal fluid (CSF), which is used as a reference for changes in brain interstitial fluid (BIF), and positron emission tomography (PET) with Aβ‐specific tracers (e.g., PiB), quantifying the insoluble amyloid in brain. PET data demonstrate an increase of insoluble amyloid, whereas CSF Aβ42 concentration is decreased,4 hypothetically due to the absorption of free Aβ into insoluble plaques. As AD progresses, the insoluble Aβ forms seem to reach a steady‐state value.5 Soluble and insoluble amyloid neurotoxicity6, 7, 8 has been explored in vitro.9 Nucleation‐dependent polymerization was shown to be required for neurotoxicity,10 and clinical studies reveal correlation between amyloid imaging load registered by PET and cognitive function (Adas‐cog score11). Although controversial results have been reported,12 postmortem soluble amyloid correlated highly with disease severity,13 and elevated oligomers were found in the CSF of subjects with AD.14 Conflicting hypotheses have been proposed that Aβ is not the main reason for pathology but is just the trace of other neurodegenerative events.15 Moreover, soluble Aβ may regulate memory functions by modulating synaptic plasticity at picomolar range from 20 to about 300 pM.16 Nevertheless, researchers concur that Aβ seems to be an important link in progression of AD pathology, an appropriate diagnostic biomarker and pharmacological target for treatment.17

Treatment options targeting amyloid clearance or production are under study.18 These therapies lead to the attenuation of cognitive decline in nonclinical studies,19 but have not shown significant benefit in clinical trials and have caused intolerable side effects in some cases.18 The relationship between soluble and insoluble Aβ depends on formation of polymerization seeds, rates of aggregation, and fibril breakage,20, 21, 22 and the rate of destruction of mature fibrils.22, 23 The aggregation of amyloid depends nonlinearly on the concentration of the soluble or free Aβ, which, in turn, depends on the production, degradation, and distribution between different compartments.

The detailed Aβ distribution kinetics and drug effect on the brain soluble and insoluble Aβ are explored only in nonclinical studies. Translation from the nonclinical species to humans requires estimation of the level of the similarity and difference between nonclinical AD models (primarily transgenic mouse models) and AD pathophysiology on a mechanistic level.

Prediction of specific mechanistic effects (through therapeutic targets) and clinical trial design optimization may be facilitated by a mechanistic model taking into account the aforementioned features. Human data should be used for verification, which was not done in existing models of amyloid accumulation.21, 24, 25

In this analysis, a quantitative description of Aβ aggregation was developed and integrated with pre‐existing distribution model26 to describe longitudinal changes of soluble and insoluble Aβ both for healthy subjects and subjects with AD. The model would allow for exploration into how the levels of brain Aβ are modified due to the changes in Aβ production, clearance, or other related processes during pathology or therapeutic interventions. It can be the basis for exploration of biomarker behavior and for prediction into clinical trials. Other hallmarks of AD pathology exist (e.g., neurofibrillary tau tangles, oxidative or inflammatory markers), many of which interact with Aβ. Two major variants of the empirical connection between cognition and amyloid are considered: (1) neurotoxicity (loss of function and apoptosis) is directly dependent on amyloid; and (2) pathological neurodegenerative process develops independently after Aβ achieves some threshold. Additionally, this model provides the ability to analyze key differences between the Tg2576 mouse AD model and human biology to make conclusions about the translation from preclinical studies into human disease understanding.

METHODS

Model structure

Two forms of Aβ (only key species participating in aggregation, Aβ40 and Aβ42) are considered in the model. The model of Aβ distribution is described in detail in ref. 26. It was extended by the addition of oligomerization and aggregation processes, as well as degradation of insoluble forms, which are presented in brain cells (BCs) and BIF (Figure 1). Differential equations for the aggregated amyloid have been derived in analogy with the standard approach for the polymerization description.20, 21, 22 The full explanation of the model is given in the Supplementary Figure S1. Model development and calibration stages are described in Table 1.

Figure 1

Sketch scheme of the model. (a) Description of the elementary processes of the polymerization and oligomerization model before reduction. Designations: D, dimers; M, monomers of amyloid beta (Aβ); N, nuclei; Pi, fibril of i‐th length; Tr, trimmers; Tt, tetramers. (b) Sketch of the processes in different compartments described in the model. The red arrows show the simplified scheme of the reduced aggregation model (see part A of the Supplementary for details). Compartments: BC, brain cells; BIF, brain interstitial fluid; CSF, cerebrospinal fluid; OT, other tissues; PL, – plasma.

Table 1

Stages of model fitting and validation using different datasets

#	Description	Parameters fitted	Data used	No. of data points
				Fitting	EV
F0	Parameters taken from distribution model	Distribution parameters	See ref. 26	‐//‐	‐//‐
F1	Fitting of parameters for subjects with AD using integrated model	Parameters of age dependence function (manually) (3) Aβ release, aggregation, and longitudinal change (20)	1) Soluble and insoluble Aβ in BR, Aβ in CSF, and PL for AD >60 years and healthy <60 years	29 (+4 +100b)a	0
		Tracer parameters taken from literature	2) PiB‐PET data for AD	0	83
F2	Fitting of parameters for HS, choice of the best model	Longitudinal parameters (12 different variants from 1–5 parameters)	1) Soluble and insoluble Aβ in BR for HS	10 +250b	0
			2) CSF for HS (OFV, AIC)	0	171b
			3) PiB‐PET for HS	0	527
F3	Translation of the model to Tg2576 mouse	Parameters of Age‐dep function (manually) (3)	Qualitative view of the temporal characteristics of the dataset	‐	0
		Release scaling factors (4)	Soluble and insoluble in BR, Aβ in CSF and PL before 6 months	10	0
		Longitudinal parameters (2)	Soluble and insoluble in BR, Aβ in CSF, and PL	70	0
F4	Translation of the model to WT mouse	Release scaling factors (3)	Steady‐state data for soluble and insoluble Aβ in BR, Aβ in CSF, and PL	4	0

Aβ, amyloid beta; AD, Alzheimer disease; Age‐dep, age‐dependent; AIC, Akaike information criterion; BR, brain; CSF, cerebrospinal fluid; EV, external verification (comparison of data with predictions); HS, healthy subjects; OFV, objective function value; PET, positron emission tomography; PiB, Pittsburgh compound B; PL, plasma; WT, wild‐type.

Data from healthy population (before 60 years of age) used during fitting of AD; bData for fitting (validation) transformed to medians.

Model for longitudinal changes

Different hypotheses for Aβ accumulation in sporadic AD, depending on age, have been accounted for in the model: the decrease of Aβ clearance,27, 28, 29 shift of Aβ40 and Aβ42 production in the brain,30, 31 raise of nucleation,32 and the downfall of insoluble Aβ destruction.33, 34 They are described empirically as explicit time dependence of corresponding rate constants Pi with baseline values until 60 years: where with specific coefficients (amplitudes, ) for each process hypothesized to contribute to the pathology (parameter of the model). Here, Ktime, prog, and tdecline are parameters of the curve, and t is the time in years.

Description of the Aβ influence on cognition

The Aβ toxicity was described under simplifying assumptions: (1) the brain soluble Aβ42 (monomer and soluble oligomers) is the key toxic form; (2) simple linear models were chosen as the cognition dataset is sparse, thereby limiting additional assumptions on model structure; (3) there are no irreversible changes, thus, a decrease in Aβ42 would lead to an immediate recovery of cognitive function; and (4) parameters have been chosen empirically to achieve gradual change of the Adas‐cog score from 8 to ∼3511 during the overall disease progression.

A potential physiological function of Aβ42 at low concentrations is also considered in the model.16

The cognitive status depends on the quantity of the nondegenerated neurons (N) and on their functional capacity (fsol): where N decreases with Aβ42 accumulation and may decrease with Aβ42 depletion. In this case, change of cognition can be described as:

Considered toxicity hypotheses expressed as are shown in Table 2 and described in detail in the Supplementary Information, as well as derivation of function .

Table 2

Hypotheses of the amyloid toxicity.

Aβ toxicity hypothesis	Aβ is requireda	#	Comments	P value for rejection of hypothesis based on datasets
				1	2	3	4	5
Proportional	No	1.1	toxAb_42=k_dr·(AsolubleBR42−base_health), base_health=1 nM	< 0.001	< 0.001	< 0.001	< 0.001	1
	Yes	2.1	Same as 1.1, but takes into account decrease of cognition for the low Ab42 levels	1	1	1	1	1
Threshold	No	1.2	toxAb_42=k_thr∗t−t_st if AsolubleBR42>base_health	< 0.001	< 0.001	< 0.001	< 0.001	1
	Yes	2.2	Same as 1.2, but takes into account decrease of cognition for the low Ab42 levels	1	1	1	1	1
Proportional and threshold	No	1.3	Combination of 1.1 and 1.2	0.2561	< 0.001	0.1116	< 0.001	1

Simulations using hypotheses have been considered with data (simulations based on data mean and SD) using Wilcoxon one‐side criterion for unpaired observations. The null hypothesis was that model overestimates treatment effect (Adas‐cog difference) not higher than by 1 point. Datasets: Adas‐cog from avagacestat data set (Coric et al.,39 2012):1–50 mg, 12 weeks, 2–125 mg, 12 weeks, 3–50 mg, 24 weeks, 4–125 mg, 24 weeks; bapineuzumab data (Rinne et al.,40 2010), Adas‐cog for 78 weeks.

Aβ, amyloid beta.

Hypothesis accounts for positive influence of Aβ on neuronal functioning.

Mechanism‐specific simulations

Simulations for specific mechanistic targets, amyloid release inhibition (RI) and destruction of insoluble amyloid (DT), were performed in different regimens: corresponding parameter of a target process (release rate constant or destruction rate constant) was multiplied on the function of the drug effect EFF: for the constant therapeutic intervention, termed a “cutting regimen” (CR), or: for adaptive, or a “progressive regimen” (PR), with therapeutic intervention increasing gradually to compensate for the pathological change of corresponding parameter (see paragraph two of the Methods section). Therapy influenced equally Aβ42 and Aβ40 forms. To test the toxicity hypotheses, simulated mechanistic intervention (simulated therapy) was switched on at 70 (preclinical) and 75 years (mild AD), which are the approximate ages of cognitive symptom manifestations for sporadic cases.

Statistical analysis

A one‐sided nonparametric (Wilcoxon) test for unpaired observations was used to compare simulations based on hypotheses with clinical trials. Observations were simulated using mean and SD from publications. To note, due to assumption (3) of the Aβ influence on cognition, the model predicts the most optimistic scenario, tending to overestimate the treatment effect. Software R, version 3.2, was used for all statistical procedures.

RESULTS

Description of patients with AD data

The model sufficiently characterized Aβ42 and Aβ40 data in patients with AD (Figure 2 a, Supplementary Figures S5 and S6), the medians of the data are within the model prediction bands for the observed medians. The CSF Aβ42 before 60 years is not well captured by the model, because age <60 years is not attributed to sporadic AD, as considered in this model. The model predicts a decrease of CSF Aβ42 vs. healthy baseline after 60 years (Figure 2 a, upper panel), which is not observed for Aβ40 (Supplementary Figure S5). Destruction of insoluble amyloid (Supplementary Figure S2) and Aβ42 production rate have the greatest amplitude of longitudinal change among considered processes (Supplementary Table S3). The decrease of the insoluble amyloid degradation rate constant during AD in the model is consistent with microglial phagocytosis descent.23 Plasma Aβ increase (Supplementary Figure S7) is mainly determined by the decline in the Aβ clearance with the minimal contribution from the brain processes.26 The model was verified using independent PET DVR dataset for subjects with AD (Figure 2 b, left panel).

Figure 2

Model calibration and performance. (a) Results of the fitting of the Alzheimer disease (AD; left) individuals and healthy (right) individuals with 95% confidential bands for uncertainty of model parameters. For the healthy individuals, results from the model with a single parameter recalibration (see Supplementary Table S2, H14) are shown. The blue and red line lines (upper right) correspond to least square analysis from published sources. Median amyloid beta (Aβ) in cerebrospinal fluid (CSF) for AD (upper left) was taken only for points after 60 years. (b) Verification of the model for AD (left) and healthy (right) individuals by the positron emission tomography (PET) DVR data from different sources with confidence bands for model predictions determined by uncertainty of parameters. For verification, we chose data from cortex and hippocampus with cerebellum used as reference region for PiB binding estimation. Colors denote different data sources (see complete reference list in Supplementary Figures S6–S9). All y‐scales are log10‐scales. BR, brain.

Description of the healthy subject data

Different variants of parameter sets were explored to describe the healthy subject data (Supplementary Table S2). This exploration revealed key factors distinguishing healthy subjects and subjects with AD. Recalibration of the progression rate for healthy subjects (Supplementary Table S2, H14) allowed for the adequate description of the data used for fitting and verification (CSF data on Figure 2 a). Although there is significant variability across healthy subject DVR data, the 95% confidence band for model prediction contains the median of the healthy subject data (Figure 2 b, right panel). The BIF soluble Aβ42 goes down (Supplementary Figure S10) in subjects with AD, coinciding with decrease in CSF. The BIF Aβ42 concentration in healthy subjects is not lower than 0.01 nM and BC soluble Aβ42 concentration does not exceed 1 nM across multiple variants of the model parameterization. This is in contrast to subjects with AD, whereas BIF Aβ42 decreases, whereas BC Aβ increases. The BIF amyloid concentration of 0.01 nM is close to the lower limit of positive LTP regulation from in vitro experiments.16 Insoluble Aβ42 and Aβ40 concentrations change similarly in healthy and AD populations, but, unlike in the AD population, the simulated confidence band for the healthy individuals does not exceed concentration of 300 nM, corresponding to the threshold value of amyloid burden detection.35 This model variant under predicts part of the data for soluble and insoluble amyloid, but adequately describes the independent dataset on CSF Aβ42 and Aβ40 concentrations used for validation.

Other models (Supplementary Table S2, Supplementary Figure S10) predict multiple variants of amyloid behavior in the BIF and BC of healthy subjects. Some models describe higher levels of soluble Aβ42 (Supplementary Figure S10) and demonstrate better values of the fitting criterions (Supplementary Table S2). Model H14 was chosen due to its simplicity and the very good validation ratio. It may represent a healthy population, which also moves toward the AD state but at a very low rate. Plasma data were used neither for fitting nor for external verification of the model for healthy subjects as there is no difference between healthy and AD data,36 and plasma has minimal influence on the brain amyloid.26 The amplitude of the amyloid clearance decrease was assumed to be the same for both groups. Based on the model output, the decrease of plasma Aβ by 80% does not change the concentration of the brain Aβ forms (Supplementary Figure S13).

Translation of the model from human to the mouse Tg2576 and wild type mouse

The final human model parameters were used to describe amyloid accumulation over the rodent lifespan. The amyloid synthesis parameters and longitudinal effect parameters were recalibrated to allow for the satisfactory description of the rodent data (Supplementary Figures S11 and S12).

Simulations of therapy

The decrease in insoluble Aβ42 accumulation (Figure 3 a, bottom panel) during the RI therapy is preceded by an abrupt decrease of soluble Aβ42 (Figure 3) below physiological range. Transgenic mouse simulation results are similar, but BC Aβ42 concentration is still higher than for the wild type mouse, and insoluble Aβ42 does not decrease from baseline. Destruction therapy leads to the increase of the soluble forms of Aβ42 in BIF both in human and mouse models (Figure 3 , middle row). Time when an intervention is initiated (prevention at 60 years or treatment at 75 years) does not influence the new amyloid steady‐state level or time of its achievement: therapy with early start and the late therapy lead to the same state at ∼75–80 years. Steady state is determined by the intrinsic time of disease progression and intervention type. Although BIF soluble Aβ42 concentration increases during DT, BC Aβ42 level decrease when compared to the AD placebo (for preclinical AD such effect appears after 10 years). To note, this effect was not predicted in the mouse model.

Figure 3

Simulation of different types of treatment (inhibition of production by 80% maximally and activation of destruction of plaques to 600%) for subjects with Alzheimer disease (a) and Tg2576 mouse (b) according to “cutting regime” (CR) and “progressive regime” (PR). For PR, the maximal value of therapeutic intervention is the same as for the CR, age‐dependence parameters are the same as for the disease progression description (see Methods and Supplementary Table S3). Therapy simulations are compared with levels of normal values (light green area on a) and the model forecasted trajectory for normal subjects (dashed green line). Red arrows indicate the time at which the biomarker was no longer contained in the region of normal values (a). Two ages of therapy start were used for humans: 60 years (left) and 75 years (right). Black and cyan lines are overlapping on the right column panels, predicting absence of difference for two regimens of DT at this stage, the same for release inhibition (RI; orange and violet). (b) Two ages of therapy start were used for the mouse models: 6 months (left) and 15 months (right). Comparison with the wild‐type (WT) mouse (green dashed line) is given (b). Aβ, amyloid beta; BC, brain cell; BIF, brain interstitial fluid.

Attenuating soluble amyloid in the brain 10 times below normal levels resulted in insoluble Aβ42 levels higher than normal after 10 years (Figure 3 a, bottom panel). A large inhibition of Aβ production (CR) before disease progression (prevention targeted treatment) leads to almost zero difference in BIF Aβ42 concentrations vs. AD placebo after 10 years of therapy (Figure 3 a, middle, left panels), although BC Aβ42 and insoluble amyloid differ significantly from placebo (Figure 3 a, top and bottom, left panels).

The magnitude of “progressive” therapy depends on the degeneration stage of the system and, thus, makes effect on the system less radical with fewer declines from the healthy physiological state. DT (PR) with titration from zero at 60 years to 600% at 80 years (Figure 3 , left panel) yields much lower growth of soluble Aβ42 in BIF than CR (“cutting” therapy) causes. Such preventive therapy makes sense only at the early stage of the disease, when levels of the most of the biomarkers are indistinguishable from control.

Testing the hypotheses of Aβ toxicity

The Adas‐cog model performed satisfactorily on collected data (Supplementary D and Figures S14 and S15). All hypotheses allow for approximate description of AD progression (Supplementary Figure S15). Models predict that starting from the value about 20 Adas‐cog will change by ∼5 point score during the next 2 years, which is similar to the estimates reported in the literature.37, 38 Each targeted mechanism (RI or DT) shows differing trajectories depending on the disease severity (Figure 4). For the preclinical AD simulations, minimal effect, as compared to placebo, would be expected for ∼2 years, but a change in the biomarker could be observed earlier. The simulations for the mild AD population suggest a more rapid therapeutic effect depending on the hypothesized toxicity. Here, Aβ production inhibition as a possible therapy would yield a more rapid effect if the toxicity follows the proportional or threshold model, but the effect is negated if a minimal concentration of Aβ42 is required for neuronal function. Threshold hypotheses predict no positive effect unless inhibition is high enough to drop Aβ42 below the toxic threshold (50% inhibition), then during the first month of treatment a large (5–10) point difference would be expected. If amyloid is functionally important (hypotheses 2.1 and 2.2 of Table 2), then an initial cognitive decline would be expected for more intense RI therapy (simulated at 50% inhibition). This apparent deterioration would be expected to decline over time leading to a net improvement over 5 years. No negative effect is predicted at low level of destruction activation for any of the hypotheses.

Figure 4

Effect of different therapies on cognitive function and other biomarkers under different hypotheses of amyloid toxicity (see Table 2). Two variants of therapy initiation were applied: preclinical stage (70 years, left) and mild Alzheimer disease (AD) stage (75 years, right).(a) Amyloid beta (Aβ) toxicity influencing cognitive outcomes (Adas‐cog) simulations under different hypotheses. (b) Influence of different treatments on the cognitive function, cerebrospinal fluid (CSF) amyloid and positron emission tomography (PET) DVR with respect to baseline under hypothesis 2.1 (Table 2). Ab, amyloid 42; RI, release inhibition.

Simulated CSF Aβ42 dynamics differ between AD populations (Figure 4 b). No decrease of PET DVR vs. baseline is predicted with any of the proposed therapies in the preclinical AD population, although the difference from placebo is expected to grow with time. For the mild AD simulation, a 200% increase of Aβ destruction is expected to cause a 15% DVR decrease as compared to baseline within 2 years (Figure 4 b, bottom right panel) as well as cognitive improvement by ∼20%, compared to baseline (as predicted by hypothesis 2.1). Simulated CSF Aβ42 concentration change vs. baseline is small (increase by 2–10%) for mild AD group and seems to be nonmonotonic for preclinical AD (Figure 4 b, middle left panel).

We simulated avagacestat clinical trial data,39 applying levels of 20% and 50% production inhibition, corresponding to CSF Aβ42 inhibition observed at doses 50 and 125 mg, respectively. Validity of each hypothesis was tested by the calculation of P value for the difference between data and model predictions. As our model predictions tend to be optimistic due to assumptions (see Methods), the null hypothesis for each test was taken as absence of significant overprediction (more than 1 point), instead of simple coincidence of predictions vs. observations. Hypotheses without Aβ42‐positive influence (Table 2 , 1 .1–3) can be rejected as they overestimate treatment effect (P value < 0.001). Hypotheses 2.1 and 2.2 cannot be rejected (P value = 1 in Table 2), and simulations under these hypotheses are more or less close to the results of the trial (Figure 5), predicting probability of cognitive decline for higher doses. The simulations for DT (increase to 150%) under toxicity hypotheses 2.1 and 2.2 are similar to the observed results of the bapineuzumab study40 (Supplementary Figure S14), which show about 15% DVR decrease at 78 weeks (Figure 4 b, bottom right panel). In contrast to avagacestat data, the bapineuzumab trial results do not allow for rejection of any hypotheses.

Figure 5

Comparison of results of amyloid beta (Aβ) toxicity influencing cognitive outcomes (Adas‐cog) simulation by hypotheses 2.1 and 2.2 and results of avagacestat clinical trial.39 Different doses of avagacestat were given to patients of about 75 years of age, and cerebrospinal fluid (CSF) changes vs. baseline at 24 weeks and Adas‐cog score changes at 12 and 24 weeks with respect to baseline were measured. Confidence intervals for simulation calculated through the uncertainty of the aggregation and longitudinal parameters. Adas‐cog change from baseline (points) and for CSF (Aβ42% change from baseline*10−1) are displayed. For observations, SD is shown.

DISCUSSION

The model presented here describes the accumulation of Aβ in the brain and longitudinal changes in other compartments for human, healthy individuals, patients with AD, and transgenic mice Tg2576. It was constructed mechanistically, calibrated, and verified on a large number of clinical and animal data.

Our model provides a unique possibility to explore the differences between healthy subjects and subjects with AD in terms of quantities that are difficult to be measured, Aβ concentrations in BIF, and BC, which may differentially contribute to toxic mechanisms, and CSF biomarkers. Concentrations of insoluble brain Aβ and soluble Aβ42 in BC exceed 1 µM and 1 nM, respectively, in the AD population, and remain below these levels in healthy subjects. Soluble Aβ42 in the BIF of subjects with AD falls to subphysiological levels due to aggregation (Figure 3 a, middle panel). It may raise the question whether the plaques themselves, or even oligomers, are the main contributors into the neuron pathology. For example, depletion of Aβ42 soluble pool due to its absorbance to plaques may lead to the loss of positive function observed for the picomolar Aβ42.41 In addition, activity‐dependent modulation of endogenous Aβ may normally participate in a negative feedback keeping neuronal hyperactivity in check.42 One of the therapeutic options may be to stabilize Aβ levels close to a physiologically relevant value. Plaque destruction increase (to 600%) would be able to normalize insoluble Aβ (to about 300 nM) in brain and soluble Aβ in BIF during 1 year. Potentially after such stabilization of plaques, additional inhibitory therapy could be applied to decrease brain Aβ. The CSF Aβ42 would be a sufficient biomarker for such a treatment and should neither exceed the control level nor significantly fall down, whereas PET‐PiB binding should decrease.

Model predicts that the main reason of Aβ accumulation is a failure of insoluble Aβ destruction. Although destruction processes, such as enzyme proteolysis, astrocyte, and microglial phagocytosis, are not described mechanistically in our model, many hypotheses exist that the breakdown of these processes may be the driving force for aging,43 and, thus, the model insight seems reasonable. To note, this model was constructed to describe the sporadic AD population, therefore, populations with APP mutations are to be described by the shifts in production, degradation, polymerization, and oligomerization parameters.

Transgenic mouse data were fitted by parameters of the release and Aβ40 polymerization. This presumably may reflect already observed distinct physical and chemical properties of peptides generated in the mouse AD model.44

Soluble Aβ42 in Tg2576 will not decrease below physiological levels corresponding to wild‐type mice until 2–3 months of the treatment. This is expected as the initial concentrations and production of amyloid in Tg2576 strongly exceeds those in the wild‐type mouse.45

Only several designs may be used for actual confirmation or rejection of Aβ toxicity hypotheses (Figure 4 a, Table 2). Nevertheless, outcomes of other designs may still inform about particular properties of the system. For example, positive clinical outcomes (Adas‐cog decrease) at light inhibition may suggest that Aβ inputs in degeneration as a primary driver of disease, whereas negative results at low RI would attest in favor of threshold (or prion‐type) hypothesis.

Initiation of early therapy (Figure 4 a, left panel) would not be expected to demonstrate any differences vs. placebo for the first 2 years. In this preclinical AD group (in the model), cognition is still very close to normal, but soluble and insoluble Aβ42 leave the region of normal values (Figure 3 , left panel). The BIF soluble Aβ42 leaves the normal region at 72 years in the preclinical AD group according to the model, that is why hypotheses 2.1 and 2.2 (Table 2) forecast no difference between treatment and placebo for over 2 years. The main limiting assumption of our hypotheses is the absence of irreversible toxicity or long‐term cellular processes. They significantly overestimate pharmacodynamic outcomes predicting almost 20 points of Adas‐cog improvement, which means complete healing. These and other additional hypotheses could be explored using this model and new emerging datasets at a future time. Soluble and insoluble Aβ42 dynamics are qualitatively similar in AD (Figure 3 a, upper and bottom rows) and leave the region of normal values simultaneously (70 years), so toxicity descriptions, in terms of dependence on soluble (in brain cells) or insoluble forms, may be similar.

Even one chosen hypothesis 2.1 (Figure 4 b) predicts the diversity of biomarker and clinical outcome behavior depending on the treatment protocol. For example, 50% RI leading to 50% CSF Aβ decline for both preclinical and mild AD causes distinct behaviors of other biomarkers: for the preclinical AD population (Figure 4 b) after 2 years of therapy, the model predicts neither difference from baseline in DVR nor Adas‐cog. In the mild AD case, 50% RI corresponds with much greater responses of DVR, and potential decline of Adas‐cog vs. placebo with time.

Responses predicted by simulation for the hypothesis 2.1–2.2., but not the other hypotheses, are qualitatively similar to results obtained in avagacestat clinical trial,39 as well as the bapineuzumab study.40 This would suggest that hypotheses of Aβ impact on cognition can be tested. The trend for a decrease in cognition in the avagacestat (and semagacestat46) trial would also bolster the concerns given hypotheses 2.1 and 2.2, in which cognition may be impaired if Aβ concentrations were altered to subphysiological levels. Additionally, if one were to extrapolate from Figure 4 a, hypotheses 2.1 and 2.2 may be tested when Aβ production is inhibited >50% over a shorter duration trial (<1 year). In this case, the predictions do not correspond to observation at 12 weeks, in which no cognitive decline is reported in the clinic (Figure 5). This may be due to model misspecification in the longitudinal nature of the low Aβ levels, other protein/macro molecules compensating for lower Aβ, or that the in vitro findings do not translate into the in vivo setting. Other adverse events may be related to GS inhibition but are not captured in this model. Results presented here do not deal with all possible hypotheses (e.g., we did not consider insoluble or oligomeric amyloid toxicity), nor do we show the extremes of Aβ inhibition (∼100% inhibition) or degradation, as neither consider specific features of different population (based, e.g., on ApoE genotype), but we have used this model and subsequent simulations to illustrate the approach of applying quantitative systems pharmacology (QSP) models for the purpose of hypothesis testing.

CONCLUSIONS

The constructed mechanistic model, verified with a number of data, allows insights into different aspects of translational and clinical studies and reveals key features, distinguishing healthy and AD populations. Physiological ranges and results of in vitro studies analyzed in this article suggest the possible benefit achieved by maintaining physiologic concentrations of different amyloid species. The model allows for the simulation of such appropriate regimens and may also be used for exploring other combination type regimens (inhibition of soluble Aβ and degradation of insoluble Aβ) in silico.

Direction, magnitude, and rate of biomarker changes predicted by the model provide the basis for planning and analysis of ongoing trials. As clinical studies for neurodegenerative diseases become more difficult to execute (due to cost and duration), systems pharmacology tools, like the one presented here, provide a greater understanding of the actual ability of a novel study to test the underlying hypothesis.

Supporting Information

Click here for additional data file.

Supplementary Figure S1

Click here for additional data file.

Supplementary Figure S2

Click here for additional data file.

Supplementary Figure S3

Click here for additional data file.

Supplementary Figure S4

Click here for additional data file.

Supplementary Figure S5

Click here for additional data file.

Supplementary Figure S6

Click here for additional data file.

Supplementary Figure S7

Click here for additional data file.

Supplementary Figure S8

Click here for additional data file.

Supplementary Figure S9

Click here for additional data file.

Supplementary Figure S10

Click here for additional data file.

Supplementary Figure S11

Click here for additional data file.

Supplementary Figure S12

Click here for additional data file.

Supplementary Figure S13

Click here for additional data file.

Supplementary Figure S14

Click here for additional data file.

Supplementary Figure S15

Click here for additional data file.

Supplementary Figure S16

Click here for additional data file.

Supplementary Table S1 Scheme of the parameter verification: groups of processes are highlighted by the color corresponding to the scheme in Supplementary Figure S4. Designations of stages of model calibration (F0‐F4) correspond to Table 1 of the main text. Different datasets have been used for each stage of calibration (see Table 1 of the main text).

Click here for additional data file.

Supplementary Table S2 Results of fitting and validation of the model for the healthy individuals: fitting of the data set for insoluble and soluble amyloid (F2 in Table 1 of the main text), for validation data set of the cerebrospinal fluid (CSF) concentrations was used (V2.2 in Table 1). Npar – number of parameters refitted for healthy individuals, Rss_fit and Rss_val – OFV for fitting and validation respectively according to formula (8), AIC – Akaike information criterion, RSS_av_fit and RSS_av_val – OFV averaged by number of points (Np) for fitting and validation dataset, Rel – relationship between averaged OFV. In Description column names of model variants are given with references to parameters of time dependences (see formuli (B.5‐B.11) that were fitted to describe difference between healthy and AD individuals.

Click here for additional data file.

Supplementary Table S3 Model parameter summary. Data used for fitting listed in Supplementary Figures S5–S7.

Click here for additional data file.

Supplementary Table S4 Values of parameters chosen on the basis of the clinical observation31 about Adas‐cog score change.

Click here for additional data file.

Supplement to Manuscript: Studying the progression of amyloid pathology and its therapy using translational longitudinal model of accumulation and distribution of amyloid beta DBSolve model code.

Click here for additional data file.

Supplementary Figure Legends

Click here for additional data file.