Invariant parameters for the experiments described in the paper:

A systematic approach to vaccine complexity using an automaton model of the cellular and humoral response

Brynja Kohler, Roberto Puzone, Philip E. Seiden and Franco Celada

Vaccine 19, 862-876 (2000)

 

Cells:

Upon receiving the necessary signals to divide, B, Th, and Tc cells divide for three time steps so that 8 new cells are produced. For B cells an additional division step provides 8 plasma cells.

Th0 cells stimulated by APCs become Th1 or Th2 cells with equal probability. The progeny of stimulated Th1 or Th2 cells retain their type, there is no Th1/Th2 conversion.

The three-in-a-bed stimulation of APC/Th/Tc does not require simultaneity. Up to 5 time steps are allowed between the APC/Th and APC/Tc contacts.

The initial population of APC, B, Th, and Tc cells is 1000 randomly distributed over the body array. The initial population of epithelial cells is 20,000. The B cells have complete receptor diversity but the Th and Tc cells go through a perfect thymus where all self-reacting cells and completely non-reacting cells are removed.

Plasma cells produce 5 antibodies on each time step.

B-cell anergy is not turned on.

B-cell mutation is not turned on. Since the B-cell diversity is complete mutation would not have an important effect.

It takes one unit of damage to activate an APC with 50% probability. Each necrotic cell death contributes one unit of damage to the site in which the cell death takes place. No contribution by antibody-antigen complexes is considered.

Half-lives:

Entity

Time steps

APC 50
B cell 50
Th cell 50
Tc cell 50
Epithelial cell 100
Plasma cell 10
Active state of an APC 50
Active state of a Th cell 50
Active state of a Tc cell 50
Anergy state 100
Antibody 10
Damage 3

 

 

 

 

 

 

 

 

Bound MHC/peptide complexes and antigen half-lives (other than the direct influence of the immune system) are 10,000 (essentially infinite).

Epitopes, peptides and binding :

Antigen epitope 51
Antigen peptide 240
MHC1 9
MHC2 234

 

 

 

These values are the decimal equivalent of the binary bit strings.  

There are no self-epitopes defined. The only self-peptide defined is 255 (but it plays no role in the experiments described in the paper).

All receptors have 8 bits, binding is allowed for 8 bits with a probability of 1 and 7 bits with a probability of 0.05. Matches of 6 bits or less do not bind. This holds for all specific receptor bindings except peptide/MHC binding. In this case the bindings for a 4, 3, 2, 1 bit match are 1, 0.5, 0.25 0.125, respectively. APCs bind antigen with a probability of 0.001. No bystander stimulation is allowed.

Presentation of peptides by B cells endocytosing antigen is on MHC2 only. Presentation by APCs is on MHC1 and MHC2 with equal probability. Presentation of peptides by infecting antigens (virus) is on MHC1 only.

General:

The size of the body array is 16x15 sites

The diffusion rate is chosen so that entities spread uniformly over the site and its six neighbors on one time step. Epithelial cells do not diffuse.

All runs are terminated at 2000 time steps (and called chronic) if they have not previously terminated due to cure or death. Cure results when all antigen is eliminated. Death results when 50% of the epithelial cells are infected or the viral load is greater than 200,000.