Home | About | Journals | Submit | Contact Us | Français |

**|**BMC Syst Biol**|**v.6; 2012**|**PMC3483162

Formats

Article sections

- Abstract
- Background
- Results and discussion
- Conclusions
- Methods
- Competing interests
- Authors' contributions
- Supplementary Material
- References

Authors

Related links

BMC Syst Biol. 2012; 6: 74.

Published online 2012 June 22. doi: 10.1186/1752-0509-6-74

PMCID: PMC3483162

Angelyn Lao,^{1} Vanessa Schmidt,^{#}^{2} Yvonne Schmitz,^{#}^{1} Thomas E Willnow,^{}^{2} and Olaf Wolkenhauer^{}^{1,}^{3}

Angelyn Lao: ed.kcotsor-inu@oal.nylegna; Vanessa Schmidt: ed.nilreb-cdm@tdimhcs.assenav; Yvonne Schmitz: ed.kcotsor-inu.kitamrofni@ztimhcs; Thomas E Willnow: ed.nilreb-cdm@wonlliw; Olaf Wolkenhauer: ed.kcotsor-inu@reuahneklow.falo

Received 2012 January 31; Accepted 2012 May 14.

Copyright ©2012 Lao et al.; licensee BioMed Central Ltd.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (
http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

This article has been cited by other articles in PMC.

Proteolytic breakdown of the amyloid precursor protein (APP) by secretases is a complex cellular process that results in formation of neurotoxic Aβ peptides, causative of neurodegeneration in Alzheimer’s disease (AD). Processing involves monomeric and dimeric forms of APP that traffic through distinct cellular compartments where the various secretases reside. Amyloidogenic processing is also influenced by modifiers such as sorting receptor-related protein (SORLA), an inhibitor of APP breakdown and major AD risk factor.

In this study, we developed a multi-compartment model to simulate the complexity of APP processing in neurons and to accurately describe the effects of SORLA on these processes. Based on dose–response data, our study concludes that SORLA specifically impairs processing of APP dimers, the preferred secretase substrate. In addition, SORLA alters the dynamic behavior of β-secretase, the enzyme responsible for the initial step in the amyloidogenic processing cascade.

Our multi-compartment model represents a major conceptual advance over single-compartment models previously used to simulate APP processing; and it identified APP dimers and β-secretase as the two distinct targets of the inhibitory action of SORLA in Alzheimer’s disease.

The amyloid precursor protein (APP) is a type-1 membrane protein expressed in neurons, which is closely linked to the etiology and pathology of Alzheimer’s disease (AD) [1]. APP undergoes two mutually exclusive processing pathways resulting in the formation of multiple soluble and membrane-associated fragments from this precursor polypeptide. Of particular relevance to AD is the amyloidogenic pathway whereby APP is first cleaved by β-secretase and subsequently by γ-secretase to produce the amyloid-β peptide (Aβ), a 40 to 42 amino acid fragment derived from part of the extracellular and the transmembrane domains of APP. According to the amyloid hypothesis, neurotoxic oligomers and senile plaques formed by Aβ cause neuronal dysfunction and cell loss in AD [2,3]. In the alternative pathway, APP is first cleaved by α-secretase, instead of β-secretase, resulting in the destruction of the Aβ peptide sequence in APP (non-amyloidogenic pathway). Adding to the complexity of APP processing is the distinct trafficking route of the precursor through intracellular compartments where the various secretases reside [4,5]. Thus, newly synthesized APP molecules move through the constitutive secretory pathway from the trans-Golgi network (TGN) to the cell surface where most are subjected to non-amyloidogenic processing by α-secretase. However, approximately 10% of the precursors remain intact and internalize from the cell surface into endosomal compartments where amyloidogenic processing is initiated by β-secretase cleavage.

In recent years, much attention has been focused on the analysis of factors that influence APP processing and that may contribute to the elevated Aβ levels seen in patients with AD. One such modifier is SORLA, a 250-kDa type-1 membrane glycoprotein widely expressed in neurons in the brain
[6,7]. It is a member of a family of mammalian proteins that share a structural similarity with the vacuolar protein sorting 10 protein (VPS10p), a y*east* sorting receptor that transports carboxypeptidase Y from the TGN to the vacuole
[8]. SORLA is proposed to act as a retention factor for APP in the TGN, preventing the release of precursor molecules into the processing pathways. Consequently, over-expression of SORLA in neurons prevents the targeting of APP from TGN to the cell surface and to endosomes and reduces the production of Aβ peptides
[9-11]. The importance of SORLA for AD is further supported by low levels of receptor expression seen in patients suffering from the disease
[12,13].

In a first approach to simulate amyloidogenic processing in AD, Schmidt and colleagues [14] developed a single-compartment model to describe APP processing. While this model had been valuable to establish the kinetics of amyloidogenic processing and the quantitative contribution of SORLA to this pathway, this single-compartment model fell short of accurately describing the complexity of APP processing in cells. It remained unclear to what extent SORLA may affect APP monomer versus dimer processing and in what compartment of the cell its activity may be most relevant. Also, a possible influence of SORLA on the dynamics of β-secretase remained unclear. Such an effect had been postulated previously based on studies in cultured cells [15].

To answer these questions, we established a multi-compartment model that represents APP processing in both its monomeric and dimeric forms. The formalism of this model is developed to integrate experimental evidences from previous biochemical and cell biological studies [9-11,13,15-17]. We combined our multi-compartment model with the recent dose–response data of APP and soluble APP products by Schmidt and colleagues [14]. The data were used to estimate the parameter values of our model. Using our multi-compartment model, we (i) established the activity distribution of APP in various compartments, and (ii) traced the activity distribution of APP, α-secretase, β-secretase and SORLA in the monomeric and dimeric processing of APP. Our simulation results showed that the decrease in total APP processing is primarily due to the influence of SORLA on APP dimer processing. Moreover, the simulations of our multi-compartment model demonstrated how SORLA alters the dynamical behavior of β-secretase, providing new insights into the mechanism of action of this important AD risk factor.

Probably more than any other major disease entity, AD is a pathological processes influenced by subtle quantitative changes in protein concentration and activity. Thus, common approaches in experimental AD research, using protein overexpression or gene-inactivation, are inadequate to study the effects of incremental changes in target protein levels on risk of neurodegeneration.

In our previous study [14], we have undertaken the first attempt to approach risk factors in AD through quantitative modeling. To do so, we have simulated the quantitative contribution of SORLA to proteolytic processing of APP, a central pathway in AD. We have chosen SORLA as a target for simulation because it represents one of the major genetic risk factors in AD. More importantly, solid experimental evidence had established the molecular mechanism of SORLA action, acting as an intracellular sorting receptor for APP that prevents proteolytic breakdown of the precursor protein into neurotoxic Aβ peptides. In [14], we have been able mathematically confirm hypotheses, derived from prior experimental work. In particular, we have confirmed the strict linear relationship between SORLA concentrations and efficiency of APP processing, and we have uncovered the ability of SORLA to prevent dimerization of APP, preventing the formation of high-affinity substrates for secretases.

While our initial study has been met with great enthusiasm in the field, it clearly falls short of addressing major aspects of SORLA activity in the cell biology of AD. Thus, for sake of simplification, our earlier study assumed a single-compartment model for simulation of the affects of SORLA levels on APP processing rates. Accordingly, it ignores the fact that APP follows a complex intracellular trafficking pathway whereby this protein moves between the TGN, cell surface, and endosomes where the various interacting proteins reside. In fact, it has the ability to show how SORLA affects APP transport between various cell compartments in neurons that initially sparked interest in this protein.

For the present work, a single-compartment model, describing the influence of SORLA in APP processing [14], was extended into a multi-compartmental model. The extended model addresses the important aspect of the cell biology of SORLA by assuming a three-compartment model that is based on experimental data. The biochemical network illustrating this multi-compartmental model can be found in Figure Figure1.1. The notation that is used in this network is described in Additional file 1: Table S1.

It is likely that there are many other proteins contribute to the processing of APP and the generation of neurotoxic Aβ peptides. However, unlike many proposed AD risk factors, the mechanism of action for SORLA has been established in numerous studies in cell cultures, in animal models, and even in patients providing a solid base for theoretical simulations. In particular, we have specifically addressed the caveat that this model focuses on pathways related to SORLA action, and that further studies will be required to sequentially add more risk factors to this model. Such approaches will require a profound understanding of the function of such risk factors; - an endeavor that clearly exceeds the scope of the present manuscript.

The choice of the compartments considered in Figure Figure11 was based on the different locations where APP was shown to interact with SORLA, with α-, and with β-secretases. The corresponding three compartments are the TGN, the cell surface and the endosomes [9,15,20], respectively. Note that the transport of APP among these compartments indirectly interconnects these three compartments to one another. As SORLA affects the initial cleavage of APP by α- and β-secretases [11], the rate limiting steps that determine the extent of amyloidogenic processing, further processing steps involving γ-secretase were not included in this model.

In order to accommodate the monomeric and dimeric forms of APP, each compartment was further divided into two subcompartments (Figure
(Figure1):1): a “red” subcompartment for APP monomer processing and a “green” one for APP dimer processing. Notice that the monomeric forms of APP, SORLA, α-secretase, and β-secretase, within the two subcompartments, were annotated differently: *APP*_{G1}, *SORLA*_{G1}, *α*_{1}, and *β*_{1} for monomer processing, and *APP*_{G2}, *SORLA*_{G2}, *α*_{2}, and *β*_{2} for dimer processing. Even so, the components from the two subcompartments were linked to each other via *APP*_{init}, *α*_{init}, and *β*_{init}. Moreover, *APP*_{G2}, *α*_{2}, and *β*_{2} undergo dimerization before the beginning of APP dimer processing. That is, two *APP*_{G2}, *α*_{2}, and *β*_{2} monomers dimerize in order to give their corresponding dimeric forms. Conversely, these dimers can dissociate to generate their respective monomers. Note that subscript ‘1’ was assigned to the reactants and products in monomer processing while subscript ‘2’ for those in dimer processing. In addition, we used subscripts ‘G’, ‘CS’, and ‘E’ for APP in TGN, at the cell surface and in the endosomes, respectively.

Up to this point, we had described the different forms of APP, α-secretase, and β-secretase in the diverse compartments, prior to the beginning of APP processing. Because SORLA interacts with APP in a 1:1 stochiometric complex
[9,16], the model described how SORLA strictly interacts with APP-monomers (but not dimers) to form an APP-SORLA complex. Consequently, this interaction is responsible for the diminished amount of APP-monomers (*APP*_{G1}) and APP-dimers (*APP*_{G2d}) transported from the TGN to the cell surface. This interaction decreases the amount of APP-monomers (*APP*_{CS1}) and APP-dimers (*APP*_{CS2d}) ending up in the endosomes as *APP*_{E1} and *APP*_{E2d}. Moreover, in order to determine whether SORLA will have a similar influence on the monomer and dimer processing, the binding affinity assigned to *APP*_{G1}*SORLA*_{G1} in monomer processing is different to that of *APP*_{G2}*SORLA*_{G2} in dimer processing.

After the interaction of SORLA and APP in the TGN, the remaining *APP*_{G1} and *APP*_{G2d} are transported to the cell surface where APP processing begins within the non-amyloidogenic pathway. Then, a small part of *APP*_{CS1} and *APP*_{CS2d}, which is not cleaved by α-secretase, are further transported from the cell surface to the endosomes, where the amyloidogenic pathway takes over. Notably, the interaction of APP and α-secretase at the cell surface leads to the formation of non-amyloidogenic products like sAPPα and C83; whereas the interaction of APP and β-secretase in the endosomes yields to the amyloidogenic products such as sAPPβ and C99. Our model was established in such a way that the dimeric form of secretases act only on the dimeric form of APP and the monomeric form of secretases act only on the monomeric form of APP.

The biochemical network (Figure (Figure1)1) that we established, was translated into a system of ordinary differential equations (ODEs), describing temporal changes of molecular numbers for the network components as a function of interaction and cleavage processes. The model equations, their reduction, and a series of steps involved in model simulations are presented in the Materials and Methods section.

With the multi-compartmental model, we showed in Figure Figure22 the corresponding model simulations for various APP products, namely, the products produced in monomer, in dimer, and in both processing pathways. The simulations of the parameterized mathematical model are in good agreement with recent experimental data by Schmidt and colleagues [14] (Figure (Figure22A-D).

In the absence of SORLA, the sigmoidal curve that is characteristic for products produced in dimer processing (green lines in Figure Figure2A2A and B) has a strong impact on the sum of the products produced in monomer and in dimer processing pathways (black lines in Figure Figure2A2A and B). As such, it very well describes the experimental data sets for sAPPα and sAPPβ (black dots in Figure Figure2A2A and B, respectively).

Surprisingly, in the presence of SORLA, one observes from the simulations a significant decrease in the products produced in dimer processing (green lines in Figure
Figure2C2C and D) as compared to those in monomer processing (red lines in Figure
Figure2C2C and D). In particular, the analysis showed that at a high level of SORLA activity (i.e. 100% of SORLA_{Tot} where SORLA_{Tot} equals 2.43 x 10^{5} fmol), there is obviously more APP bound to SORLA in dimer processing (Figure
(Figure3B)3B) than in monomer processing (Figure
(Figure33A).

Taken together, our simulations shown in Figure Figure22 and Figure Figure3,3, strongly supported the hypothesis whereby SORLA prevents oligomerization of APP, thereby having a bigger impact on the products produced in dimer processing than in monomer processing.

Up to this point, we only showed simulations of our model in the two most extreme scenarios: with no (Figure
(Figure2A2A and B) or high levels of SORLA activity (Figure
(Figure2C2C and D). However, subtle alterations of SORLA concentration are likely to be more relevant for the determination of its influence in APP processing pathways. Accordingly, we adapted our multi-compartment model to intermediate concentrations of SORLA. As shown in Figures
Figures44,
,55,
,66 and
and77 the simulations are all in dependence of three intermediate SORLA expression levels, namely, 3%, 12%, and 30% of SORLA_{Tot}.

Remarkably, we observed in Figure Figure44 that the simulations in dependence of the three intermediate SORLA expression levels are either “spread” (as in Figure Figure4A4A and Figure Figure4D)4D) or “clustered” (as in Figure Figure4B4B and Figure Figure4C)4C) into the two most extreme scenarios of SORLA concentration. This came as a surprised because the dose–response kinetics of total sAPPα production in dependence of the intermediate SORLA expression levels (Figure (Figure4A)4A) is expected to be “clustered” like that of sAPPβ (Figure (Figure4B).4B). Likewise in the case of the amount of APP bound to SORLA in monomer (Figure (Figure4C)4C) and in dimer processing (Figure (Figure4D).4D). We say that the simulations are “clustered” when

$${\text{X}}_{Y}\approx ({X}_{100\%}-{X}_{0\%})\cdot Y+{X}_{0\%}$$

where *Y*={3%, 12%, 30%}, and *X* denotes the amount of concentration at a given percentage value of SORLA_{Tot} that is specified by its subscript. Otherwise, we say that the simulations are “spread”.

Next, we investigated what leads to the observation made in Figure Figure4,4, in dependence of the intermediate SORLA expression levels.

First, we analyzed the simulations of the influence of intermediate levels of SORLA on APP processing on the amount of α-secretase (Figure (Figure5A-F)5A-F) and β-secretase (Figure (Figure5)5) concentration. In Figure Figure5,5, the term “used” refers to the complex formation of the secretases and APP, while the term “free” refers to the secretases that are not bound in a complex.

The total amount of α-secretase and the total amount of β-secretase were assumed to be constant (depicted by the black lines in Figure
Figure5E-F5E-F and Figure
Figure5K-L,5K-L, respectively). Due to the conservation law assumption, the total amount of each secretase in each subcompartment is conserved (i.e. α_{monomer} and β_{monomer} depicted by red lines in Figure
Figure5F5F and Figure
Figure5L;5L; α_{dimer} and β_{dimer} depicted by green lines in Figure
Figure5F5F and Figure
Figure5L).5L). Consequently, the total amount of each secretase in the whole system was thus also conserved (α_{Tot} and β_{Tot} shown by the black lines in Figure
Figure5E-F5E-F and Figure
Figure5K-L,5K-L, respectively).

The simulations of the influence of intermediate levels of SORLA on APP processing on the amount of α-secretase (Figure (Figure5A-F)5A-F) concentration showed that (i) there are more α-secretases that were used (Figure (Figure5C)5C) than left free (Figure (Figure5A)5A) in monomer processing, (ii) there are more α-secretases that are left free (Figure (Figure5B)5B) than used (Figure (Figure5D)5D) in the dimer processing, (iii) the total amount of α-secretase that is free and used (blue and orange lines in Figure Figure5E,5E, respectively) is dominated by the corresponding amount of α-secretase concentration in dimer (Figure (Figure5B)5B) and in monomer processing (Figure (Figure5C),5C), (iv) SORLA influences the amount of α-secretase concentration in dimer processing (Figure (Figure5B5B and and5D),5D), but not those in monomer processing (Figure (Figure5A5A and and5C),5C), and (v) its simulations in dependence of the three intermediate SORLA expression levels (Figure (Figure5D)5D) is consistent to that of dose–response kinetics of total sAPPα production (Figure (Figure44A).

The significant difference in the free (Figure (Figure5B)5B) and used (Figure (Figure5D)5D) amounts of α-secretase in dimer processing is a consequence of the large amount of α-secretase used in monomer processing (shown in Figure Figure5C).5C). As the total amount of the APP concentration increases (from 0 nM to 400 nM), the amount of α-secretase, free in dimer processing, (Figure (Figure5B)5B) decreases, while the amount of α-secretase used in monomer processing (Figure (Figure5C)5C) increases. As the amount of SORLA concentration increases, the curves representing the secretases move from solid to dashed lines. SORLA does affect α-secretase in dimer processing (Figure (Figure5B5B and and5D):5D): those used in dimer processing decreases (Figure (Figure5D),5D), while those that are free in dimer processing increases (Figure (Figure5B).5B). In the later figure, the increase is not obvious because the amount of change is so small as compared to the concentration values of α-secretase.

As for the influence of intermediate levels of SORLA on APP processing on the amount of β-secretase (Figure (Figure5G-L)5G-L) concentration, the simulations showed that (i) there are more β-secretases that are left free (Figure (Figure5G)5G) than used (Figure (Figure5I)5I) in monomer processing, (ii) SORLA has no influence on β-secretase in monomer processing (Figure (Figure5G5G and Figure Figure5I),5I), (iii) SORLA alters the dynamical behaviors of β-secretase in dimer processing (Figure (Figure5H5H and Figure Figure5J),5J), (iv) the total amount of β-secretase that is free and used (blue and orange lines in Figure Figure5K,5K, respectively) is dominated by the amount of β-secretase concentration in dimer processing (Figure (Figure5H5H and Figure Figure5J,5J, correspondingly), and (v) its simulations in dependence of the three intermediate SORLA expression levels (Figure (Figure5H5H and Figure Figure5J)5J) is consistent to that of dose–response kinetics of total sAPPβ production (Figure (Figure4B).4B). The curves for beta-secretase with SORLA (dashed lines in Figure Figure5H)5H) are greater in values as compared to those without SORLA (solid line in Figure Figure5H),5H), as a consequence of SORLA’s influence on beta-secretase that is used in dimer processing (Figure (Figure55J).

When a comparison is made between the total amount of α- and β-secretase concentration that is free (blue lines in Figure Figure5E5E and Figure Figure5K)5K) and used (orange lines in Figure Figure5E5E and Figure Figure5K)5K) in dependence of the three intermediate SORLA expression levels, we observed that the total amount of β-secretase concentration for both free and used deviated (Figure (Figure5K),5K), which was not the case for α-secretase (Figure (Figure5E).5E). This observation suggested that SORLA is indirectly affecting the dynamics of β-secretase but not that of α-secretase. This result supports the hypothesis presented by Schmidt el al. [14]: “the global–local estimation of the parameter values in the model suggested a yet unidentified biological process whereby SORLA might indirectly affect the β-secretase, but not with the α-secretase”. The present result therefore clarifies what was unidentified in our previous study [14].

With SORLA concentration greater than the estimated total amount of SORLA concentration (i.e. SORLA_{Tot}=2.43 x 10^{5} fmol), we arrived at Figure S1 shown in Additional file
1: Figures S1D and S1J show that for a very large amount of SORLA_{Tot} (greater than 1 x SORLA_{Tot} for α-secretase and greater than 10 x SORLA_{Tot} for β-secretase), the amount of α- and β-secretase are barely “used”. Consequently, the amount of α- (Figure B) and β-secretase (Figure H) are all “free” in dimer processing, and there will be no sAPP products produced in dimer processing.

We also investigated the amount of APP concentrations that is either free or used, in monomer or in dimer processing, and which is in the TGN, at the cell surface or in the endosomes (Figure (Figure6).6). The term “used” refers to the complex formation of (i) APP and SORLA in the TGN, (ii) APP and α-secretase at the cell surface, and (iii) APP and β-secretase in the endosomes. Wherein, the term “free” refers to the APP that is not bound in the respective compartments.

First, we showed the simulations of the amount of APP concentrations that is free or used in monomer and in dimer processing. The simulations under dimer processing showed that the amount of APP concentrations that is free or used in each compartment were significantly affected by the presence of SORLA (last two columns of Figure Figure6:6: Figure Figure6D-E,6D-E, Figure Figure6H-I,6H-I, Figure Figure6L-M,6L-M, and Figure Figure6P-Q),6P-Q), as compared to those under monomer processing (first three columns of Figure Figure6:6: Figure Figure6A-C,6A-C, Figure Figure6F-G,6F-G, Figure Figure6J-K,6J-K, and Figure Figure6N-O).6N-O). In particular, one observes from the simulations that the amount of APP concentrations that is used to bind with SORLA in dimer processing of the TGN tremendously increases from 0 M to at most 300 nM (Figure (Figure6E),6E), wherein those in monomer processing are so small that they can be neglected (Figure (Figure6C).6C). Consequently, SORLA decreases the amount of APP concentrations that is free or used at the cell surface and in the endosomes (Figure (Figure6H-I6H-I and Figure Figure6L-M,6L-M, respectively). Also, the total amount of APP concentrations in dimer processing is dominated by the total amount of free APP in the absence of SORLA and by the total amount of used APP in the presence of SORLA (depicted by the two outermost lines in Figure Figure6P6P and Figure Figure66Q).

Next, in each compartment, the simulations for the total amount of APP concentrations that is free, used, in monomer processing, or in dimer processing, are shown in Figure
Figure7.7. Consistent to our previous observation (Figure
(Figure6),6), the simulations for total amount of APP concentrations in monomer processing for the three different compartments (3^{rd} column of Figure
Figure7)7) were not influenced by SORLA, while those in dimer processing were affected by the presence of SORLA (4^{th} column of Figure
Figure7).7). Moreover, the simulations, in the first two columns of Figure
Figure7,7, also showed that the presence of SORLA in the TGN decreases the total amount of free APP (Figure
(Figure7P),7P), and increases the total amount of used APP (Figure
(Figure7Q).7Q). In particular to the total amount of used APP under the influence of SORLA, it is (i) enormously increased in the TGN (Figure
(Figure7B),7B), (ii) not affected at the cell surface (Figure
(Figure7G),7G), and (iii) reduced by at most half in the endosomes (Figure
(Figure7L).7L). Taken together, the presence of SORLA increases the total amount of APP concentrations in the TGN (Figure
(Figure7E),7E), and subsequently decreases the total amount of APP concentrations at the cell surface (Figure
(Figure7J)7J) and in the endosomes (Figure
(Figure77O).

The simulations for the total amount of APP concentrations in monomer processing (Figure (Figure7R),7R), in dimer processing (Figure (Figure7S),7S), and in both monomer and dimer processing (Figure (Figure7T)7T) show that a conservation law was assumed for APP in monomer and in dimer processing. Above all, one observes that there are more APP concentrations in dimer processing (Figure (Figure7S)7S) than in monomer processing (Figure (Figure77R).

As noted in the subsection, *Intermediate levels of SORLA*, the simulations show that SORLA expression levels are either “spread” (Figure
(Figure4A)4A) or “clustered” (Figure
(Figure4B).4B). This is most likely due to the effect of SORLA on the processing of APP dimer. With respect to the total amount of APP, the amount of APP concentrations (Figure
(Figure6I)6I) and α-secretase concentrations (Figure
(Figure5D)5D) that are “used” at the cell surface in dimer processing “spread”. Considering the relevance of APP and α-secretase at the cell surface to the production of sAPPα, the observations thus suggest the “spread” observed in Figure
Figure4A4A for sAPPα. Similarly for the “clustering” observed in Figure
Figure4B4B for sAPPβ, it is a consequence of the “clustering” that is observed on APP (Figure
(Figure6M)6M) and β-secretase (Figure
(Figure5J)5J) that are “used” in the endosome in dimer processing, which are relevant in producing sAPPβ. Moreover, the change from “spread” at the cell surface (Figure
(Figure4A)4A) to “clustered” in the endosome (Figure
(Figure4B)4B) is probably due to the indirect influence of SORLA on the dynamical behavior of β-secretase that is observed in Figure
Figure55.

Lastly, in Figure
Figure8,8, are given simulations of the influence of SORLA on APP processing into sAPPα (Figure
(Figure8A8A and Figure
Figure8C)8C) and sAPPβ (Figure
(Figure8B8B and Figure
Figure8D).8D). The simulations show that the switch from preferred dimer-to-monomer processing is observed at 25% of SORLA_{Tot} for α-secretase (Figure
(Figure8A)8A) and at 3% of SORLA_{Tot} for β-secretase (Figure
(Figure8B),8B), where SORLA_{Tot} equals 2.43 x 10^{5} fmol. In agreement with the study performed by Schmidt and colleagues
[14] previously, we therefore find that the switch from cooperative (dimer) to less efficient non-cooperative (monomer) processing occurs at small amount of SORLA concentration. Moreover, the end product obtained from monomer processing dominates the total amount of end product at 145% of SORLA_{Tot} for α-secretase (Figure
(Figure8C)8C) and at 150% of SORLA_{Tot} for β-secretase (Figure
(Figure8D).8D). In connection to what we observed in Figure
Figure44 for the simulations of the influence of intermediate levels of SORLA on APP processing into sAPPα (Figure
(Figure4A)4A) and sAPPβ (Figure
(Figure4B),4B), these two sets of results (Figure
(Figure44 and Figure
Figure8)8) suggest that SORLA reduces the products produced in non-amyloidogenic and amyloidogenic pathways of APP processing at different rate.

Our multi-compartment model is an extension of the single-compartment model that was established by Schmidt and colleagues [14] previously. To our knowledge, this is the first multi-compartmental model developed to analyze APP processing in the context of Alzheimer’s disease. In addition, our model represents the regulated trafficking of APP by SORLA through the intracellular compartments, which critically affects amyloidogenic and non-amyloidogenic processing pathways [20]. Our model was established to answer questions that arose from a study based on a single-compartment model [14].

The first question that emerged concerned the relative contributions of SORLA to monomer and in dimer processing. In the study conducted by Schmidt and colleagues [14], we showed that SORLA influenced the combined products obtained by monomer and dimer processing. However, limited by the structure of the single-compartment model, we were not able to investigate the relative contribution of SORLA in monomer and in dimer processing. Herein, using our multi-compartment model, we showed that the decrease in total amounts of sAPPα and sAPPβ is mainly due to the influence of SORLA in dimer processing. This observation confirms previous hypothesis that SORLA, prevent oligomerization of APP, eliminating the preferred substrates for secretases.

The second question was how does SORLA affect the dynamics of β-secretase? In the study conducted by Schmidt and colleagues [14], it was suggested that there is an indirect effect of the SORLA receptor on the β-secretase, which contributes to the regulation of amyloidogenic processing in the context of an intact cell. However, in order for the single-compartment model to closely resemble the experimental data, the model required a local parameter estimate for β-secretase activity in the presence or absence of SORLA. Through our multi-compartment model, where all parameters are estimated globally, we now confirmed that SORLA affects the interaction between APP and β-secretase, but not that of APP with α-secretase. A previous study suggested that SORLA directly interacts with β-secretase, preventing access of the enzyme to its substrate APP [15]. While our simulations confirm an important influence of SORLA on β-secretase, this influence may also be indirect, for example by effecting trafficking of cofactors essential for enzyme activity. An indirect effect of SORLA is in line with findings that the receptor does not impair β-secretase activity in cell-free assays [14].

In addition, we investigated the regulated trafficking of APP by SORLA in monomer and dimer processing, considering several cellular compartments, including TGN, cell surface, and endosomes. Simulations of our multi-compartment model showed that SORLA increases the total amount of APP concentrations in the TGN (Figure (Figure7E)7E) and subsequently decreases the total amount of APP at the cell surface (Figure (Figure7J)7J) and endosomes (Figure (Figure7O).7O). In agreement with Andersen and Willnow [21], this result suggests that an over-expression of SORLA prevents the localization of APP from the TGN to the cell surface and to the endosomes, whereby an over-expression of SORLA decreases the products produced in the amyloidogenic and non-amyloidogenic pathways of APP processing. Furthermore, our study confirmed that SORLA is more influential in dimer than in monomer processing. This observation is in line with our previous model that APP dimers represent the preferred substrate for α- and β-secretase as they enable cooperativity in substrate binding [14]. Taken together, data obtained both in single and in multi-compartment models strongly suggested that depletion of APP dimer processing represents a major molecular mechanism in the pathology of Alzheimer’s disease.

Our multi-compartment model was used to simulate pathological situations involving APP under different level of SORLA concentration. Our model can also be used as a kinetic-dynamic model to study the effects of SORLA on α- and β- secretase. Moreover, we observed that as the amount of SORLA concentration increases, there is a relatively large decrease in the production rate of sAPPα as compared to that of sAPPβ (Figure (Figure44 and Figure Figure88).

Using this refined model together with the chosen set of estimated parameter values (shown in Additional file
1: Table S4), our results suggest the following biological implications of SORLA: (1) Decrease in total amounts of sAPP products is mainly due to the large amount of SORLA concentration in dimer processing (2.43 x 10^{5} fmol), and not to the small amount of SORLA concentration in monomer processing (1.23 x 10^{1} fmol). (2) SORLA indirectly affects the dynamical behavior of the β-secretase but not that of α-secretase. The receptor targets β-secretase, the enzyme responsible for initial amyloidogenic cleavage. This finding represents a major conceptual advance in our understanding of the complex processes in APP processing and supports initial biochemical data that SORLA can bind to β-secretase
[15]. (3) SORLA is more influential in dimer processing that in monomer processing, which confirmed our initial hypothesis that blockade of APP dimerization is an important aspect of SORLA action on AD.

In future studies, we will extend this model by including additional cleavage activity by γ-secretase in monomer and in dimer processing. Cleavage of APP by γ-secretase leads to the formation of Aβ peptides, which is the main pathology of Alzheimer’s disease. The ultimate goal is therefore to establish a model that will test the potential effects of SORLA on APP processing in the context of AD therapy.

Based on the biochemical network (Figure (Figure1),1), we established ODEs that describe temporal changes of molecular numbers for the network components as a function of interaction and cleavage processes, such that the changes with large numbers of molecules can be assumed to be smooth. The complete formulation of the model, realized as set of ODEs, can be found in the Additional file 1. Also, the notations used in the equations are described in Additional file 1: Table S2.

Recall that the components in our multi-compartment model, separated by the three different compartments and two different subcompartments, are annotated differently. Distinctly labeling them allows differentiation and comparison of the results generated from these multi-compartment. Consequently, we could distinguish the level of influence of SORLA between the monomeric and the dimeric form of APP processing, and also determine in which form of APP processing is SORLA more influential.

Herein, we showed a series of assumptions that allowed the reduction of the equations.

Firstly, we started with the transportation of APP, α-, and β-secretase among the three compartments and between the two subcompartments. Recall that the (i) APP-monomers and APP-dimers are transported from the TGN to the cell surface, and are then further transported to the endosomes, (ii) monomeric forms of APP, SORLA, α-secretase, and β-secretase, within the two subcompartments, are annotated differently, and (iii) components within the two subcompartments are linked to each other via *APP*_{init}, *α*_{init}, and *β*_{init}. These properties are reflected by introducing the following biochemical reactions: For APP,

$$\begin{array}{c}AP{P}_{\mathit{init}}\stackrel{{K}_{G1}}{\to}AP{P}_{G1}\\ AP{P}_{\mathit{init}}\stackrel{{K}_{G2}}{\to}AP{P}_{G2}\\ AP{P}_{G1}\stackrel{{K}_{CS1}}{\to}AP{P}_{CS1}\\ AP{P}_{G2d}\stackrel{{K}_{CS2}}{\to}AP{P}_{CS2d}\end{array}$$

and for the secretases,

$$\begin{array}{c}{\alpha}_{\mathit{init}}\stackrel{{K}_{C1}}{\to}{\alpha}_{1}\\ {\alpha}_{\mathit{init}}\stackrel{{K}_{C2}}{\to}{\alpha}_{2}\\ {\beta}_{\mathit{init}}\stackrel{{K}_{B1}}{\to}{\beta}_{1}\\ {\beta}_{\mathit{init}}\stackrel{{K}_{B2}}{\to}{\beta}_{2}\end{array}$$

where the second reactant is assumed to be in quasi-equilibrium with the first reactant. Without loss of generality, the concentration of the second reactant is related to the first reactant by an ordinary equilibrium expression, such as

$$\left(\begin{array}{c}AP{P}_{G1}={K}_{G1}\cdot AP{P}_{\mathit{init}}\\ AP{P}_{G2}={K}_{G2}\cdot AP{P}_{\mathit{init}}\\ AP{P}_{CS1}={K}_{CS1}\cdot AP{P}_{G1}\\ AP{P}_{CS2d}={K}_{CS2}\cdot AP{P}_{G2d}\\ {\alpha}_{1}={K}_{C1}\cdot {\alpha}_{\mathit{init}}\\ {\alpha}_{2}={K}_{C2}\cdot {\alpha}_{\mathit{init}}\\ {\beta}_{1}={K}_{B1}\cdot {\beta}_{\mathit{init}}\\ {\beta}_{2}={K}_{B2}\cdot {\beta}_{\mathit{init}}\end{array}\right\}$$

(1)

where ${K}_{G1}={k}_{g1}/{k}_{-g1}$, ${K}_{G2}={k}_{g2}/{k}_{-g2}$, ${K}_{CS1}={k}_{cs1}/{k}_{-cs1}$, ${K}_{CS2}={k}_{cs2}/{k}_{-cs2}$, ${K}_{C1}={k}_{c1}/{k}_{-c1}$, ${K}_{C2}={k}_{c2}/{k}_{-c2}$, ${K}_{B1}={k}_{b1}/{k}_{-b1}$, and ${K}_{B2}={k}_{b2}/{k}_{-b2}$.

As for

$$\begin{array}{c}AP{P}_{CS1}\begin{array}{c}\hfill \stackrel{{k}_{e1}}{\to}\hfill \\ \hfill \underset{{k}_{-e1}}{\leftarrow}\hfill \end{array}AP{P}_{E1}\\ AP{P}_{CS2d}\begin{array}{c}\hfill \stackrel{{k}_{e2}}{\to}\hfill \\ \hfill \underset{{k}_{-e2}}{\leftarrow}\hfill \end{array}AP{P}_{E2d}\end{array}$$

the ratio of the association constant is taken into consideration. This assumption permits the two reactants to be related to each other by an association constant, such that

$$\left(\begin{array}{c}AP{P}_{E1}={K}_{E1}\cdot AP{P}_{CS1}\\ AP{P}_{E2d}={K}_{E2}\cdot AP{P}_{CS2d}\end{array}\right\}$$

(2)

where ${K}_{E1}={k}_{e1}/{k}_{-e1}$ and ${K}_{E2}={k}_{e2}/{k}_{-e2}$.

Secondly, the monomeric forms of APP, α-, and β-secretase in dimer processing undergo dimerization. In other words, two monomeric forms of APP, α-, or β-secretase are dimerized, and a dimeric form of APP, α-, or β-secretase is dissociated, as shown by the reactions below:

$$\begin{array}{cc}\hfill AP{P}_{G2}+AP{P}_{G2}\hfill & \hfill \begin{array}{c}\hfill \stackrel{{k}_{g3}}{\to}\hfill \\ \hfill \underset{{k}_{-g3}}{\leftarrow}\hfill \end{array}AP{P}_{G2d}\hfill \\ \hfill {\alpha}_{2}+{\alpha}_{2}\hfill & \hfill \begin{array}{c}\hfill \stackrel{{k}_{c3}}{\to}\hfill \\ \hfill \underset{{k}_{-c3}}{\leftarrow}\hfill \end{array}{\alpha}_{2d}\hfill \\ \hfill {\beta}_{2}+{\beta}_{2}\hfill & \hfill \begin{array}{c}\hfill \stackrel{{k}_{b3}}{\to}\hfill \\ \hfill \underset{{k}_{-b3}}{\leftarrow}\hfill \end{array}{\beta}_{2d}\hfill \end{array}$$

As we took into account the association constant of each reaction above, it allows us to have the following representation of the equation:

$$\left(\begin{array}{c}AP{P}_{G2d}={K}_{G3}\cdot AP{P}_{G2}^{2}\\ {\alpha}_{2d}={K}_{C3}\cdot {\alpha}_{2}^{2}\\ {\beta}_{2d}={K}_{B3}\cdot {\beta}_{2}^{2}\end{array}\right\}$$

(3)

where ${K}_{G3}={k}_{g3}/{k}_{-g3}$, ${K}_{C3}={k}_{c3}/{k}_{-c3}$, and ${K}_{B3}={k}_{b3}/{k}_{-b3}$.

Up to this point, we had established equations that describe the relationship of APP, α-, or β-secretase, which are defined in the three different compartments.

Thirdly, at the cell surface, APP interacts with α-secretase through the following reactions

$$\begin{array}{c}AP{P}_{CS1}+{\alpha}_{1}\begin{array}{c}\hfill \stackrel{{k}_{5}}{\to}\hfill \\ \hfill \underset{{k}_{-5}}{\leftarrow}\hfill \end{array}{C}_{APPCS1\alpha 1}\stackrel{{k}_{6}}{\to}sAPP{\alpha}_{1}+C{83}_{1}+{\alpha}_{1}\\ AP{P}_{CS2d}+{\alpha}_{2d}\begin{array}{c}\hfill \stackrel{{k}_{51}}{\to}\hfill \\ \hfill \underset{{k}_{-51}}{\leftarrow}\hfill \end{array}{C}_{APPCS2d\alpha 2d}\stackrel{{k}_{61}}{\to}2\cdot sAPP{\alpha}_{2}+C{83}_{2d}+{\alpha}_{2d}\end{array}$$

Whereas, in the endosomes, the following reactions take place:

$$\begin{array}{c}AP{P}_{E1}+{\beta}_{1}\begin{array}{c}\hfill \stackrel{{k}_{3}}{\to}\hfill \\ \hfill \underset{{k}_{-3}}{\leftarrow}\hfill \end{array}{C}_{APPE1\beta 1}\stackrel{{k}_{4}}{\to}sAPP{\beta}_{1}+C{99}_{1}+{\beta}_{1}\\ AP{P}_{E2d}+{\beta}_{2d}\begin{array}{c}\hfill \stackrel{{k}_{31}}{\to}\hfill \\ \hfill \underset{{k}_{-31}}{\leftarrow}\hfill \end{array}{C}_{APPE2d\beta 2d}\stackrel{{k}_{41}}{\to}2\cdot sAPP{\beta}_{2}+C{99}_{2d}+{\beta}_{2d}\end{array}$$

For the ODEs of the complexes above (found in Additional file 1), a quasi-steady state can be assumed. This allows the complexes to be represented as:

$$\left(\begin{array}{c}{C}_{APPCS1\alpha 1}={\alpha}_{1}\cdot AP{P}_{CS1}/{K}_{M\alpha 1}\\ {C}_{APPCS2d\alpha 2d}={\alpha}_{2d}\cdot AP{P}_{CS2d}/{K}_{M\alpha 2d}\\ {C}_{APPE1\beta 1}={\beta}_{1}\cdot AP{P}_{E1}/{K}_{M\beta 1}\\ {C}_{APPE2d\beta 2d}={\beta}_{2d}\cdot AP{P}_{E2d}/{K}_{M\beta 2d}\end{array}\right\}$$

(4)

where ${K}_{M\alpha 1}=({k}_{-5}+{k}_{6})/{k}_{5}$, ${K}_{M\alpha 2d}=({k}_{-51}+{k}_{61})/{k}_{51}$, ${K}_{M\beta 1}=({k}_{-3}+{k}_{4})/{k}_{3}$, and ${K}_{M\beta 2d}=({k}_{-31}+{k}_{41})/{k}_{31}$. These equations for the complexes are turned into the ODEs that describe the formation of the sAPP products: the ODEs of the sAPP products in the monomeric form of APP processing can be rewritten as

$$\left(\begin{array}{c}\stackrel{\cdot}{sAPP{\alpha}_{1}}=\frac{{k}_{6}\cdot {\alpha}_{1}\cdot AP{P}_{CS1}}{{K}_{M\alpha 1}}\\ \stackrel{\cdot}{sAPP{\beta}_{1}}=\frac{{k}_{4}\cdot {\beta}_{1}\cdot AP{P}_{E1}}{{K}_{M\beta 1}}\end{array}\right\}$$

(5)

whereas those in the dimeric form of APP processing can be written in the following form

$$\left(\begin{array}{c}\stackrel{\cdot}{sAPP{\alpha}_{2}}=\frac{{k}_{61}\cdot {\alpha}_{2d}\cdot AP{P}_{CS2d}}{{K}_{M\alpha 2d}}\\ \stackrel{\cdot}{sAPP{\beta}_{2}}=\frac{{k}_{41}\cdot {\beta}_{2d}\cdot AP{P}_{E2d}}{{K}_{M\beta 2d}}\end{array}\right\}$$

(6)

From Equations (5) and (6), we obtain

$$\left(\begin{array}{c}\stackrel{\cdot}{sAPP{\alpha}_{\mathit{Tot}}}=\stackrel{\cdot}{sAPP{\alpha}_{1}}+\stackrel{\cdot}{sAPP{\alpha}_{2}}\\ \stackrel{\cdot}{sAPP{\beta}_{\mathit{Tot}}}=\stackrel{\cdot}{sAPP{\beta}_{1}}+\stackrel{\cdot}{sAPP{\beta}_{2}}\end{array}\right\}$$

(7)

Fourthly, recall that *APP*_{G1} binds to *SORLA*_{G1} with a binding affinity of *K*_{S1} in monomer processing, whereas *APP*_{G2} to *SORLA*_{G2} with a different binding affinity *K*_{S2} in dimer processing. Note that
${K}_{S1}={k}_{s1}/{k}_{-s1}$ and
${K}_{S2}={k}_{s2}/{k}_{-s2}$. These properties are reflected in the biochemical reactions below:

$$\begin{array}{l}AP{P}_{G1}+SORL{A}_{G1}\begin{array}{c}\hfill \stackrel{{k}_{s1}}{\to}\hfill \\ \hfill \underset{{k}_{-s1}}{\leftarrow}\hfill \end{array}{C}_{APPG1SORLAG1}\\ AP{P}_{G2}+SORL{A}_{G2}\begin{array}{c}\hfill \stackrel{{k}_{s2}}{\to}\hfill \\ \hfill \underset{{k}_{-s2}}{\leftarrow}\hfill \end{array}{C}_{APPG2SORLAG2}\end{array}$$

We then took into consideration the rapid-equilibrium assumption for the *C*_{APPG1SORLAG1} and *C*_{APPG2SORLAG2} complexes, which gives us

$$\left(\begin{array}{c}{C}_{APPG1SORLAG1}={K}_{S1}\cdot AP{P}_{G1}\cdot SORL{A}_{G1}\\ {C}_{APPG2SORLAG2}={K}_{S2}\cdot AP{P}_{G2}\cdot SORL{A}_{G2}\end{array}\phantom{\rule{1em}{0ex}}\right\}$$

(8)

Lastly, the ODEs include conservation laws for the molecule numbers of enzymes and substrates. Herein, we took into account Equations (1) to (4) and (8) that are shown previously. For α-secretase and β-secretase, conservation law assumption leads to

$$\left(\begin{array}{l}{\alpha}_{\mathit{Tot}}={\alpha}_{\mathit{init}}+{\alpha}_{\mathit{monomer}}+{\alpha}_{\text{dimer}}\\ {\beta}_{\mathit{Tot}}={\beta}_{\mathit{init}}+{\beta}_{\mathit{monomer}}+{\beta}_{\text{dimer}}\end{array}\right\}$$

(9)

where

$$\left(\begin{array}{c}{\alpha}_{\mathit{monomer}}={\alpha}_{1}+{C}_{APPCS1\alpha 1}\\ {\alpha}_{\text{dimer}}={\alpha}_{2}+{\alpha}_{2d}+{C}_{APPCS2d\alpha 2d}\\ {\beta}_{\mathit{monomer}}={\beta}_{1}+{C}_{APPE1\beta 1}\\ {\beta}_{\text{dimer}}={\beta}_{2}+{\beta}_{2d}+{C}_{APPE2d\beta 2d}\end{array}\right\}$$

(10)

For SORLA,

$${\text{SORLA}}_{\text{Tot}}={\text{SORLA}}_{\text{monomer}}+{\text{SORLA}}_{\text{dimer}}$$

(11)

where

$$\left(\begin{array}{c}SORL{A}_{\mathit{monomer}}=SORL{A}_{G1}+{C}_{APPG1SORLAG1}\\ SORL{A}_{dimer}=SORL{A}_{G2}+{C}_{APPG2SORLAG2}\end{array}\phantom{\rule{1em}{0ex}}\right\}$$

(12)

Moreover, the ODEs also include conservation of the APP substrate. This leads to the following representation:

$$AP{P}_{\mathit{Tot}}=AP{P}_{\mathit{init}}+AP{P}_{\mathit{monomer}}+AP{P}_{\text{dimer}}$$

(13)

where

$$\left(\begin{array}{c}AP{P}_{\mathit{monomer}}(AP{P}_{\mathit{init}},{\alpha}_{\mathit{init}},{\beta}_{\mathit{init}})=AP{P}_{G1}+AP{P}_{CS1}+AP{P}_{E1}\hfill \\ +{C}_{APPCS1\alpha 1}+{C}_{APPE1\beta 1}\\ AP{P}_{\mathit{\text{dimer}}}(AP{P}_{\mathit{\text{init}}},{\alpha}_{\mathit{init}},{\beta}_{\mathit{\text{init}}})=AP{P}_{G2}+AP{P}_{G2d}+AP{P}_{CS2d}+AP{P}_{E2d}\\ +{C}_{APPCS2d\alpha 2d}+{C}_{APPE2d\beta 2d}\end{array}\phantom{\rule{1em}{0ex}}\right\}$$

(14)

in the absence of SORLA, whereas

$$\left(\begin{array}{c}AP{P}_{\mathit{monomer}}(AP{P}_{G1},{\alpha}_{1},{\beta}_{1})=AP{P}_{G1}+AP{P}_{CS1}+AP{P}_{E1}+\hfill \\ {C}_{APPCS1\alpha 1}+{C}_{APPE1\beta 1}+{C}_{APPG1SORLAG1}\\ AP{P}_{\mathit{\text{dimer}}}(AP{P}_{G2},{\alpha}_{2},{\beta}_{2})=AP{P}_{G2}+AP{P}_{G2d}+AP{P}_{CS2d}+AP{P}_{E2d}+\\ {C}_{APPCS2d\alpha 2d}+{C}_{APPE2d\beta 2d}+{C}_{APPG2SORLAG2}\end{array}\phantom{\rule{1em}{0ex}}\right\}$$

(15)

in the presence of SORLA.

Below, we discuss in more detail the properties behind equations (9)-(15): (I) each conserved equation is a function of free reactant and of reactant bound in the complexes. Note that the reactant can be the APP, SORLA, α-, or β-secretase. (II) Each *Reactant*_{Tot} function is represented differently in the presence and in the absence of SORLA. In particular, each *Reactant*_{Tot} function without SORLA is transcribed as function of *APP*_{init}, *α*_{init}*,* and *β*_{init}, whereas those with SORLA are transcribed as function of *APP*_{G1}*, APP*_{G2}*, α*_{1}*, α*_{2}*, β*_{1}*,* and *β*_{2}. (III) The amount of the *Reactant*_{init} in each *Reactant*_{Tot} with SORLA is set to be the same as that calculated from the corresponding *Reactant*_{Tot} without SORLA. (IV) Similarly, the total amount of *Reactant*_{monomer} and *Reactant*_{dimer} in each *Reactant*_{Tot} with SORLA are equivalent to those calculated from the corresponding *Reactant*_{Tot} without SORLA. (V) The amount of *SORLA*_{Tot}, *α*_{Tot}, and *β*_{Tot} are assumed to be constant for different amount of *APP*_{Tot} concentration. This assumption is based on the experimental design applied on the series of dose–response data
[14] that are used in this study. (VI) Without loss of generality, *SORLA*_{monomer} and *SORLA*_{dimer} are also assumed to be constant for different amount of *APP*_{Tot} concentration.

The properties defined above are, in particular, necessary and important. Without those properties, the presence of SORLA in monomer processing will not only affect the monomeric form of APP processing, it will also indirectly influence the dimeric form of APP processing, and vice versa. As such, it defeated the main purpose of this study, which is to differentiate the level of influence of SORLA in monomer and in dimer processing.

The development of the model described in the previous section reduced the number of free parameters from 77 to 27. The reduced number of parameter values of the model were estimated by nonlinear optimization such that the model simulations fit four biological independent dose–response series without SORLA (a total of *N*=64 experimental data points) and five biological independent dose–response series with SORLA (also a total of *N*=64 experimental data points). We looked for a set of parameter values that minimizes the weighted least squares function of *APP*_{Tot} with SORLA, and *sAṖPα*_{Tot} and *sAṖPβ*_{Tot} in regardless of SORLA (Equation (7)). On account of the different orders of magnitude of the experimental values of APP, sAPPα, and sAPPβ, weights were assigned such that the influence of each data set in the process of optimization will be equal. The weights are defined as

$$\begin{array}{c}{w}_{a}=\frac{{\sum}_{k=1}^{N}\stackrel{.}{sAPP{\alpha}_{k}^{E}}}{N},\\ {w}_{b}=\frac{{\sum}_{k=1}^{N}\stackrel{.}{sAPP{\beta}_{k}^{E}}}{N},\\ {w}_{\mathit{aS}}=\frac{{\sum}_{k=1}^{N}\stackrel{.}{sAPP{\alpha}_{S,k}^{E}}}{N},\\ {w}_{\mathit{bS}}=\frac{{\sum}_{k=1}^{N}\stackrel{.}{sAPP{\beta}_{S,k}^{E}}}{N},\\ {w}_{\mathit{appS}}=\frac{{\sum}_{k=1}^{N}AP{P}_{S,k}^{E}}{N}\end{array}$$

(16)

where the superscript ‘*E*’ and the subscript ‘*S*’ denotes experimental data points and the influence of SORLA, respectively. The goodness of fit was quantified by calculating the residual value, i.e. the sum of the squared differences between the data and model, divided by a respective weight:

$$residual=min\sum _{k=1}^{N}\left(\begin{array}{l}\frac{{(AP{P}_{S,k}^{E}-AP{P}_{Tot,S,k})}^{2}}{{w}_{\mathit{appS}}}+\frac{{(\stackrel{.}{sAPP{\alpha}_{k}^{E}}-\stackrel{.}{sAPP{\alpha}_{Tot,k}})}^{2}}{{w}_{a}}+\frac{{(\stackrel{.}{sAPP{\beta}_{k}^{E}}-\stackrel{.}{sAPP{\beta}_{Tot,k}})}^{2}}{{w}_{b}}\\ +\frac{{(\stackrel{.}{sAPP{\alpha}_{S,k}^{E}}-\stackrel{.}{sAPP{\alpha}_{Tot,S,k}})}^{2}}{{w}_{\mathit{aS}}}+\frac{{(\stackrel{.}{sAPP{\beta}_{S,k}^{E}}-\stackrel{.}{sAPP{\beta}_{Tot,S,k}})}^{2}}{{w}_{\mathit{bS}}}\end{array}\phantom{\rule{1em}{0ex}}\right)\text{.}$$

(17)

We used the lsqnonlin and fzero functions in the MATLAB optimization toolbox [22] to estimate unknown parameter values. The estimation of parameter values was performed by the steps elaborated in Additional file 1: Table S3.

We performed 500 global estimates, satisfying the condition that all parameter values are positive. Note that none of the parameter values are taken from the literature due to the differences in the experimental methods applied. Most kinetic data available in the literature on α/β-secretase activity were obtained in cell free assays with purified enzyme and artificial peptide substrate. This is in contrast to our model that relies in quantitative data obtained on APP processing in intact cell. Furthermore, the parameter values estimated for our multi-compartment model are expected to differ from that of the single-compartment model by Schmidt and colleagues [14]. Out of the 500 simulation runs, we took the set of estimated parameter values that has the smallest residual value, as shown in Additional file 1: Table S4.

Details about the protocol and assay procedures of the experimental data used in this study can be found in the paper by Schmidt and colleagues [14].

The authors declare that they have no conflict of interest.

AL initiated the ideas and concepts of the study, carried out the calculations, and performed the simulations. VS and YS participated in the design of the study, formulation of the model, analysis of the results, and editing the manuscript. TW and OW supervised the work and contributed to the writing of the final manuscript. All authors read and approved the final manuscript.

**Mathematical Modeling.** The complete formulation of the model, realized as set of ODEs. **Table S1 – Variables in the biochemical network.** The table contains the description of the variables that are used in the different compartments of the biochemical network (shown in Figure
Figure1).1). **Table S2 – Variables and parameters in the mathematical model**. The table contains the unit and the description of the variables and parameters used in the mathematical model. **Table S3 – Simulation steps.** The steps performed for the estimation of parameter values are elaborated here. **Table S4 - Estimated parameter values for Figure**Figure2.2. The table shows a set of estimated parameter values, which has the lowest residual value out of 500 simulation runs. **Figure S1 – Concentration values of the secretases with higher SORLA**_{Tot}**values.** The figure shows simulations of the influence of intermediate levels of SORLA that are greater that the value of SORLA_{Tot}, on the amount of α-secretase and β-secretase concentrations on APP processing.

Click here for file^{(1010K, pdf)}

We are indebted to Katja Rateitschak, Jana Wolf, and Thomas Millat for valuable discussion. Studies in this project were funded in part by a MSBN grant from the Helmholtz-Association (to TW and OW).

- Haass C, Selkoe DJ. Soluble protein oligomers in neurodegeneration: lessons from the Alzheimer’s amyloid beta-peptide. Nat Rev Mol Cell Biol. 2007;8:101–112. doi: 10.1038/nrm2101. [PubMed] [Cross Ref]
- Hardy J, Selkoe DJ. The amyloid hypothesis of Alzheimer’s disease: progress and problems on the road to therapeutics. Science. 2002;297:353–356. doi: 10.1126/science.1072994. [PubMed] [Cross Ref]
- Hung LW, Ciccotosto GD, Giannakis E, Tew DJ, Perez K, Masters CL, Cappai R, Wade JD, Barnham KJ. Amyloid-{beta} Peptide (A{beta}) Neurotoxicity Is Modulated by the Rate of Peptide Aggregation: A{beta} Dimers and Trimers Correlate with Neurotoxicity. J Neurosci. 2008;28:11950–11958. doi: 10.1523/JNEUROSCI.3916-08.2008. [PubMed] [Cross Ref]
- Kins S, Lauther N, Szodorai A, Beyreuther K. Subcellular trafficking of the amyloid precursor protein gene family and its pathogenic role in Alzheimer’s disease. Neurodegener Dis. 2006;3:218–226. doi: 10.1159/000095259. [PubMed] [Cross Ref]
- Thinakaran G, Koo EH. Amyloid precursor protein trafficking, processing, and function. J Biol Chem. 2008;283:29615–29619. doi: 10.1074/jbc.R800019200. [PMC free article] [PubMed] [Cross Ref]
- Jacobsen L, Madsen P, Moestrup SK, Lund AH, Tommerup N, Nykjær A, Sottrup-Jensen L, Gliemann J, Petersen CM. Molecular Characterization of a Novel Human Hybrid-type Receptor That Binds the α2-Macroglobulin Receptor-associated Protein. J Biol Chem. 1996;271:31379–31383. doi: 10.1074/jbc.271.49.31379. [PubMed] [Cross Ref]
- Yamazaki H, Bujo H, Kusunoki J, Seimiya K, Kanaki T, Morisaki N, Schneider WJ, Saito Y. Elements of Neural Adhesion Molecules and a Yeast Vacuolar Protein Sorting Receptor Are Present in a Novel Mammalian Low Density Lipoprotein Receptor Family Member. J Biol Chem. 1996;271:24761–24768. doi: 10.1074/jbc.271.40.24761. [PubMed] [Cross Ref]
- Marcusson EG, Horazdovsky BF, Cereghino JL, Gharakhanian E, Emr SD. The sorting receptor for yeast vacuolar carboxypeptidase Y is encoded by the VPS10 gene. Cell. 1994;77:579–586. doi: 10.1016/0092-8674(94)90219-4. [PubMed] [Cross Ref]
- Andersen OM, Reiche J, Schmidt V, Gotthardt M, Spoelgen R, Behlke J, von Arnim CAF, Breiderhoff T, Jansen P, Wu X, Bales KR, Cappai R, Masters CL, Gliemann J, Mufson EJ, Hyman BT, Paul SM, Nykjær A, Willnow TE. Neuronal sorting protein-related receptor sorLA/LR11 regulates processing of the amyloid precursor protein. Proc Natl Acad Sci U S A. 2005;102:13461–13466. doi: 10.1073/pnas.0503689102. [PubMed] [Cross Ref]
- Offe K, Dodson SE, Shoemaker JT, Fritz JJ, Gearing M, Levey AI, Lah JJ. The lipoprotein receptor LR11 regulates amyloid beta production and amyloid precursor protein traffic in endosomal compartments. J Neurosci. 2006;26:1596–1603. doi: 10.1523/JNEUROSCI.4946-05.2006. [PMC free article] [PubMed] [Cross Ref]
- Schmidt V, Sporbert A, Rohe M, Reimer T, Rehm A, Andersen OM, Willnow TE. SorLA/LR11 Regulates Processing of Amyloid Precursor Protein via Interaction with Adaptors GGA and PACS-1. J Biol Chem. 2007;282:32956–32964. doi: 10.1074/jbc.M705073200. [PubMed] [Cross Ref]
- Scherzer CR, Offe K, Gearing M, Rees HD, Fang G, Heilman CJ, Schaller C, Bujo H, Levey AI, Lah JJ. Loss of apolipoprotein E receptor LR11 in Alzheimer disease. Arch Neurol. 2004;61:1200–1205. doi: 10.1001/archneur.61.8.1200. [PubMed] [Cross Ref]
- Rogaeva E, Meng Y, Lee JH, Gu Y, Kawarai T, Zou F, Katayama T, Baldwin CT, Cheng R, Hasegawa H, Chen F, Shibata N, Lunetta KL, Pardossi-Piquard R, Bohm C, Wakutani Y, Cupples LA, Cuenco KT, Green RC, Pinessi L, Rainero I, Sorbi S, Bruni A, Duara R, Friedland RP, Inzelberg R, Hampe W, Bujo H, Song Y-Q, Andersen OM, Willnow TE, Graff-Radford N, Petersen RC, Dickson D, Der SD, Fraser PE, Schmitt-Ulms G, Younkin S, Mayeux R, Farrer LA, St George-Hyslop P. The neuronal sortilin-related receptor SORL1 is genetically associated with Alzheimer disease. Nat Genet. 2007;39:168–177. doi: 10.1038/ng1943. [PMC free article] [PubMed] [Cross Ref]
- Schmidt V, Baum K, Lao A, Rateitschak K, Schmitz Y, Teichmann A, Wiesner B, Petersen CM, Nykjaer A, Wolf J, Wolkenhauer O, Willnow TE. Quantitative modelling of amyloidogenic processing and its influence by SORLA in Alzheimer’s disease. EMBO J. 2012;31:187–200. [PubMed]
- Spoelgen R, von Arnim CAF, Thomas AV, Peltan ID, Koker M, Deng A, Irizarry MC, Andersen OM, Willnow TE, Hyman BT. Interaction of the Cytosolic Domains of sorLA/LR11 with the Amyloid Precursor Protein (APP) and beta-Secretase beta-Site APP-Cleaving Enzyme. J Neurosci. 2006;26:418–428. doi: 10.1523/JNEUROSCI.3882-05.2006. [PubMed] [Cross Ref]
- Andersen OM, Schmidt V, Spoelgen R, Gliemann J, Behlke J, Galatis D, McKinstry WJ, Parker MW, Masters CL, Hyman BT, Cappai R, Willnow TE. Molecular Dissection of the Interaction between Amyloid Precursor Protein and Its Neuronal Trafficking Receptor SorLA/LR11† Biochemistry. 2006;45:2618–2628. doi: 10.1021/bi052120v. [PubMed] [Cross Ref]
- Dodson SE, Andersen OM, Karmali V, Fritz JJ, Cheng D, Peng J, Levey AI, Willnow TE, Lah JJ. Loss of LR11/SORLA enhances early pathology in a mouse model of amyloidosis: evidence for a proximal role in Alzheimer’s disease. J Neurosci. 2008;28:12877–12886. doi: 10.1523/JNEUROSCI.4582-08.2008. [PMC free article] [PubMed] [Cross Ref]
- Kitano H, Funahashi A, Matsuoka Y, Oda K. Using process diagrams for the graphical representation of biological networks. Nat Biotech. 2005;23:961–966. doi: 10.1038/nbt1111. [PubMed] [Cross Ref]
- Cell Designer. The Systems Biology Institute, Tokyo, Japan; 2008.
- Willnow TE, Petersen CM, Nykjaer A. VPS10P-domain receptors - regulators of neuronal viability and function. Nat Rev Neurosci. 2008;9:899–909. doi: 10.1038/nrg2454. [PubMed] [Cross Ref]
- Andersen OM, Willnow TE. Lipoprotein receptors in Alzheimer’s disease. Trends Neurosci. 2006;29:687–694. doi: 10.1016/j.tins.2006.09.002. [PubMed] [Cross Ref]
- MATLAB. Natick, Massachusetts: Math Works Inc; 2009.

Articles from BMC Systems Biology are provided here courtesy of **BioMed Central**

PubMed Central Canada is a service of the Canadian Institutes of Health Research (CIHR) working in partnership with the National Research Council's national science library in cooperation with the National Center for Biotechnology Information at the U.S. National Library of Medicine(NCBI/NLM). It includes content provided to the PubMed Central International archive by participating publishers. |