Effect of mitochondrial circulation on mitochondrial age density distribution

Recent publications report that although the mitochondria population in an axon can be quickly replaced by a combination of retrograde and anterograde axonal transport (often within less than 24 hours), the axon contains much older mitochondria. This suggests that not all mitochondria that reach the soma are degraded and that some are recirculating back into the axon. To explain this, we developed a model that simulates mitochondria distribution when a portion of mitochondria that return to the soma are redirected back to the axon rather than being destroyed in somatic lysosomes. Utilizing the developed model, we studied how the percentage of returning mitochondria affects the mean age and age density distributions of mitochondria at different distances from the soma. We also investigated whether turning off the mitochondrial anchoring switch can reduce the mean age of mitochondria. For this purpose, we studied the effect of reducing the value of a parameter that characterizes the probability of mitochondria transition to the stationary (anchored) state. The reduction in mitochondria mean age observed when the anchoring probability is reduced suggests that some injured neurons may be saved if the percentage of stationary mitochondria is decreased. The replacement of possibly damaged stationary mitochondria with newly synthesized ones may restore the energy supply in an injured axon. We also performed a sensitivity study of the mean age of stationary mitochondria to the parameter that determines what portion of mitochondria re‐enter the axon and the parameter that determines the probability of mitochondria transition to the stationary state. The sensitivity of the mean age of stationary mitochondria to the mitochondria stopping probability increases linearly with the number of compartments in the axon. High stopping probability in long axons can significantly increase mitochondrial age.


| INTRODUCTION
The primary purpose of mitochondria in the cells is to generate easily accessible chemical energy in the form of ATP.In addition, mitochondria serve various other functions, such as calcium buffering.They are also known for their involvement in apoptosis, [1][2][3] particularly in the context of stroke-induced injury. 4ince mitochondrial proteins have a limited half-life, these proteins (or the entire mitochondrion) need to be periodically replaced. 5Old and damaged mitochondria are usually degraded in somatic lysosomes. 6,7Mitochondria are continuously transported in axons. 8The representative length of mitochondria is between 0.5 and 10 μm. 9 Hence, mitochondria are too large to exhibit any significant diffusivity.Mitochondria transport is accomplished by molecular motors, kinesin-1 (and possibly kinesin-3) in the anterograde direction and by cytoplasmic dynein in the retrograde direction. 10,11Mitochondria move with an average velocity of 0.4-0.8μm/s. 12][15] In neurons, most of the degradation of old and damaged mitochondria occurs in the soma.7][18] However, Reference 19 reported that after retrogradely transported mitochondria depart from the axon terminal, a portion of them remain in the axon for days.The persistence of mitochondria after leaving the terminal for a prolonged duration may indicate that they will never reach the soma, as they spend much time in the stationary state.Alternatively, this may suggest that retrogradely transported mitochondria undergo recycling in the soma.Based on this evidence, we hypothesize that not all mitochondria that enter the soma from the axon undergoing retrograde movement are degraded in the soma; a portion return to the axon.This hypothesis is supported by unpublished observations reported in Reference 3.These observations suggest that retrogradely moving mitochondria, which are traditionally expected to be segregated for degradation through mitophagy upon entering the soma, instead merge with the existing mitochondria in the soma.This fusion leads to the exchanging of components between the returning mitochondria and the newly synthesized mitochondria in the soma, producing repaired mitochondria that are ready to return to the axon.
A high level of energy is required for neurons to survive an injury.Since the distance to which ATP can be transported by diffusion is limited, replacing damaged mitochondria with new ones is critical for neuron survival. 20Axonal regeneration after injury or ischemic stress 21 can be promoted by turning off mitochondria anchoring.This leads to faster replacement of damaged mitochondria with healthy mitochondria. 22Our model simulates the reduction of mitochondrial anchoring by decreasing a value of the parameter that characterizes the probability of mitochondrial anchoring, p s .
The return of mitochondria to the axon must be associated with changing the type of molecular motors that move the mitochondria.Since mitochondria are driven toward the soma by dynein motors, to turn around and re-enter the axon they must switch from dynein motors to kinesin motors.
Reference 23 used mean-field equations to study the delivery of motor-cargo complexes to branched axons.A quantity called mitochondrial health was introduced in Reference 24 and treated as a proxy for mitochondrial membrane potential.Models that treat mitochondrial health as a decaying component were developed in Reference 25.This paper adopts a different approach by using a compartment-based model.Compartments do not necessarily need to be physical volumes with actual bottlenecks between them.A compartmental model, for example, could be used to simulate a tracer in tissue, which binds to various specific target receptor sites. 26,27In this paper we need to simulate mitochondria transport between various mitochondrial demand sites.Each demand site is represented by three compartments containing populations of anterograde, stationary, and retrograde mitochondria.We developed this model in our previous publications and utilized it for determining the mean age of mitochondria along the axon length 28 and for computing the mean age and age density distributions of mitochondria in a branching axon. 29he main question that we ask in this paper is how the mean age and age density distributions of mitochondria in the compartments located at different distances from the soma are affected by mitochondria returning to the circulation.According to our hypothesis, some of the mitochondria that exit the axon change the motors that propel them to kinesin and return to the circulation.This would explain the presence of older mitochondria in the axon, reported, for example, in Reference 19.Intuitively, this should increase the mean age of mitochondria in the axon as a result of the described process in which older mitochondria re-enter the axon.Another question we attempt to answer in this paper is whether making mitochondria less likely to stop would contribute to a faster renewal of mitochondria (which in our model will appear as a reduction of the mean age of mitochondria).This would improve chances for neuronal recovery after ischemic stress or injury. 22

| Problem setup
Mitochondria usually dock near energy demand sites in the axon. 30We simulated the axon as consisting of N demand sites (Figure 1).Mitochondria in the axon can be divided into stationary, anterograde, and retrograde, with transitions possible between these three pools. 13We therefore represented each demand site as that consisting of three compartments which contain stationary, anterograde, and retrograde mitochondria.The total number of compartments in our multi-compartment model [31][32][33] is thus 3 N, see Figure 2 (we used N = 10).Following Reference 25, we define the demand sites as narrow sections of the axon where stationary mitochondria are localized.The selection of the number of the demand sites should be interpreted as a means of simulating axonal morphology rather than discretization.It is important to acknowledge that in a real axon, the number of the demand sites will be substantially larger than 10.A narrow demand site is assumed to be located at the center of each compartment, and the compartment is assumed to surround it.
The primary motivation for employing a multi-compartment model is our desire to use the methodology for computing the particle age density distribution in the compartments developed in References 34, 35.We used a similar approach to simulate transport of dense core vesicles in Drosophila motoneurons. 36,37The demand sites in the axon terminal were numbered from the most proximal (site 1) to the most distal (site N) (Figure 1).
To characterize mitochondria concentration in the axon, we used the total length of mitochondria in a particular kinetic state (Figure 2) per unit axon length, measured in (μm of mitochondria)/μm.The dependent variables utilized in the model are given in Table S1 in Supporting Information.

|
Governing equations and the model for the mean age and age density distributions of mitochondria in the demand sites The total length of mitochondria is a conserved quantity.The equations stating the conservation of the total length of mitochondria in 3 N compartments (Figure 2) are stated using the methodology developed in Reference 31.In a vector form, these equations are where t is the time and C Ã is the state vector, following the terminology used in Reference 34.C Ã is a column vector; its components are defined as follows: It is convenient to introduce vector l defined by its components where l i is the total length of mitochondria per unit length of the axon in stationary, anterograde, or retrograde states in the compartment by the ith demand site (Table S1).
Model parameters are summarized in Table S2.In Equation (2) L i (i = 1,…, N) is the length of an axonal compartment around the ith demand site.Since the axonal length L ax was split into N compartments of equal length, L i ¼ L ax =N.
The first N components of the state vector represent the total length of mitochondria in stationary states, the next N components (N + 1,…, 2 N) represent the total length of mitochondria in anterograde states, and the last N components (2 N + 1,…, 3 N) represent the total length of mitochondria in retrograde states.The definition of vector u t ð Þ is given below in Equations ( 45) and (46).Matrix B (3 N, 3 N) is composed following Reference 31.It accounts for the internal mitochondria fluxes between the compartments, the external mitochondria flux entering the axon from the soma, and the mitochondria flux leaving the axon, part of which returns to the soma for degradation and part re-enters the axon (Figure 2).The analysis of mitochondria fluxes to and from the stationary, anterograde, and retrograde compartments surrounding the most proximal demand site gives equations for the following elements of matrix B: where k w is the kinetic constant characterizing the rate of mitochondria re-release from the stationary state into the mobile states.
where j denotes the flux of mitochondria length between the compartments and p s is the probability of mitochondria transitioning from a moving to the stationary state.
where j soma!1 is the total flux of anterograde mitochondria entering the most proximal demand site, which consists of mitochondria newly synthesized in the soma and mitochondria returning from their journey in the axon.Our model assumes that a portion of mitochondrial flux exiting the axon, αj 1!soma , returns to the axon while the other portion of the flux, 1 where ε is the portion of mitochondria that re-enter the anterograde component and 1 À ε ð Þis the portion of mitochondria that re-enter the retrograde component (Figure 2, Table S2).
By analyzing mitochondria fluxes to and from the compartments surrounding demand site i (i = 2,…, NÀ1), equations for the following elements of matrix B are obtained: The analysis of mitochondria fluxes to and from the compartments surrounding the most distal demand site gives equations for the following elements of matrix B: Other elements of matrix B, except for those given by Equations ( 4)- (34), are equal to zero.Equations ( 4)- (34) were supplemented by the following equations that simulate fluxes of mobile mitochondria between the demand sites.Equations for anterograde fluxes (Figure 2) are given by the following equations: Equations for retrograde fluxes (Figure 2) are given by the following equations: The N + 1th element of column vector u is given by the following equation: In Equation (45) j soma!1 is the total flux of anterograde mitochondria entering the most proximal demand site.It consists of mitochondria newly synthesized in the soma and mitochondria that just left the axon and re-enter the axon anterogradely.Equation (45) means that the flux of newly synthesized mitochondria entering the axon from the soma, u Nþ1 , is adjusted in our model depending on the flux of mitochondria that re-entered the axon after leaving it retrogradely, αj 1!soma .This is done to keep the total flux of mitochondria entering the axon (newly synthesized plus returning) constant and independent of α.
The only flux entering the compartmental system displayed in Figure 2 is the flux of newly synthesized mitochondria entering the N + 1th compartment.Therefore, all other elements of vector u are equal to zero: The aim of using Equations ( 1)-( 46) was to utilize advanced methods developed for compartmental systems 34,35,38,39 to analyze the average age and distribution of age density of mitochondria in demand sites at varying distances from the soma.A different approach to expressing the conservation equations for the overall length of mitochondria in the axon, as presented in Reference 28, is described in Section S1 of Supporting Information.Equations in Section S1 directly state mitochondria conservation in the demand site, rather than combining stationary, anterograde, and retrograde states into a single vector.Consequently, they are easier to interpret physically.However, they are not suitable for calculating the mean age and age density distributions of mitochondria, as those require stating the conservation equations in the vector form given by Equation (1).
We used the formulas reported in Reference 35 to obtain results for the steady-state situation (when the left-hand side of Equation ( 1) is equal to zero).The solution of Equation ( 1) for a steady-state situation is given by the following equation: In Equation ( 47), the superscript À1 on a matrix denotes the inverse of the matrix and the subscript ss denotes steady-state.
The mean ages and age density distributions of mitochondria in various demand sites were computed using the method described in References 34, 35.The mean ages of mitochondria in various compartments displayed in Figure 2 at steady-state are found from the following equation: In Equation (48), is the diagonal matrix with the components of the state vector, defined in Equation ( 2), on the main diagonal.These components (C Ã 1,ss , C Ã 2,ss , etc.) are calculated at steady-state.The left-hand side of Equation ( 48) is a column vector composed of the mean ages of mitochondria in various compartments: where superscript T denotes the transposed matrix.
The age densities of mitochondria in various compartments are given by the elements of the following vector: where e is the matrix exponential.Equations ( 47)-( 51) allow finding steady-state solutions directly.We implemented these equations using standard MATLAB operators, such as matrix inverse and matrix exponential.
We defined the following vectors, each of size N, which represent the mean ages of mitochondria in the stationary, anterograde, and retrograde states in the demand sites, respectively: 3 | RESULTS 3.1 | Total length of mitochondria per unit axon length in the demand sites and mitochondria fluxes between the demand sites Model parameters were chosen using the computational study of Reference 25, which found parameter values that optimized mitochondrial health.The parameter values used in our research are summarized in Table S2.
For α ¼ 0, p s ¼ 0:4, the distribution of the total length of mitochondria per unit length of the axon is uniform along the length of the axon (Figure S1a).The uniform distribution is expected for all parameters, regardless of α.This is because the model assumes that no mitochondria are destroyed while in the axon, although 1 À α ð Þportion of all mitochondria are assumed to be destroyed in the soma when they re-enter the soma from the axon.The uniform distribution displayed in Figure S1a is consistent with Reference 40 which noted that in axons the mitochondrial density is independent of the distance from the soma (Figure 1C in Reference 40).A constant mitochondrial density of 5 mitochondria/100 μm was also reported for growing neuronal processes in C. elegans. 41he return of 30% (α ¼ 0:3, p s ¼ 0:4) of exiting mitochondria back to the axon results in the same uniform distribution of mitochondria along the axon length (Figure S1a).This is because we assumed that the total flux of mitochondria entering the first demand site is the same for all cases computed in this paper.The flux is independent of what portion of returning mitochondria re-enter the axon (Equation ( 45)).If more mitochondria re-enter the axon (which corresponds to a greater value of α), the flux of newly synthesized mitochondria entering the first demand site in our model is reduced to keep the total flux of mitochondria entering the first demand site (newly synthesized plus returning, Figure 2) the same.
The reduction of parameter p s , which characterizes the rate of mitochondria transition to the stationary state (Figure 2), to 0.1, resulted in a reduction of the total length of mitochondria per unit length in the stationary state by approximately a factor of 4, which is due to fewer mitochondria transitioning to the stationary state (Figure S1b).The increase of the value of α to 0.3 did not change the distributions of the total length of mitochondria per unit axonal length (Figure S1b).This is because, according to our assumptions, the change of α does not affect the flux of mitochondria entering the most proximal demand site.
Anterograde and retrograde fluxes between the compartments are also independent of the values of α and p s , and are uniform along the axon length (Figure S2a,b).

| Mean age and age density distributions of mitochondria
In Figure 3, we use line graphs to depict the mean age of mitochondria versus the demand site number, which characterizes the distance from the soma (Figure 1).It should be noted, however, that per Reference 25, we simulate the demand sites as highly confined segments within the axon where stationary mitochondria are present.For p s ¼ 0:4 and α ¼ 0 the mean age of mitochondria increases gradually from the most proximal to the most distal demand site.The mean age of anterograde mitochondria is the smallest, the mean age of stationary mitochondria is in the middle, and the mean age of retrograde mitochondria is the greatest.The mean age of stationary mitochondria in the most distal demand site is approximately 21 hours (Figure 3a).Since the compartmental model does not consider the time it takes for mitochondria to transition between the compartments, 28 it is necessary to prove that the age of mitochondria predicted by the model is greater than the time it takes for mitochondria to travel from the soma to the axon tip without stops, L ax =v: This means that a compartmental model can be used to simulate the transport of mitochondria in relatively short axons, as one simulated in this paper, of the length of 10 4 μm.
The mean age for mitochondria for p s ¼ 0:4 and α ¼ 0 roughly follows the same trend (an increase from the most proximal to the most distal demand site), but it is now much greater.The mean age of stationary mitochondria in the most distal demand site is now approximately 36 hours (Figure 3b).This is because for α ¼ 0:3 the flux of mitochondria entering the axon includes mitochondria that leave the axon retrogradely (Figure 2).We thus successfully simulated the situation described in Reference 19 that suggests that older mitochondria appear in the axon due to return of mitochondria leaving the axon, which turn around and return to the axon in the anterograde component.It should be   noted that the real situation is even more complicated, since mitochondria returning to the axon are most likely repaired in the soma by fusion with newly synthesized mitochondria.Furthermore, moving mitochondria can fuse with stationary mitochondria residing in the axon and exchange mitochondrial proteins. 5,25,42This process will need to be simulated in future versions of the model.For a smaller value of p s (p s ¼ 0:1, α ¼ 0, Figure 3C), the mean age of mitochondria is significantly reduced.The mean age of stationary mitochondria in the most distal demand site is now approximately 10 h.This is because the probability of transitioning between moving states (anterograde or retrograde) and the stationary state (Figure 2) is now much smaller.It is interesting that although the mean age of anterograde and stationary mitochondria increases from the most proximal to the most distal demand site, the age of retrograde mitochondria now decreases from the most proximal to the most distal demand site (Figure 3C).This is because the exchange between moving and stationary mitochondria is now much less frequent (due to a smaller value of p s ), and mitochondria mostly remain in the same kinetic state (anterograde, stationary, or retrograde).Retrograde mitochondria age as they return from the most distal demand site (site 10) to the most proximal demand site (site 1), which explains why the age of retrograde mitochondria is the largest in the most proximal demand The increase of α (p s ¼ 0:1, α ¼ 0:3, Figure 3D) results in the increase of the mean age of mitochondria in all demand sites, but the trends stay the same.The mean age of the stationary mitochondria in the most distal demand site is now   approximately 20 h.The increase in the mean age is due to the return of 30% of mitochondria exiting the most proximal demand site, which turn around and come back to the axon (Figure 2).The age density of mitochondria is important to characterize the ranges of mitochondria ages in different kinetic states (anterograde, stationary, and retrograde) in different demand sites.Figure 4A compares the age densities of mitochondria computed for the base case (p s ¼ 0:4, α ¼ 0, red lines) with age densities of mitochondria computed for the case of p s ¼ 0:4, α ¼ 0:3, blue lines.The age densities of anterograde mitochondria are more narrow (Figure 4A, see also Figure 4B that shows a magnified view of Figure 4A).Notably, the distribution of age densities widens for stationary mitochondria (Figure 4C) and becomes even wider for retrograde mitochondria (Figure 4D).The shift of the age density curves towards older mitochondria, observed when α is increased from 0 to 0.3 (Figure 4D), can be attributed to the reentry of previously aged retrograde mitochondria back into the axon.The presence of older mitochondria is consistent with the results reported in Reference 19.The peaks on the curves displaying the age density distributions shifted toward older mitochondria for more distal demand sites (Figure 4).
The reduction of the rate of mitochondria transition to the stationary state (characterized by p s ) makes the age density distribution of mitochondria (p s ¼ 0:1, α ¼ 0, blue curves, Figure 5) quite complex.The most noticeable change is the bimodal distribution that is visible in Figure 5D in the solid blue curve displaying the age density of retrograde mitochondria in the most proximal demand site (site 1).The first peak occurs at approximately 1 h.It is explained by mitochondria that recently entered the axon, but then transitioned to the stationary state, and then to the retrograde state (Figure 2).The second peak occurs at approximately 11 h.It is explained by older mitochondria that travelled to the axon tip and then returned to the most proximal site.
The increase of the value of α (more mitochondria return to recirculation in the axon) while keeping p s low (small transition rate to the stationary state) makes the tails in age density distributions longer, which means the presence of older mitochondria in the axon, see Reference 35 (p s ¼ 0:1, α ¼ 0:3, blue curves in Figure S3d).Interestingly, the first peak on the solid blue curve in Figure S3d, displaying the age density of retrograde mitochondria in site 1 (at approximately 1 h), has almost disappeared; it degraded nearly to the horizontal slope of the curve (compare Figure S3d with Figure 5D).All age density distributions in Figures 4, 5, S3, and in particular in Figure S3d, are skewed right.This indicates that although the mean ages of mitochondria are of the order of 20 h (see Figure 3D, which shows mean ages of mitochondria that correspond to the age density distributions depicted by blue lines in Figure S3d), there are also much older mitochondria present in the axon, which stay in the axon several times longer than 20 h, probably completing several circulations or staying a long time in the stationary state before returning to the mobile pool.
We have confirmed that the bimodal age density distribution of retrograde mitochondria is not exclusive to 10 demand sites.In Figures S4 and S5 in Supporting Information, we analyzed the age density distributions of mitochondria in an axon consisting of five demand sites (N = 5).The bimodal distributions observed for retrograde mitochondria in the most proximal demand site are visible in Figures S4d and S5d.Our hypothesis suggests that the bimodal distribution arises from the presence of retrograde mitochondria with two distinct ages: (i) relatively young mitochondria that stop and become retrograde, and (ii) relatively old mitochondria that travel to the end of the axon, reverse their direction, and return.The difference in ages between these two types of mitochondria is most prominent in the most proximal demand site, resulting in the occurrence of a bimodal distribution.In more distant demand sites, the bimodal pattern may become less pronounced.

| Sensitivity of the mean age of mitochondria in demand sites to p s and α
4][45][46] For example, the sensitivity coefficient of the mean age of resident mitochondria in demand sites to parameter α at steady-state (ss) was calculated as follows: where Δα ¼ 10 À1 α (we tested the accuracy by using various step sizes).
To make sensitivity coefficients independent of the magnitude of the parameter whose sensitivity was tested, we calculated non-dimensional relative sensitivity coefficients, 44,47 which were defined as (for example): The dimensionless sensitivity of the mean age of mitochondria to the probability that mobile mitochondria would transition to a stationary state at a demand site, p s , is positive in all demand sites (Figure 6A).This is because mitochondria spend more time in the stationary state when the stopping probability is larger.
The dimensionless sensitivity of the mean age of mitochondria to the portion of mitochondria that return to the axon after exiting the axon, α, is also positive in all demand sites (Figure 6B).This is because older mitochondria, which already circulated in the axon, re-enter the axon when the value of α is larger than zero.The number of older mitochondria in the axon thus increases.This is consistent with the results depicted in Figure 3 (compare the results in Figure 3A and Figure 3B, and also the results in Figure 3C and Figure 3D).It is interesting to note that the mean age of stationary mitochondria is approximately 10 times more sensitive to α than to p s (compare Figure 6B with Figure 6A).We examined whether this conclusion remains consistent regardless of the number of demand sites.The findings depicted in Figure S6, computed for five demand sites (N = 5), reveal that the sensitivity to α is independent of the number of compartments, while the sensitivity to p s scales linearly with the number of compartments in the axon.This outcome arises from the fact that p s represents the stopping probability in each compartment.Notably, this result implies that for a very long axon with many compartments, the age of mitochondria may be highly sensitive to the stopping probability.

| DISCUSSION, LIMITATIONS OF THE MODEL, AND FUTURE DIRECTIONS
We developed a model that accounts for the return of a portion of mitochondria (characterized by parameter α) back to the axon after mitochondria complete circulation in the axon.We investigated how the mean age and age density distribution of mitochondria depend on α.We also investigated the dependence of the same quantities on p s .
Bimodal age density distributions are found for a smaller value of p s (p s ¼ 0:1Þ.The peak corresponding to older mitochondria on the curve displaying the age density distribution of retrograde mitochondria is explained by mitochondria that traveled to the tip of the axon and then traveled back in the retrograde component.The peak corresponding to younger mitochondria on the curve displaying the age density distribution of retrograde mitochondria is explained by mitochondria that transitioned from the anterograde state to the state occupied by stationary mitochondria and then to the retrograde state.We found that the mean age of stationary mitochondria is very sensitive to parameter α.It is an order of magnitude less sensitive to parameter p s .We also investigated how the sensitivity of the mean age of stationary mitochondria depends on the number of compartments in the axon.The sensitivity to the stopping probability, p s , increases proportionally with the number of compartments in the axon.This suggests that in an extremely long axon with a large number of compartments, the age of mitochondria could be highly influenced by the stopping probability.
The prediction of a broader age density distribution for retrograde mitochondria compared to anterograde mitochondria could in principle be accessible with experimental measurements.
Future research should aim to explore the feasibility of analytically solving equations of our model to obtain explicit expressions for the age densities of mitochondria in anterograde, stationary, and retrograde kinetic states.Future modeling work should prove that the time-scale for transitioning between compartments in our model is comparable to the time-scale for mixing within an individual compartment.To investigate the time-scale for mixing within an individual compartment, it may be necessary to develop a model that would treat mitochondria as individual particles. 25,48he utilization of a discrete model would provide much more detailed information on mitochondria trajectories and age.A promising extension of the current approach would also involve incorporating the stochastic nature of mitochondrial transport.Furthermore, in future research, a model of fusion/fission of mitochondria should be developed.This is especially important because fusion with newly synthesized mitochondria contributes to the replenishment of older mitochondria with new mitochondrial proteins. 19Mitochondria removal via mitophagy 49,50 also needs to be incorporated in future models.Some mitochondria degradation may occur in axons locally. 51,52The effect of this possibility should be incorporated into future models.Future research should also investigate the effect of domain discretization (using a different number of compartments) on the predictions of the developed model.

F
I G U R E 1 (a) Schematic diagram of a neuron with an axon that is assumed to contain N energy demand sites.Demand sites are numbered starting with the most proximal (site 1) to the most distal (site N).Figure is generated with the aid of Servier Medical Art, licensed under a creative commons attribution 3.0 generic license, http://Smart.servier.com.F I G U R E 2 A diagram of a compartmental model showing transport in the compartments surrounding the demand sites in the axon.Arrows show mitochondria exchange between the transiting (anterograde and retrograde) states in adjacent demand sites, mitochondria capture into the stationary state and re-release from the stationary state.A portion of mitochondria re-enter the axon, and the rest return to the soma for degradation.Newly synthesized mitochondria entering the axon from the soma are assumed to have zero age.

F
I G U R E 4 (a) Age density of anterogradely moving mitochondria in various demand sites.(b) Similar to Figure 4A, this graph magnifies a specific age range of [0 10 hours] of the x-axis.(c) Age density of mitochondria in the stationary state in various demand sites.(d) Age density of retrogradely moving mitochondria in various demand sites.Two sets of parameter values: p s ¼ 0:4, α ¼ 0 (base case) and p s ¼ 0:4, α ¼ 0:3.Number of demand sites N = 10.

F
I G U R E 5 (a) Age density of anterogradely moving mitochondria in various demand sites.(b) Similar to Figure 5A, this graph magnifies a specific age range of [0 10 hours] of the x-axis.(c) Age density of mitochondria in the stationary state in various demand sites.(d) Age density of retrogradely moving mitochondria in various demand sites.Two sets of parameter values: p s ¼ 0:4, α ¼ 0 (base case) and p s ¼ 0:1, α ¼ 0. Number of demand sites N = 10.

F
I G U R E 6 (a) Dimensionless sensitivity of the mean age of mitochondria to the probability that mobile mitochondria would transition to a stationary state at a demand site, p s , versus the site number.Computations were performed with Δp s ¼ 10 À1 p s .A close result was obtained for Δp s ¼ 10 À2 p s .(b) Dimensionless sensitivity of the mean age of mitochondria to the portion of mitochondria that return to the axon after exiting the axon, α, versus the site number.Computations were performed with Δα ¼ 10 À1 α.A close result was obtained for Δα ¼ 10 À2 α.Sensitivities are analyzed around p s ¼ 0:1 s À1 , α ¼ 0:3.Number of demand sites N = 10.