alexander2

The SBML for this model was obtained from the BioModels database (BioModels ID: BIOMD0000000290) Biomodels notes: Reproduction of the results in the first row of Fig. 2 of the original publication. The time course was integrated using Copasi 4.6. JWS Online curation: This model was curated by reproducing the figures as described in the BioModels Notes. No additional changes were made.

None

None

None

None

None

None

Self-tolerance and autoimmunity in a regulatory T cell model.

  • HK Alexander
  • Lindi M Wahl
Bull. Math. Biol. 2011; 73 (1): 33
Abstract
The class of immunosuppressive lymphocytes known as regulatory T cells (Tregs) has been identified as a key component in preventing autoimmune diseases. Although Tregs have been incorporated previously in mathematical models of autoimmunity, we take a novel approach which emphasizes the importance of professional antigen presenting cells (pAPCs). We examine three possible mechanisms of Treg action (each in isolation) through ordinary differential equation (ODE) models. The immune response against a particular autoantigen is suppressed both by Tregs specific for that antigen and by Tregs of arbitrary specificities, through their action on either maturing or already mature pAPCs or on autoreactive effector T cells. In this deterministic approach, we find that qualitative long-term behaviour is predicted by the basic reproductive ratio R(0) for each system. When R(0)<1, only the trivial equilibrium exists and is stable; when R(0)>1, this equilibrium loses its stability and a stable non-trivial equilibrium appears. We interpret the absence of self-damaging populations at the trivial equilibrium to imply a state of self-tolerance, and their presence at the non-trivial equilibrium to imply a state of chronic autoimmunity. Irrespective of mechanism, our model predicts that Tregs specific for the autoantigen in question play no role in the system's qualitative long-term behaviour, but have quantitative effects that could potentially reduce an autoimmune response to sub-clinical levels. Our results also suggest an important role for Tregs of arbitrary specificities in modulating the qualitative outcome. A stochastic treatment of the same model demonstrates that the probability of developing a chronic autoimmune response increases with the initial exposure to self antigen or autoreactive effector T cells. The three different mechanisms we consider, while leading to a number of similar predictions, also exhibit key differences in both transient dynamics (ODE approach) and the probability of chronic autoimmunity (stochastic approach).

Unit definitions have no effect on the numerical analysis of the model. It remains the responsibility of the modeler to ensure the internal numerical consistency of the model. If units are provided, however, the consistency of the model units will be checked.

Name Definition
1.0 item
86400.0 second
1.1574074074074073e-05 second^(-1.0)
1.1574074074074073e-05 second^(-1.0) item^(-1.0)
1.1574074074074073e-05 second^(-1.0) item
1.0 dimensionless
Id Name Spatial dimensions Size
body 3.0 1.0
Id Name Initial quantity Compartment Fixed
A 1.0 body
A_im 0.0 body
E 0.0 body
G 100000000.0 body
R 0.0 body

Initial assignments are expressions that are evaluated at time=0. It is not recommended to create initial assignments for all model entities. Restrict the use of initial assignments to cases where a value is expressed in terms of values or sizes of other model entities. Note that it is not permitted to have both an initial assignment and an assignment rule for a single model entity.

Definition
Id Name Objective coefficient Reaction Equation and Kinetic Law Flux bounds
r10 r10: A suppression by Tregs of other specificity A > ∅

b1 * A
r11 r11: A suppression by R A > ∅

sigma1 * A * R
r1a r1a: self-antigen uptake G > ∅

v_max / (k + G) * G
r1b r1b: pAPC maturation A_im > A

f * (v_max / (k + G)) * G
r2 r2: self-antigen release triggered by E ∅ > G

gamma * E
r3 r3: R activation by A ∅ > R

beta * A
r4 r4: R activation by A and E ∅ > R

pi1 * E * A
r5 r5: E generation by A ∅ > E

lambdaE * A
r6 r6: A death A > ∅

muA * A
r7 r7: R death R > ∅

muR * R
r8 r8: E death E > ∅

muE * E
r9 r9: G clearance G > ∅

muG * G

Global parameters

Id Value
R0 <assignment rule> dimensionless
b1 0.25 per_day
beta 200.0 per_day
f 0.0001 dimensionless
gamma 2000.0 per_day
k 50000000.0
lambdaE 1000.0 per_day
mA <assignment rule> per_day
mG <assignment rule> per_day
muA 0.25 per_day
muE 0.25 per_day
muG 5.0 per_day
muR 0.25 per_day
pi1 0.016 per_day_per_item
sigma1 3e-06 per_day_per_item
v_max 125000.0 items_per_day

Local parameters

Id Value Reaction

Assignment rules

Definition
R0 = f * (v_max / k) * lambdaE * gamma / (mG * mA * muE)
mG = muG + v_max / k
mA = b1 + muA

Rate rules

Definition

Algebraic rules

Definition
Trigger Assignments