Nonlocal models in the analysis of brain neurodegenerative protein dynamics with application to Alzheimer's disease.

It is well known that today nearly one in six of the world's population has to deal with neurodegenerative disorders. While a number of medical devices have been developed for the detection, prevention, and treatments of such disorders, some fundamentals of the progression of associated diseases are in urgent need of further clarification. In this paper, we focus on Alzheimer's disease, where it is believed that the concentration changes in amyloid-beta and tau proteins play a central role in its onset and development. A multiscale model is proposed to analyze the propagation of these concentrations in the brain connectome. In particular, we consider a modified heterodimer model for the protein-protein interactions. Higher toxic concentrations of amyloid-beta and tau proteins destroy the brain cell. We have studied these propagations for the primary and secondary and their mixed tauopathies. We model the damage of a brain cell by the nonlocal contributions of these toxic loads present in the brain cells. With the help of rigorous analysis, we check the stability behaviour of the stationary points corresponding to the homogeneous system. After integrating the brain connectome data into the developed model, we see that the spreading patterns of the toxic concentrations for the whole brain are the same, but their concentrations are different in different regions. Also, the time to propagate the damage in each region of the brain connectome is different.

Introduction

Alzheimer’s disease (AD) is the most prevalent neurodegenerative disorder that slowly increases degrees of dementia over time[1]. In most people with AD, the first symptoms appear after 60 years of age, whereas some persons reported AD before 60 years of age, and it is not frequent. According to the Alzheimer’s Association, five years ago 5.5 million Americans in the U.S. and 44 million worldwide dealt with this disease in all ranges of ages[2]. At the same time, of the population over 65 years of age is suffering, and it is predicted to be a public health crisis in the coming decades caused by this disease. Baumgart et al. estimated that to increase to 131 million people by 2050[3]. The history of the disease goes back well over 100 years ago, when Dr. Alois Alzheimer in 1906 noticed that a woman died with an unusual mental illness, and the symptoms were memory loss, language problems, etc. He examined her dead brain and identified many amyloid plaques and neurofibrillary tangles (NFTs) in the brain tissue[4-6]. Today a wealth of medical and bioengineering tools exist for the detection, prevention, and treatments of such disorders. Researchers believe that the toxic changes in the brain destroy the communications between neurons[2,7].

The research involving amyloid-beta (A) in AD has progressed more quickly than tau protein (P). Amyloid-beta is a small piece of a larger protein amyloid precursor protein (APP)[2]. It is still largely unknown how amyloid-beta accumulates in the brain and initiates AD. According to “the amyloid hypothesis”, it is believed that flaws in the processes governing the production, accumulation, or disposal of amyloid-beta are the primary cause of AD[2,8]. Recently, the Food and Drug Administration (FDA) approved the first drug “Aducanumab” for the treatment of some cases of Alzheimer’s disease[9]. This drug is a human antibody that targets the amyloid-beta protein and helps to reduce amyloid plaques associated with AD. It is still in the trial stage and there is no royal road to overcome AD[10,11]. Researchers continue to study whether this medication affects a person’s rate of cognitive decline over time. In the AD-affected brain, A develops slowly and has many stages of building. Initially, A forms small clusters (oligomers), then a chain of oligomers (fibrils), and then “mats” of fibrils called beta-sheets. In the final stage, beta-sheets clump with other substances, form plaques. These amyloid plaques block the blood vessels and disrupt cell-to-cell communications (blood circulations) and active immune cells. These immune cells trigger inflammation and destroy the brain cells. Researchers believe that, along with the A, there may be some other factor in the AD onset and progression.

There is a lot of evidence of the spread of tau protein with AD progressing[12-14]. Along with the accumulation of A in the brain, scientists have also focussed on the importance of tau protein in AD. The tau proteins are neuronal microtubule (MT)-associated proteins found in the axon. Under physiological conditions, tau acts like a highly soluble and unfolded protein. It interacts with tubulin and promotes its assembly into MTs, which helps to stabilize their structure[12,15]. Also, the tau protein plays an essential role in the balance between MT-dependent axonal transport of organelles and biomolecules[16,17]. At the first stage of AD progression, tau protein continuously accumulates in the entorhinal cortex and hippocampus regions. Then tau protein spreads into extracellular space in the form of oligomers[18]. Oligomers travel through a structure known as the synapse. It allows passing electrical or chemical signals from one brain cell to another through diffusion[13]. This abnormal tau distribution over the brain causes further disease progression.

Many cellular mechanisms and their interactions between two different groups of proteins are not fully understood. In particular, the interactions between A and are also, at a great degree unknown. Nowadays, researchers have focussed on finding the possible protein–protein interactions between two or more groups of proteins in AD[14,19,20]. Following Walker et al., for AD, “the amyloid-beta-tau nexus is central to disease-specific diagnosis, prevention and treatment”[21]. Thompson et al. observed some crucial interactions of these proteins[22]: (i) amyloid-beta enhances the new seeding of toxic tau protein, (ii) toxic concentration of amyloid-beta depends on the concentration of tau protein, and (iii) amyloid-beta and tau protein enhance each other’s toxicity.

A mathematical model is an excellent computational tool to work with these interactions. Many researchers considered these protein–protein interactions in mathematical models related to AD[23-28]. However, some questions regarding AD remain open, and more refined models need to be developed to address them. Experimental studies suggest non-uniform distributions of the proteins in the brain connectome[29]. In this work, our primary goal is to capture such types of non-uniform distributions of the proteins in the brain through a mathematical model. In the mathematical model, we apply a multiscale model that describes the interaction of two different protein families: A and P. Each of the proteins follows the dynamics of a modified heterodimer model with local and nonlocal interactions and a coupling parameter between two proteins[22,30-33]. Motivated by Meisl et al.[34], we replace exponential growth with logistic growth in both the proteins A and P. Due to toxic variants of these proteins present in the brain cell, neuronal damage occurs in the brain. We consider a model with local and nonlocal interactions for capturing such type of neuronal damage in the brain cells[22].

The rest of this manuscript is organized as follows. We consider a continuous nonlocal model for protein–protein interactions in Sect. 2. A nonlocal version of the damage of the cell is also described here. We have derived a nonlocal network model that can be used for integrating the brain connectome data. In Sect. 3, we discuss the temporal behaviour of the stationary points of the continuous mathematical model and provide details of the associated network model. The propagation of the toxic concentrations in the brain connectome, related to primary, secondary, and mixed tauopathies, is discussed in Sect. 4. Finally, the summary of the work is given in Sect. 5.

Mathematical model

Continuous model

We consider a spatial domain in . For and time , and represent the concentrations of healthy and toxic A respectively. Similarly, and represent the concentrations of healthy and toxic P respectively. Following the ideas of[22,30,32,33], the concentration evolution is given by the set of coupled integro-differential equations: where the first two equations correspond to the healthy and toxic variants of the protein A and the last two equations play the same role for P. Model (1) is an example of nonlocal models which become increasingly important in diverse areas of applications[35]. Here, and are the mean production rates of healthy proteins, and are the mean clearance rates of healthy and toxic proteins, and and represent the mean conversion rates of healthy proteins to toxic proteins. The parameter is the coupling between the two proteins A and P, and it is considered due to the fact that the A enhances the seeding of new toxic P concentration[22]. The parameters and have units of the reciprocal concentrations of the healthy proteins of A and P, respectively. The diffusion tensors , , and characterize the spreading of each proteins. For any time instance, the convolution term at the spatial point is given byIn a similar way, we define the other convolutions in (1). Here, is the kernel function and it describes the conversion efficiencies between the spatial points and . We assume that the kernel function is non-negative, even function, and has a compact support in . Also, satisfies a standard normalized condition:Based on our application area, it is reasonable to assume that all variables and initial conditions are non-negative, and also that all the parameters are strictly positive.

Following Meisl et al.[34], the fundamental model that includes the replication takes into account the logistic growth rather than exponential growth in the reaction-diffusion model. Therefore, we modify the exponential growth of healthy A and healthy P of the model (1) by the logistic growth that leads to the following system: with appropriate initial and boundary conditions for all the components.

The above system (2) dictates the spread of two healthy and toxic variants of the proteins throughout the domain . The increase in the density of toxic proteins at a spatial point in disrupts the extracellular environment near and causes intracellular activities. The impact of these toxic variants is not fully understood. However, Thompson et al. observed their correlation and defined a gross measure of regional neuronal damage by a function at a spatial point at time t[22]. We say the neuron is healthy (functional) at a spatial point if . On the other hand, if , then the neuron is no longer active, i.e., the neuron reached a fully damaged asymptotic state. The evolution of the damage is described aswith the initial condition , where and account for the damaging effect of toxic A and P, respectively. The coefficient represents the damage due to the combined presence of both toxic loads. Finally, is the rate of transneuronal damage propagation, and it reflects aggregate neuronal damage from regional neighbours. Here, all the ’s are non-negative, and denotes the convolution, defined as before.

Network model

Having defined the nonlocal continuous model, we are now in a position to develop a nonlocal network mathematical model that would account for the brain network data. The development of a coarse-grain mathematical model based on a continuous system is a key to defining a new network model for the brain connectome[36,37]. Let be a weighted brain graph with N number of nodes and E be the number of edges defined in the domain . In the model, each node corresponds to a small area (1–1.5 cm) of the gray matter, called Regions of Interest, and the edges of the network are the axonal bundles in white matter tracts[38,39]. Suppose, is an adjacency matrix corresponding to the weighted graph . With the help of the adjacency matrix , we construct the Laplacian for the graph. For , we denote the elements of aswhere is the mean fiber number and is the mean length squared between the nodes i and j. Here, the mean fiber length is measured in millimetres[40]. The Laplacian of the graph is a square matrix of order N with the elementswhere is the diffusion coefficient and are the elements of the diagonal weighted degree matrix.

Now, we begin to define the convolution term in the network model similar to the continuous model (2) for each of the nodes . In the continuous integro-differential model, the convolution term at a spatial point is an integral in the neighbourhood (inside the domain of definition) of the spatial point with Gaussian weight function[33]. Again, for the continuous model, all the neighbourhood points of are connected. We apply the same connectedness technique for finding the convolution term at each node in the network model.

First, following[33], we fix a node j and define a set of nodes containing all directly connected nodes to the node j. Then, we find another set of nodes containing all the immediate connected nodes, connected with the nodes in the set . We continue this process in the whole graph until we find the complete set of connected nodes (directly or indirectly). In this procedure, we obtain “complete sets” of nodes connected with the node j and let be the union of these sets of nodes. Now, we sort all the nodes in according to the “shortest distance” from the node j. Suppose, is the total number of nodes connected with the node j including the self node j and are the labels of the sorted connected nodes starting with the self node j, i.e., . We define a set of elements of weights aswhere denotes the shortest distance from the node to the node j along the edges.

We normalize as , whereNow, we define a row vector of N number of elements with non-zero elements being at the indices respectively. Finally, we are ready to define the convolution at the node j in the graph asIf the node j is not connected with any other nodes, then . Similarly, we define the other convolutions.

Taking all the factors (Laplacian and convolutions), we build a nonlocal network mathematical model on the brain graph . Let, be the concentrations of healthy and toxic A and be the concentrations of healthy and toxic P at the node j. Then for , the network equations corresponding to the continuous model (2) is a system of first order differential equations and it is given by with non-negative initial conditions.

In a similar fashion, we define the network model for the nonlocal version of the damage equation (3). Let be the neuronal damage of cell located at the node j. Then for all the nodes , the network damage equation corresponding to the network model is given bywith the initial condition . Similar to the nonlocal interaction defined in (4), we use the same technique to find the convolution for the neuronal damage corresponding to the node j, but with a different controlling parameter in the place of .

Analysis of the continuous model

Before analyzing the full reaction-diffusion model consisting of a coupled set of differential equations (2), we first study the temporal dynamics (e.g., equilibrium points and their stabilities) of the model. For studying the temporal dynamics of a reaction-diffusion model, we ignore all the spatial dependencies in the model. In this case, the diffusion term becomes zero and the convolution term , for any w. Therefore, under these simplifications, the temporal model corresponding to the nonlocal model (2) is given by with non-negative initial conditions. The system (7) has some feature points, such that characterize its behaviour but do not change over time, called equilibrium points. Note that these temporal dynamics apply to both continuous and network models. The stationary states or the solutions of the continuous and network models are the equilibrium points of the system (7).

The system (7) always has a trivial equilibrium point . This equilibrium point corresponds to the the absence of all the ingredients in the system. Generally, it does not occur in the living brain but occurs only in the dead brain cell. Depending on the parameter values, we have a maximum of seven equilibrium points of the system (7) where at least one component of the equilibrium point is zero. These are called the axial equilibrium points. The possible axial equilibrium points areThe concentration in each of the components cannot be negative. Therefore, for the existence of the equilibrium points and , we must have and . Similarly, for the feasibility of the equilibrium points and , we need and .

We call a stationary state healthy A (healthy P) if the second (fourth) component of the equilibrium point is zero; otherwise, it is toxic A (toxic P). The equilibrium points and contain either healthy A or healthy P concentrations, and these are healthy A and healthy P stationary states, respectively. is a “healthy A–healthy P” stationary state since it does not have any of the toxic loads A or P, i.e., no amyloid plaques or neurofibrillary tau tangles. Also, we call this equilibrium point the disease-free stationary state. The other equilibrium points are brain disease states since these brain states have either amyloid plaques or neurofibrillary tau tangles. Now, we are interested in finding an equilibrium point where both toxic loads are present.

Suppose, is a positive equilibrium point. Then we have , , , where satisfies the quadratic equationFor the feasibility of the first two components of the equilibrium point , we must have and . Along with these conditions, we also require , where is a positive root of the Eq. (8). In this case, the equilibrium point becomes “toxic A–toxic P” stationary state as it contains both amyloid plaques and neurofibrillary tau tangles.

Now, we briefly discuss the stability of the equilibrium points of the system (7). The stability behaviour of an equilibrium point of the system depends only on the eigenvalues of the Jacobian matrix calculated at that point. For a fixed equilibrium point, if the real parts of all the eigenvalues of the Jacobian matrix are negative, the equilibrium point is stable; otherwise, it is unstable. The Jacobian matrix of the system (7) about any equilibrium point is given bywhere , and .

Therefore, the eigenvalues of the Jacobian matrix are , where and . Depending on the parameter values, the number of stationary points is different. Hence, in the next section, we discuss the stability of the stationary points after fixing the parameter values, motivated by[22].

Results and discussion

Tau is a microtubule-associated protein predominantly expressed in nerve cells promoting microtubule assembly and stabilization. It is a cytosolic protein mainly present in axons and involved in anterograde axonal transport. Tau protein alters its metabolism in several neurodegenerative diseases (e.g., AD). Thus alterations in the amount of the tau protein, missense mutations, post-transcriptional modifications, aberrant tau aggregation, or a different expression of some of its isoforms could provoke pathological effects resulting in the appearance of neuronal disorders known as tauopathies. Researchers are making a substantial effort to generate tau oligomers by purified recombinant protein strategies for studying oligomer conformations and their toxicity. Still, there is no specific toxic tau protein has been identified[41]. However, some cellular biosensor technology has been discovered to monitor the formation of tau oligomers and aggregates in the living cells[41], e.g., fluorescence resonance energy transfer (FRET), bimolecular fluorescence complementation (BiFC), and split luciferase complementation (SLC).

Based on the pathology, there are two groups of disease propagation which are classified as: primary, when tau is the main lesion, and secondary, when tau is associated with other pathology[42,43]. Here, for the primary tauopathy, the non-zero concentration of toxic P exists independent of the concentration of A. On the other hand, this dependency needs for the secondary tauopathy. After observing several characteristics of the models, we move to an interesting scenario in the brain network model where some regions in the brain satisfy the primary tauopathy, and the other regions satisfy the secondary tauopathy. This is called a mixed tauopathy, and it is more realistic than primary or secondary[44]. We also study the damage dynamics in each of the tauopathies.

In this work, we use actual brain connectome data[38-40] for the network model available on the website https://braingraph.org. We choose high-resolution data consisting of 1,015 nodes and a different number of edges. We compared different data sets in the brain connectome to check the numerical artifact, and the results are consistent. We have used the C programming language and Matlab for all the computations and simulations.

As the first step in the disease propagation analysis, we set seeding sites of toxic amyloid-beta and toxic tau proteins in the brain region(s) associated with Alzheimer’s disease. In the AD-affected brain, researchers observed that the amyloid-beta plaques start their propagation from temporobasal and frontomedial regions in the brain[19,22,27,45,46]. All the nodes in these two regions are the initial seeding sites for the toxic A [see Fig. 1]. On the other hand, locus coeruleus and transentorhinal associated regions are the initial tau staging sites for the AD[34,47] [see Fig. 1]. Seeding sites are highlighted in red, whereas the gray color corresponds to the zero toxic loads in the brain connectome. Initially, we set all the nodes to a healthy state (i.e., ) except the seeding nodes. We add small toxic loads and in the toxic A and toxic P concentrations for the seeding sites, respectively. In modelling other scenarios, one may choose some other percentages of the toxic loads for the computations. We use the same initial conditions every time for finding the solution of the network model.

Figure 1

(Color online) Initial seeding sites for Alzheimer’s disease: (left) toxic A and (right) toxic P.

Primary and secondary tauopathies

For the analysis of disease progression in the brain initially, we take a globally-constant synthetic parameter values throughout all the regions in the brain connectome. The first example is corresponding to the primary tauopathy. All the uniform parameter values for the primary tauopathy are in Table 1 except the coupling parameter and the parameters and . The coupling parameter is not a passive facet of disease phenomenology but plays a much more integral role in secondary tauopathy[22]. In the primary tauopathy, first we fix and study the effect of the parameters and on the local network model.

Table 1

Fixed parameter values for the primary tauopathy.

Healthy Aβ	Toxic Aβ	Healthy τP	Toxic τP
ρ=1	ρ=1	ρ=1	ρ=1
a_0=1.035	a~_1=0.828	b_0=0.69	b~_1=0.552
a_1=1.38	a_2=1.38	b_1=1.38	b_2=1.38

Firstly, we take . All possible stationary states exist (described in the earlier section) for the system with these parameter values. Here, all the stationary states are unstable except for the positive stationary state . Therefore, the disease propagates from the initial seeding sites and converges to the positive steady state of the local system. We plot the spatial averages of the solutions of the local network model in Fig. 2a. Similarly, the spatial averages of the solutions of the local network model for two other values of the parameters and are shown in Fig. 2b,c. With an increase in the parameters and , the required time to propagate the disease in the whole brain connectome increases.

Figure 2

(Color online) Spatial averages of the solutions of the local model with different values of and : (a) , (b) and (c) . The curves (solid), (dashed-dashed), (dashed-dotted) and (dotted) are corresponding to the spatial averages of u, , v and , respectively.

We see some spikes in the spatial-average solution of toxic P in Fig. 2. Here, two eigenvalues of the Jacobian matrix around the positive stationary states are negative real numbers, but the other two are complex numbers with negative real parts. Therefore, the solution converges to the stable positive stationary state, but oscillations happen in a small neighbourhood of the point before converging to the solution. With an increase in the parameters and , the amplitude of the oscillations decrease [see Fig. 2b,c]. Generally, at the initial stage, toxic P grows at the seeding sites and distributes the concentration throughout the brain connectome. After this spreading stage, toxic P concentration rises in the whole brain. So the average density of toxic P grows very fast, and it crosses the stationary state component corresponding to toxic P. Generally, the cerebrovascular system reduces the toxic density in the brain and keeps it free from harmful agents. Due to the high accumulation rate of toxic P all over the brain, the cerebrovascular system can not control the whole toxic density. At last, it reduces some toxic levels of P, not the total concentrations, and saturates to the disease state “toxic A–toxic P”.

Now, we analyze the general dynamics at each node located in different regions in the brain connectome. We provide the parameter values in Table 1 along with and . The distributions of toxic A and toxic P in the brain connectome are shown in Fig. 3. Here, non-uniform density distributions of all four components exist in the brain connectome, and they occur due to the initial seeding of toxic loads in some particular regions. To verify the non-uniformity, we plot the solution profiles for toxic A and toxic P corresponding to each of the nodes in Fig. 4. Here, we see three different solution behaviours for toxic A and toxic P. For the case of toxic A, the dashed-dotted, solid, and dotted curves correspond to the seeding connected nodes (temporobasal and frontomedial regions), disconnected nodes, and the rest of the nodes, respectively. The same types of curves are plotted for toxic P.

Figure 3

(Color online) Toxic A and toxic P concentrations in the brain connectome at different times (left to right: and 100) for the local model corresponding to the primary tauopathy. Upper and lower panels are corresponding to toxic A and toxic P, respectively.

Figure 4

(Color online) Node-wise toxic propagations over time: (a) toxic A and (a) toxic P.

The solutions of toxic A and toxic P converge to three different concentrations as demonstrated by Fig. 4. But, the system has only one stable stationary point. This means the solution converges to some other stationary point(s) . If a node is connected with other nodes, then the solutions converge to the positive stationary point . On the other hand, suppose a node is not connected to the other nodes. Now, if this node is inside the seeding sites of toxic A or toxic P, then the solution converges to or , respectively; otherwise, it remains in “healthy A–healthy P” state, i.e., converges to . This happens due to the existence of stable manifolds of the stationary points , , and , called semi-stable stationary points. These results are distinct from those obtained earlier[22,32,33]. Indeed, if we consider the exponential growth in A and P, i.e., the network model corresponding to the system (1), then we can not get semi-stable stationary points[22].

Now, we are in a position to analyze the neuronal damage of the nodes in the brain connectome in further detail. We solve the local and nonlocal versions of the neuronal damage model (6) along with the local interaction of the model (5). In this case, we take the temporal parameter values from Table 1 along with , . We choose less than to see the influence of toxic P neurofibrillary tangles on neural damage[48-50] and larger than , for enhancing the toxic effect of P in the presence of toxic A[50-54]. Following these parameter restrictions, we fix , and . For , we plot the neuronal damage solutions for the local and nonlocal models () in Fig. 5. Similar to Fig. 4, four types of damage profile exist for the neuronal damage model. Solid, dashed-dotted, dashed-dashed and dotted curves are corresponding to the neuronal damages of the nodes which are converging to the stationary points , , and , respectively. The time taking to damage the nodes for the nonlocal model is larger compared to that of the local model. Also, an increase in the parameter causes faster neuronal damage in the brain connectome [see Fig. 6].

Figure 5

(Color online) Node-wise damage dynamics over time for local and nonlocal interactions in the neuronal damage model: (a) local and (b) nonlocal with .

Figure 6

(Color online) Spatial average of the damage of all the nodes with respect to time for different values of : (dotted), (dashed-dotted) and (solid).

For the case of secondary tauopathy, the stationary points and do not exist. Therefore, if a node is connected with other nodes or it is located in the seeding sites of toxic P, solutions converge to the positive stationary point . Now, if a node is disconnected and it is located inside the seeding sites of toxic A, then the solution converges to ; otherwise, it remains in a healthy stationary state. In this case and are semi-stable stationary points.

The toxic propagations and the concentrations for the local and nonlocal network models are the same, corresponding to both the tauopathies. This happens due to the uniform parameter values all over the regions in the brain connectome. However, a qualitative change can be seen in the nonlocal interactions for non-uniform parameter values in the brain connectome. We study this behaviour in the following subsection.

Mixed tauopathy

We have studied the local and nonlocal network models with global constant parameter values. In the analysis of the whole living brain, this consideration of uniform parameters would need to be generalized further. The positron emission tomography (PET) imaging studies of A and P radiotracer uptake provide us with a better insight[55]. In AD, the distribution of PET- SUVR intensities is biased in different parts of the brain. We incorporate this idea in the network model along with the AD patient data. We consider a sample parameter data provided in the Alzheimer’s Disease Neuroimaging Initiative (ADNI) database (adni.loni.usc.edu)[22].

For the mixed tauopathy, we study a mixed version of primary and secondary tauopathies, i.e., some parts in the brain connectome satisfy secondary tauopathy, and the rest satisfies primary tauopathy. A mixed neuropathology with amyloid-beta plaques and tau pathology is present in nearly of patients and are associated with an increased rate of dementia and a decreased survival time[56]. First, we set all the nodes in the brain connectome to a secondary tauopathy by considering the parameters in Table 2, and . According to the ADNI data, a modified value of the coupling parameter in some of the regions, given in Table 3, and the rest regions . The modification is symmetric, i.e., has the same value for both the left and right hemispheres in the corresponding regions. Motivated by[22], we also modify and in some of the regions [see Table 4], so that the system has a state of primary tauopathy.

Table 2

General synthetic parameter values.

Healthy Aβ	Toxic Aβ	Healthy τP	Toxic τP
ρ=1.38	ρ=0.138	ρ=1.38	ρ=0.014
a_0=1.035	a~_1=0.828	b_0=0.69	b~_1=0.552
a_1=1.38	a_2=1.38	b_1=1.38	b_2=1.035

Table 3

Modified values.

Brain region ID and modified b_3 value
Pars Opercularis	7.452	Rostral middle frontal gyrus	6.707
Superior frontal gyrus	7.452	Caudal middle frontal gyrus	7.452
Precentral gyrus	5.589	Postcentral gyrus	3.726
Lateral orbitofrontal cortex	6.486	Medial orbitofrontal cortex	6.486
Pars triangularis	5.520e-6	Rostral anterior cingulate	6.210e-6
Posterior cingulate cortex	3.45	Inferior temporal cortex	13.11
Middle temporal gyrus	11.04	Superior temporal sulcus	8.97
Superior temporal gyrus	8.28	Superior parietal lobule	12.42
Cuneus	13.8	Pericalcarine cortex	13.8
Inferior parietal lobule	11.73	Lateral occipital sulcus	15.18
Lingual gyrus	13.8	Fusiform gyrus	7.59
Parahippocampal gyrus	11.04	Temporal pole	1.104e-5

Table 4

Modified and values.

Brain region	Entorhinal cortex	Pallidum	Locus coeruleus	Putamen	Precuneus
b_2	3.125	2.76	1.38	3.795	3.105
b_3	1.104e−5	2.76	1.38	3.795	3.105

We use the same seeding sites as the initial condition for the mixed tauopathy. The solutions of the local and nonlocal models corresponding to the toxic P are shown in Fig. 7. Due to the mixed parameter values in the brain region, the concentrations (all four components) are different in each node in the brain connectome. With decreasing the parameter values of , we see a higher toxic density accumulation in some nodes. The targeted nodes depend on the degree of the node, i.e., the number of nodes connected with the node. Therefore, the spreading pattern of the toxic concentrations for the mixed tauopathy is different from primary or secondary tauopathies.

Figure 7

(Color online) Toxic-P distributions in the brain over time (left to right: and 100) for the mixed tauopathy. Upper panel is corresponding to the local model, while middle and lower panels are corresponding to the nonlocal model with and , respectively.

We plot the dynamics of toxic P in each of the regions in Fig. 8. The spreading profiles of the toxic concentration in different nodes are different. Therefore, the time required to damage the nodes is different [see Fig. 9].

Figure 8

(Color online) Node-wise toxic-P distributions over time for the nonlocal model () with mixed tauopathy: (a) different curves corresponding to different brain regions; (b) zoom version of (a).

Figure 9

(Color online) Region-wise damage for the nonlocal model with mixed tauopathy.

Summary

A modified nonlocal coupled heterodimer multiscale model has been developed and applied for a better understanding of the dynamics of Alzheimer’s disease. We have derived a nonlocal network model corresponding to the model based on a coupled system of integro-differential equations and incorporated the brain connectome data into this new model. Further, we have considered a nonlocal interaction in the damage dynamics. The stability behaviour of each of the stationary points of the system has been analyzed by using linear stability analysis. Two types of tauopathies (primary and secondary) have been discussed in detail, and their combination has been examined.

For the primary tauopathy, the spreading patterns of the toxic concentration of A and P in the brain connectome are different in different regions. This occurs due to the initial seeding sites in the brain connectome. In this case, a total of four different spreading patterns are possible. On the other hand, three different types of spreading patterns exist for the secondary tauopathy. Our analysis has also revealed that many types of spreading patterns exist in the case of mixed tauopathy. This non-uniform toxic tau protein distribution in the brain connectome confirms good agreement with experimental results[29]. The damage dynamics follow the same number of spreading patterns. Irrespective of the tauopathies, the nonlocal model takes longer to propagate the disease than the local model.

Given that tau is an MT-building protein, its more refined description can be achieved by including corresponding microtubules into consideration. Earlier developed models for the analysis of the microtubules[57-59] can prove to be useful for further refinement of the models presented here. Further, it is now known that AD may affect not only the brain but the entire body, including the cardiovascular system[60]. Amyloid-beta protein fragments that form plaques in the brain of AD patients might stiffen their heart muscles due to the deposition of sticky amyloid-beta protein. Researchers discovered that Alzheimer’s is caused by amyloid-beta proteins building up in the spaces between brain cells. While this causes noticeable symptoms in the brain, this same protein plaque can build up around the heart. The advancement of nonlocal models in this direction represents an important avenue of future research, where the methodology developed here may prove to be advantageous, along with the integration of such models with coupled modelling approaches accounting for the role of astrocytes[61,62].