PMCCPMCCPMCC

Search tips
Search criteria 

Advanced

 
Logo of plosonePLoS OneView this ArticleSubmit to PLoSGet E-mail AlertsContact UsPublic Library of Science (PLoS)
 
PLoS One. 2012; 7(12): e51909.
Published online 2012 December 14. doi:  10.1371/journal.pone.0051909
PMCID: PMC3522595

Integrated Model of Metabolism and Autoimmune Response in β-Cell Death and Progression to Type 1 Diabetes

Andrew R. Dalby, Editor

Abstract

Progression to type 1 diabetes is characterized by complex interactions of environmental, metabolic and immune system factors, involving both degenerative pathways leading to loss of pancreatic β-cells as well as protective pathways. The interplay between the degenerative and protective pathways may hold the key to disease outcomes, but no models have so far captured the two together. Here we propose a mathematical framework, an ordinary differential equation (ODE) model, which integrates metabolism and the immune system in early stages of disease process. We hypothesize that depending on the degree of regulation, autoimmunity may also play a protective role in the initial response to stressors. We assume that β-cell destruction follows two paths of loss: degenerative and autoimmune-induced loss. The two paths are mutually competing, leading to termination of the degenerative loss and further to elimination of the stress signal and the autoimmune response, and ultimately stopping the β-cell loss. The model describes well our observations from clinical and non-clinical studies and allows exploration of how the rate of β-cell loss depends on the amplitude and duration of autoimmune response.

Introduction

Type 1 diabetes (T1D) is an autoimmune disease characterized by a relatively long symptom-free period that precedes the flair-up of clinical signs of the disease. T1D is caused by progressive loss of insulin-secreting capacity of the β-cells and, finally, by selective death of these cells in the over 1 million islets of Langerhans in the human pancreas. The incidence of T1D among children and adolescents has increased markedly in the Western countries over the recent decades, with the annual incidence continuing to rise at an accelerating pace [1]. Although approximately 80% of subjects with T1D carry defined risk-associated genotypes at the HLA locus, only 3–7% of the carriers of such genetic risk markers develop the disease [2]. Seroconversion to islet autoantibody positivity has been the first detectable signal implicating initiation of autoimmunity and progression towards diabetes [3]. However, although seroconversion to autoantibody positivity precedes the clinical disease by months to several years, the time point at which seroconversion occurs may already be too late for therapeutic approaches aimed at preventing progression to overt diabetes. As long as the initiators of the autoimmune response remain unknown and the mechanisms supporting progression towards β-cell failure are poorly understood, the estimation of disease risk, time of disease presentation in genetically susceptible individuals as well as discovery of effective prevention will remain a challenge.

It is generally accepted that the primary role of the immune system is in protection of the body against foreign pathogens which requires the inhibition of immune response against self. The classical view known as “clonal selection theory” states that potentially pathogenic T cells such as the cells that react specifically to autoantigens should not be present in healthy individuals [4], [5]. However, there is also evidence to the contrary; suggesting autoimmunity may also have a protective role in specific circumstances. For example, it has been shown that autoimmunity may also be a physiological response with a protective role following an insult to the central nervous system [6], [7], [8].

Our metabolomics studies in children who later progressed to T1D (progressors) identified specific metabolic profile preceding the first autoimmune response [9]. Surprisingly, following the autoimmune response this metabolic profile was largely corrected towards the normal levels found in children who did not progress to the disease. This led us to hypothesize that (a) autoimmunity results from physiological adaptation triggered by metabolic stress, and that (b) the subsequent T1D in some individuals is due to defective regulation of the immune response (e.g., its amplitude or timing). Subsequent studies in non-obese diabetic (NOD) mice supported this view as well as identified specific “protective pathways” in mice who seroconverted to insulin autoantibody (IAA) positivity but did not develop autoimmune diabetes [10].

Given the time-dependent complexity of T1D pathogenesis, including interactions across multiple physiological systems, there is a clear case for mathematical modeling in order to capture the complex interplay of the factors involved in the disease process. However, only few attempts to model the disease process have been made. The so-called Copenhagen model [11], [12], [13] focuses on the interplay between activated macrophages, T-helper (Th)-lymphocytes and target cells in the early phases of T1D pathogenesis. According to this model, the mechanism of β-cell destruction as induced by cytokines occurs prior to β-cell destruction by T-lymphocyte mediated mechanisms. This simple model already displayed remarkable complexity and led the authors to conclude that “onset of T1D is due to a collective, dynamical instability, rather than being caused by a single etiological factor” [13]. However, the model does not include protective pathways, or metabolism in general, as factors protecting from or contributing to the disease.

Nevo et al. introduced a model of autoimmune response where autoimmunity is viewed as a defense mechanism against tissue-specific degenerative processes [14]. The degenerative process of tissue loss is triggered by the primary insult and it further activates alarm signals which provoke an autoimmune response [15], [16]. The model suggests that anti-self response after being triggered competes with self-perpetuating degenerative tissue loss. Depending on timing and intensity of the autoimmune response, there are three different outcomes of the competition: (a) if the autoimmune response occurs too late or is too weak, it is not able to stop the degenerative process, consequently the self-destruction dominates and tissue degeneration continues; (b) if the autoimmune response fails to stop, regardless of its intensity, the tissue loss continues and ultimately leads to autoimmune disease; (c) autoimmune response can also be protective if it starts early enough and shuts off at the right time, leading to diminishment of the self-degenerative process.

Herein we propose an integrated model of metabolism and autoimmune response to study how they together contribute to progressive β-cell loss and T1D. The model builds on earlier models of β-cell loss [13] and autoimmune response [14] as well as on our longitudinal metabolomics and autoantibody data from children who progressed to T1D in a prospective birth cohort [9] and the data from NOD mouse [10].

Results and Discussion

Our primary aim is to model the interaction between metabolism and immune system in the early stages of T1D. In summary, β-cell destruction follows two paths of loss: degenerative loss and autoimmune-induced loss. First, the degenerative process of β-cell loss is triggered, e.g., by an infection [17], leading to activation of autoimmune response via secretion of a danger signal. The autoimmune response stimulates the autoimmune-induced β-cell loss. The two paths of loss are mutually competing, leading to termination of the degenerative path of loss and further to elimination of the danger signal and the autoimmune response. As a consequence, the autoimmune induced path of loss is terminated as well, stopping the β-cell loss.

A simplified Copenhagen model [13] is assumed for the degenerative path. The process of degenerative loss starts with release of proteins triggered by the environmental factors attacking the β-cells. These proteins are detected by macrophages in the islets of Langerhans, leading to formation of activated macrophages. Activated macrophages release signal molecules (e.g., cytokines) and immunogenic danger signals. Danger (alerting) signals activate the autoimmune response which further triggers the autoimmune-induced process of β-cell loss. Production of cytokines by activated macrophages is on the other hand associated with release of β-cell antigenic proteins. Protective pathway is supposed to be the one that is associated with the seroconversion to autoantibody positivity [9], [10]. Serum lysophosphatidylcholine (lysoPC) is selected to represent such pathway because it is inversely related to insulin autoantibody positivity (IAA+) in NOD progressors and non-progressors to autoimmune diabetes [10], thus reflecting altered association of metabolic and immune system status in disease pathogenesis.

Model variables are a combination of variables from Copenhagen model [13] and from protective autoimmunity model [14]: An external file that holds a picture, illustration, etc.
Object name is pone.0051909.e001.jpg is the amount of macrophages, An external file that holds a picture, illustration, etc.
Object name is pone.0051909.e002.jpg is the amount of activated macrophages, An external file that holds a picture, illustration, etc.
Object name is pone.0051909.e003.jpg is the amount of beta cell antigenic proteins, as in [13]. An external file that holds a picture, illustration, etc.
Object name is pone.0051909.e004.jpg is the strength of autoimmune response and it is introduced in the same sense as autoimmune response variable An external file that holds a picture, illustration, etc.
Object name is pone.0051909.e005.jpg in [14]. An external file that holds a picture, illustration, etc.
Object name is pone.0051909.e006.jpg is the concentration of metabolite of protective pathway which is related to autoimmune induced path of loss. This variable has the same meaning as variable An external file that holds a picture, illustration, etc.
Object name is pone.0051909.e007.jpg from [14] which represents the population density of cells that undergo the immune mediated (i.e. positive) path of loss. We suggest biological interpretation to this path as a protective metabolic pathway which due to its protectiveness has a positive (defensive) role in beta cell loss. Schematic representation of the model is given in Figure 1, with the model variables and parameters shown in Table 1. Mathematical model is described in detail in Materials and Methods section.

Figure 1
Schematic representation of the model.
Table 1
Description of the model variables and parameters.

Since the experimental longitudinal data on β-cell mass are not available, we use experimental longitudinal data of glucose concentrations from NOD mice [10] and apply the model by [18] in order to estimate β-cell mass and select the model parameters (Figures 2A and 2B). The mice are divided into four groups: (a) progressors to autoimmune diabetes who were IAA+ when measured at age of 8 weeks, (b) non-progressors who were IAA+ at 8 weeks of age (“protected” mice), (c) progressors who were IAA- at 8 weeks of age, and (d) non-progressors who were IAA- at 8 weeks of age. As a verification of our qualitatively estimated parameters model simulations are performed with parameters values of An external file that holds a picture, illustration, etc.
Object name is pone.0051909.e008.jpg (recruitment rate of macrophages by activated i.e. inflow rate induced by activated macrophages) and An external file that holds a picture, illustration, etc.
Object name is pone.0051909.e009.jpg (beta cell apoptosis induced by activated macrophages i.e. rate of creation of beta cell antigenic proteins) for NOD mice which were estimated based on experiments and literature [19]. Simulated profiles obtained with our qualitatively estimated set of parameters for the case of IAA+ progressors fit well the simulations reproduced with two estimated parameters from [19]. A comparison of beta cell profiles is shown in Figure 2C.

Figure 2
Model predictions of age-dependent glucose profiles in disease progressors and non-progressors and model validation based on prediction of β-cell mass.

The differential equations determine the values of the simulation results that depend on altogether nine parameters. Among these, An external file that holds a picture, illustration, etc.
Object name is pone.0051909.e012.jpg, An external file that holds a picture, illustration, etc.
Object name is pone.0051909.e013.jpg and An external file that holds a picture, illustration, etc.
Object name is pone.0051909.e014.jpg affect the shape of the simulated glucose curve (see equations 6.1–6.3 in Materials and Methods). Values for these parameters were set by fitting the model outcomes to the experimental data from [10] using multidimensional unconstrained nonlinear minimization based on Nelder-Mead method [20]. Since the number of experimental points in available datasets is quite small (17 points for IAA+ progressors, 30 points for IAA- progressors, 21 points for IAA+ non-progressors and 31 points for IAA- non-progressors) one cannot properly draw conclusions about the exact behavior of data (i.e. about the shapes of the data curves). We rather limit ourselves to interpret and to fit the trends of data: how fast the concentration of glucose tends to increase depending on conditions of IAA positivity and progressivity or non to the disease. Small oscillations in non progressors data are neglected and the trend in these data sets are interpreted as slowly increasing (almost constant). Fitting with such a small but only available data set is not precise enough which can be also seen from the goodness-of-fit An external file that holds a picture, illustration, etc.
Object name is pone.0051909.e015.jpg values (e.g. IAA- non-progressors An external file that holds a picture, illustration, etc.
Object name is pone.0051909.e016.jpg is quite low). According to our model the amount of macrophages increases while the amount of activated macrophages decreases with time (Figure 3). This is in agreement with the Copenhagen model [13], which explains that after an initial time during which the number of macrophages drops down because of antigen uptake, it increases later on because of activation-induced inflow and deactivation. The number of activated macrophages according to Copenhagen model after an initial period during which it grows because antigens are taken up, it declines later on when the antigens have been swept away.

Figure 3
Simulated amounts of macrophages and of activated macrophages.

As part of qualitative evaluation of the model, intensity profiles of autoimmune response were compared to the ones reported in [14]. Up to the authors knowledge these simulated profiles of autoimmunity are the only ones found in literature and there are no experimental profiles available. The protective autoimmune response simulated by our model appears to agree with the simulated protective autoimmune response reported in [14]. Comparison between simulated autoimmune responses for different values of parameters An external file that holds a picture, illustration, etc.
Object name is pone.0051909.e017.jpg and An external file that holds a picture, illustration, etc.
Object name is pone.0051909.e018.jpg (equation 2 in Materials and Methods) is shown in Figure 4. Autoimmune response simulated by setting parameter An external file that holds a picture, illustration, etc.
Object name is pone.0051909.e019.jpg to be higher than the one in protective autoimmunity can be interpreted as autoimmune response with premature shut-off (Figure 4B). Autoimmune response simulated by setting parameter An external file that holds a picture, illustration, etc.
Object name is pone.0051909.e020.jpg to be lower than the one in protective autoimmunity corresponds to delayed-onset autoimmunity (Figure 4C).

Figure 4
Different time courses of autoimmune response.

Next, we studied time profiles of the protective pathway and β-cell loss (Figure 5). The influence of intensity of autoimmune response on activation of protective pathway is shown in Figure 5A. Once autoimmune response is triggered, protective pathway behaves differently depending on the intensity of autoimmune response. Autoimmune response which is strong enough leads to creation of effective protective pathway and weak autoimmune response provokes defective protective pathway. When autoimmune response is stronger, the concentration of metabolite of protective pathway is higher. This is in agreement with the changes of serum lysoPC in NOD mice in the context of IAA positivity and progression to autoimmune diabetes [10], where it has been observed that IAA+ positive mice with high lysoPC were protected from autoimmune diabetes. When comparing autoimmune progressors (defective protective pathway) and non-progressors (effective protective pathway), the model predicted that rate of β-cell loss is markedly slowed in the presence of protective pathway (Figure 5B). It is assumed that all the underlying metabolic mechanisms are implicitly included in the protective pathway so as the role of T cells in the pathogenesis of T1D.

Figure 5
Protective pathway and progression to autoimmune diabetes.

While inclusion of the protective pathway in the model of β-cell loss shows promise for the study of mechanisms behind protection from T1D, at present it is generalized by taking one enzymatic reaction for the whole protective pathway. This is clearly an over-simplification. In the future, the protective path equations need to be developed for a realistic pathway and to take into account explicitly the role of T cells in progression to T1D. For example, in our earlier study of NOD mice [10] we identified specific inflammatory pathways in isolated pancreatic islets which were up-regulated in IAA+ non-progressors. Among those, product of IL-4 pathway is in fact known to prevent diabetes in NOD mouse [21]. The serum lysoPC used in the present study is an inflammatory marker [22] and its up-regulation in IAA+ non-progressors may therefore reflect the activation of protective inflammatory pathways at specific stage following the autoimmune response. As another candidate protective pathway, islet pathways related to mitochondrial function such as TCA cycle, branched chain amino acid catabolism, beta oxidation and oxidative phosphorylation, were down-regulated in IAA+ non-progressors. These pathways lead to reduced production of reactive oxygen species (ROS) [10]. Since ROS generated by mitochondria plays an important role in the release of pro-apoptotic proteins which can trigger apoptosis [23], decreased ROS production could thus be linked to prevention of loss of β-cell functionality in IAA+ non-progressors. Incorporation of more detailed pathways such as those described above into the specific parts of the model may identify specific molecular targets behind the control of the protective pathway.

In summary, herein we presented the first model of β-cell loss in progression to T1D which also includes the protective pathway; with the model parameters fitted using the experimental data from NOD mouse. The model predicts that the autoimmune response, if properly tuned, can also have a protective role, leading to the reduced rate of β-cell loss. At present, the model is qualitative but it already allows exploration of how the rate of β-cell loss depends on the amplitude and duration of autoimmune response, which in turns depends on the activation of protective pathways. As such, it is also a useful conceptual tool for the study of how immune system interacts with other physiological systems in progression to autoimmune disease.

Materials and Methods

Mathematical Model

In the degenerative path of loss, () equations from Copenhagen model [13] are applied. These equations describe changes in numbers of macrophages (An external file that holds a picture, illustration, etc.
Object name is pone.0051909.e023.jpg), activated macrophages (An external file that holds a picture, illustration, etc.
Object name is pone.0051909.e024.jpg) and β-cell antigenic proteins (An external file that holds a picture, illustration, etc.
Object name is pone.0051909.e025.jpg), where number of β-cell antigenic proteins corresponds to the amount of β-cells. Since the death of β-cells leads to release of proteins, it can be assumed that the rate of death of β-cells is proportional to the amount of activated macrophages [13]. There is therefore no need to explicitly include the β-cell population in the equations because release of proteins is directly related to the number of activated macrophages. Three ordinary differential equations for this path are:

equation image
(1.1)
equation image
(1.2)
equation image
(1.3)

Parameters An external file that holds a picture, illustration, etc.
Object name is pone.0051909.e029.jpgare rate parameters of different mechanisms included in progressive loss of β-cells. Their description is given in Table 1 and it follows the ones given in [13].

Equation for autoimmune response is inspired by the equation of the mean level of immune stimulation [14] which introduces a sigmoid function of the spatially accumulated alerting signal with a treshold value below which the induction of immune stimulation is negligible. Variables in our model are assumed as functions of time only. Danger (alerting) signal is not explicitly included but it is assumed that activated macrophages induce autoimmune response via danger signals. Therefore instead of sigmoid function of spatially accumulated danger signal we introduce a simplified function of activated macrophages An external file that holds a picture, illustration, etc.
Object name is pone.0051909.e030.jpg to the equation for autoimmune response. This function accounts for two limit values An external file that holds a picture, illustration, etc.
Object name is pone.0051909.e031.jpg (activation value) and An external file that holds a picture, illustration, etc.
Object name is pone.0051909.e032.jpg (shut-off value). The function is chose so that below An external file that holds a picture, illustration, etc.
Object name is pone.0051909.e033.jpg autoimmune response is not activated (analogously to the threshold value of the sigmoid function from [14]) and above An external file that holds a picture, illustration, etc.
Object name is pone.0051909.e034.jpg autoimmune response is closed (Figure 6). This equation is given as

equation image
(2)
Figure 6
Schematic presentation of function used in autoimmune response equation.

Because of the assumption that autoimmune response is created from activated macrophages and the hypothesis that rate of death of β-cell depends on the strength of autoimmune response, the amplitude of the autoimmune response is taken to be proportional to the product of parameters An external file that holds a picture, illustration, etc.
Object name is pone.0051909.e036.jpg (rate of activation of macrophages) and An external file that holds a picture, illustration, etc.
Object name is pone.0051909.e037.jpg(rate of creation of β-cell antigenic proteins) in the form An external file that holds a picture, illustration, etc.
Object name is pone.0051909.e038.jpg. An external file that holds a picture, illustration, etc.
Object name is pone.0051909.e039.jpg is decay of autoimmune response. Danger signal is not introduced as a separate variable but it is implicitly included in the equation which describes autoimmune response.

Protective pathway is simplified by assuming only one enzymatic reaction and one relevant metabolite for the whole metabolic pathway. The corresponding equation for the protective pathway is motivated by differential equations for the concentrations of metabolites given in [24] and the equation of the rate of transition of positive path [14]. This equation therefore takes the following form

equation image
(3)

where An external file that holds a picture, illustration, etc.
Object name is pone.0051909.e043.jpg is the reaction rate constant of the enzymatic reaction of the protective pathway,An external file that holds a picture, illustration, etc.
Object name is pone.0051909.e044.jpg is the concentration of the initial metabolite in the reaction and An external file that holds a picture, illustration, etc.
Object name is pone.0051909.e045.jpg is the concentration of the observed metabolite which represents the protective pathway and An external file that holds a picture, illustration, etc.
Object name is pone.0051909.e046.jpg is the decay rate. The influence of autoimmune response immune-mediated path of loss is expressed through multiplication with An external file that holds a picture, illustration, etc.
Object name is pone.0051909.e047.jpg.Immune-mediated path is triggered after activation of autoimmune response i.e. above activation threshold value An external file that holds a picture, illustration, etc.
Object name is pone.0051909.e048.jpg. Protective metabolic pathway is represented by serum lysoPC metabolite, which is upregulated in IAA+ non-progressor NOD mice [10].

Complete system of equations for the model is resumed below:

equation image
(4.1)
equation image
(4.2)
equation image
(4.3)
equation image
(4.4)
equation image
(4.5)

Description of all the parameters and variables of the system of five ordinary differential equations is given in Table 1.

Connecting β-cell Mass with Glucose Levels

The only available experimental data associated with β-cell mass are glucose concentration data. In order to connect β-cell mass with glucose levels an additional part of the model is developed that is based on combination between Copenhagen model equations (equations 4.1–4.3) and β-cell mass, insulin and glucose kinetics model equations [18]:

equation image
(5.1)
equation image
(5.2)
equation image
(5.3)

where An external file that holds a picture, illustration, etc.
Object name is pone.0051909.e057.jpg is the glucose concentration, An external file that holds a picture, illustration, etc.
Object name is pone.0051909.e058.jpg is the insulin concentration and An external file that holds a picture, illustration, etc.
Object name is pone.0051909.e059.jpg is the mass of β-cells. Short description of variables and parameters used in equations 5 according to [18] can be found in Table 2. According to [18] glucose and insulin dynamics are fast relative to An external file that holds a picture, illustration, etc.
Object name is pone.0051909.e060.jpgcell mass dynamics. Therefore the An external file that holds a picture, illustration, etc.
Object name is pone.0051909.e061.jpg cell mass, insulin and glucose model can be decomposed into fast (An external file that holds a picture, illustration, etc.
Object name is pone.0051909.e062.jpg,An external file that holds a picture, illustration, etc.
Object name is pone.0051909.e063.jpg) and slow (An external file that holds a picture, illustration, etc.
Object name is pone.0051909.e064.jpg cell mass) subsystems. To study the slow subsystem it is assumed that the fast subsystem is at steady state. The right side of equation 4.3., where parameter A represents β-cell mass, is equivalent with the right side of equation 5.3, resulting in expressing variable An external file that holds a picture, illustration, etc.
Object name is pone.0051909.e065.jpg (activated macrophages) in terms of variable An external file that holds a picture, illustration, etc.
Object name is pone.0051909.e066.jpg, equation (6.2). When evolving equation (6.2) the steady states of glucose and insulin are assumed according to the assumption made in [18], as explained above. Consequently, resulting system of three ordinary differential equations is obtained:

equation image
(6.1)
equation image
(6.2)
equation image
(6.3)
Table 2
Description of the variables and parameters from the β-cell model.

For more detail about the original model of coupled β-cell mass, insulin and glucose dynamics one should see [18]. We find such model suitable for needs of our model because it gives a way to obtain shape of β-cell mass from the profile of glucose which is reproduced from experimental data.

The whole model thus consists of two parts: the “main” part expressed through equations 4.1–4.5 and the “auxiliary” part expressed through equations 6.1–6.3. There are 8 ordinary differential equations in total and 6 variables. Variables An external file that holds a picture, illustration, etc.
Object name is pone.0051909.e070.jpg and An external file that holds a picture, illustration, etc.
Object name is pone.0051909.e071.jpg are connecting the two parts of the model.

Funding Statement

The study was supported by Academy of Finland (Centre of Excellence in Molecular Systems Immunology and Physiology Research 2012–2017, Decision No. 250114, http://www.symmys.fi/). The funder had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

References

1. Patterson CC, Dahlquist GG, Gyürüs E, Green A, Soltész G (2009) Incidence trends for childhood type 1 diabetes in Europe during 1989–2003 and predicted new cases 2005–20: a multicentre prospective registration study. Lancet 373: 2027–2033. [PubMed]
2. Achenbach P, Bonifacio E, Koczwara K, Ziegler A-G (2005) Natural history of type 1 diabetes. Diabetes 54: S25–31. [PubMed]
3. Kimpimaki T, Kulmala P, Savola K, Kupila A, Korhonen S, et al. (2002) Natural history of {beta}-cell autoimmunity in young children with increased genetic susceptibility to type 1 diabetes recruited from the general population. J Clin Endocrinol Metab 87: 4572–4579. [PubMed]
4. Burnet M (1959) Auto-immune disease. II. Pathology of the immune response. Br Med J 2: 720–725. [PMC free article] [PubMed]
5. Burnet M (1959) Auto-immune disease. I. Modern immunological concepts. Br Med J 2: 645–650. [PMC free article] [PubMed]
6. Yoles E, Hauben E, Palgi O, Agranov E, Gothilf A, et al. (2001) Protective autoimmunity is a physiological response to CNS trauma. J Neurosci 21: 3740–3748. [PubMed]
7. Schwartz M, Shaked I, Fisher J, Mizrahi T, Schori H (2003) Protective autoimmunity against the enemy within: fighting glutamate toxicity. Trends Neurosci 26: 297–302. [PubMed]
8. Moalem G, Leibowitz-Amit R, Yoles E, Mor F, Cohen IR, et al. (1999) Autoimmune T cells protect neurons from secondary degeneration after central nervous system axotomy. Nat Med 5: 49–55. [PubMed]
9. Oresic M, Simell S, Sysi-Aho M, Nanto-Salonen K, Seppanen-Laakso T, et al. (2008) Dysregulation of lipid and amino acid metabolism precedes islet autoimmunity in children who later progress to type 1 diabetes. J Exp Med 205: 2975–2984. [PMC free article] [PubMed]
10. Sysi-Aho M, Ermolov A, Gopalacharyulu PV, Tripathi A, Seppanen-Laakso T, et al. (2011) Metabolic regulation in progression to autoimmune diabetes. PLoS Comput Biol 7: e1002257. [PMC free article] [PubMed]
11. Nerup J, Mandrup-Poulsen T, Helqvist S, Andersen HU, Pociot F, et al. (1994) On the pathogenesis of IDDM. Diabetologia 37 Suppl 2S82–89. [PubMed]
12. Nerup J, Mandrup-Poulsen T, Molvig J, Helqvist S, Wogensen L, et al. (1988) Mechanisms of pancreatic beta-cell destruction in type I diabetes. Diabetes Care 11 Suppl 116–23. [PubMed]
13. Freiesleben De Blasio B, Bak P, Pociot F, Karlsen AE, Nerup J (1999) Onset of type 1 diabetes: a dynamical instability. Diabetes 48: 1677–1685. [PubMed]
14. Nevo U, Golding I, Neumann AU, Schwartz M, Akselrod S (2004) Autoimmunity as an immune defense against degenerative processes: a primary mathematical model illustrating the bright side of autoimmunity. J Theor Biol 227: 583–592. [PubMed]
15. Nevo U, Kipnis J, Golding I, Shaked I, Neumann A, et al. (2003) Autoimmunity as a special case of immunity: removing threats from within. Trends Mol Med 9: 88–93. [PubMed]
16. Matzinger P (1998) An innate sense of danger. Semin Immunol 10: 399–415. [PubMed]
17. Knip M, Veijola R, Virtanen SM, Hyöty H, Vaarala O, et al. (2005) Environmental triggers and determinants of type 1 diabetes. Diabetes 54: S125–136. [PubMed]
18. Topp B, Promislow K, deVries G, Miura RM, Finegood DT (2000) A model of beta-cell mass, insulin, and glucose kinetics: pathways to diabetes. J Theor Biol 206: 605–619. [PubMed]
19. Maree AF, Kublik R, Finegood DT, Edelstein-Keshet L (2006) Modelling the onset of Type 1 diabetes: can impaired macrophage phagocytosis make the difference between health and disease? Philos Transact A Math Phys Eng Sci 364: 1267–1282. [PubMed]
20. Nelder JA, Mead R (1965) A simplex method for function minimization. Comput J 7: 308–313.
21. Rapoport MJ, Jaramillo A, Zipris D, Lazarus AH, Serreze DV, et al. (1993) Interleukin 4 reverses T cell proliferative unresponsiveness and prevents the onset of diabetes in nonobese diabetic mice. J Exp Med 178: 87–99. [PMC free article] [PubMed]
22. Mehta D (2005) Lysophosphatidylcholine: an enigmatic lysolipid. Am J Physiol Lung Cell Mol Physiol 289: L174–175. [PubMed]
23. Ott M, Gogvadze V, Orrenius S, Zhivotovsky B (2007) Mitochondria, oxidative stress and cell death. Apoptosis 12: 913–922. [PubMed]
24. Selivanov VA, Votyakova TV, Pivtoraiko VN, Zeak J, Sukhomlin T, et al. (2011) Reactive oxygen species production by forward and reverse electron fluxes in the mitochondrial respiratory chain. PLoS Comput Biol 7: e1001115. [PMC free article] [PubMed]

Articles from PLoS ONE are provided here courtesy of Public Library of Science