 Research
 Open Access
 Published:
Modeling rejection immunity
Theoretical Biology and Medical Modelling volume 9, Article number: 18 (2012)
Abstract
Background
Transplantation is often the only way to treat a number of diseases leading to organ failure. To overcome rejection towards the transplanted organ (graft), immunosuppression therapies are used, which have considerable sideeffects and expose patients to opportunistic infections. The development of a model to complement the physician’s experience in specifying therapeutic regimens is therefore desirable. The present work proposes an Ordinary Differential Equations model accounting for immune cell proliferation in response to the sudden entry of graft antigens, through different activation mechanisms. The model considers the effect of a single immunosuppressive medication (e.g. cyclosporine), subject to firstorder linear kinetics and acting by modifying, in a saturable concentrationdependent fashion, the proliferation coefficient. The latter has been determined experimentally. All other model parameter values have been set so as to reproduce reported state variable timecourses, and to maintain consistency with one another and with the experimentally derived proliferation coefficient.
Results
The proposed model substantially simplifies the chain of events potentially leading to organ rejection. It is however able to simulate quantitatively the time course of graftrelated antigen and competent immunoreactive cell populations, showing the longterm alternative outcomes of rejection, tolerance or tolerance at a reduced functional tissue mass. In particular, the model shows that it may be difficult to attain tolerance at full tissue mass with acceptably low doses of a single immunosuppressant, in accord with clinical experience.
Conclusions
The introduced model is mathematically consistent with known physiology and can reproduce variations in immune status and allograft survival after transplantation. The model can be adapted to represent different therapeutic schemes and may offer useful indications for the optimization of therapy protocols in the transplanted patient.
Background
It is unfortunately not rare in medical practice that some diseases lead to organ failure, which may eventually require organ transplantation. The liver, the kidney and the heart are the most frequently transplanted organs. Diseases leading to organ transplantation span a wide spectrum of medical conditions: cancer, infections, autoimmune and degenerative diseases. Transplantation into the recipient of a foreign organ (graft), even from an individual of the same species (allograft), if left to itself causes rejection, a strong response by the recipient’s immune system leading to irreversible damage of the graft. Depending on the timeframe over which rejection occurs, “acute” rejections are differentiated from “chronic” ones. Acute rejection develops in the first few weeks or months after transplantation and is produced by cellular and molecular mechanisms, which may partially differ from those leading, over the course of many months or years, to chronic rejection.
After many decades of experimentation on animals and cells, and of the development of pharmacological tools, organ transplantation has evolved into a common therapeutic procedure. The success of a solid organ transplant relies in equal measure on the technical aspects of the implant and on the recipient’s acceptance or tolerance of the implanted graft. This last phenomenon is clinically induced by the administration of immunosuppressive drugs, which specifically decrease the recipient’s reactivity towards the graft, thus allowing the maintenance of the functional activity of the organ. Currently available drugs belong to several classes (calcineurine inhibitors (CNI), antimetabolites, target of rapamycin (TOR) inhibitors, steroids, and monoclonal antibodies)[1]. Canonical combinations of these drugs are typically used by the attending physician in a rather standardized fashion, attempting to maintain measurable drug plasma concentrations within established limits. Episodic and emergency use of immunosuppressive agents is then performed if signs of rejection become clinically evident.
Acute rejection response has been extensively studied in vivo and in vitro. From available studies, indications on the action of several drugs in controlling acute rejection and maintaining vitality of the graft have been obtained. However, even if the recipient initially accepts the graft, the more insidious and slow phenomenon of chronic rejection often ensues. The pathogenesis of chronic rejection is less well known and may depend not only on a partially different set of immune response mechanisms, but also on actual drug toxicity, recommending the use of the minimal clinically effective dose of medication.
Rejection has been widely described from a medical and biological viewpoint, but there have so far been no mathematical models describing this process. Mathematical models have been used to describe immunological behavior for a long time, beginning with the classical SIR model (SusceptibleInfectiousRecovered) first created for investigating the progress of an epidemic[2]. Several models describing the immune response in a number of diseases exist (HIV infection, tuberculosis, tumors..)[3–7] but, at present, only one model representing immune system dynamics during transplantation[8] has been published, which has the important aim to investigate Tcell population growth mechanisms, using thymus transplantation to follow the development of Tcells and their regulatory signals. The goal of the present work is somewhat different, in that we attempt to describe the main features of the immune system dynamics during general solid organ rejection. This will allow in the future the description of possible consequences of different therapeutic regimens. From the immunologic viewpoint, rejection mechanisms are substantially different from other immune responses (e.g. towards HIV, tuberculosis, or tumors), and the development of a specific model seems therefore warranted. Such a model should help transplantation clinicians and allied health care personnel in forecasting and treating rejection, without relying solely on empirical protocols. A good mathematical model should eventually allow the physician to consider in real time the several interrelated aspects of the immune response to transplantation, while jointly incorporating the known pharmacokinetics of the many potentially useful available drugs. The ultimate goal would be to help bridge the gap between the pharmacology and the biology of transplantation, explaining or at least representing the temporal relations between drug efficacy, possible drug adverse effects and the development of immune tolerance or graft acceptance.
In the present work we propose a tentative mathematical model of the rejection towards a solid organ transplant (kidney, liver, pancreas, heart). This model describes the evolution of the main cellular immune response as well as the kinetics and action of a single representative drug (e.g. cyclosporine). In order to have some physiological support for parameter assessment, a Mixed Lymphocyte Reaction (MLR) experiment was performed, on the basis of which the clonal expansion rate of Tcells could be determined. This experiment was chosen because the parameter of interest for the model was directly computable from the experimental data. Simulations with MatlabⒸ2010b have been performed and, on the basis of the experimentally determined clonal expansion rate and of relevant published material, the other model parameters were calibrated. While the current model is relatively simple, it introduces the main elements needed for the eventual description of more detailed response dynamics. Three case scenarios are discussed: nonimmunosuppressed transplantation and two immunosuppressive therapeutical regimens, one with moderate drug dose and the other with high drug dose. Additional simulations were made in order to explore different hypothetical therapy scenarios. Also, three state variables were selected as being most clinically relevant, and a sensitivity analysis was performed on these at two time points (one and ten years posttransplant).
Methods
Relevant physiology
After transplantation, the majority of exogenous molecules on the allograft are recognized by the immune system as selfantigens, as they are the same in both the donor and recipient. The main donor molecules which induce the immune response are those coded from the Major Histocompatibility Complex (MHC) genomic locus, which is the most polymorphic one on the human genome (it is essentially impossible to find two persons with the same gene cluster, except for monozygotic twins). Proteins coded by the MHC are located on the cell surface and bind antigen epitopes, creating a complex, which is recognized by the TCell Receptor (TCR) on Tlymphocytes, inducing the immune response. Cells which carry the MHC/epitope complexes are called Antigen Presenting Cells (APCs)[9].
In the conventional immune response, selfAPCs carry external antigen epitopes bound to selfMHC molecules. APCs normally circulate in the body and are found in every organ  during transplantation many APCs from the donor are introduced in the recipient, and they carry the donor’s MHCs. In fact, MHCs on donor cells are the major target of the rejection immune response[10], as Tcells can recognize complexes formed by allogenic MHCs[11]. This mechanism is called “direct activation”, since there is no processing of the antigen, and it is faster than the indirect activation, where recipient’s APCs have to process donor’s antigen[12]. In Figure1 a schematic representation of the direct and indirect mechanisms is shown.
Every different MHC on the surface of an APC can be recognized by and activate a different Tcell clone to proliferate (“clonal expansion”). It should be noticed that in the indirect mechanism many of the foreign epitopes involved are donor MHC fragments, due to the above mentioned polymorphism[13]. However, indirect activation can also occur in response to peptides derived from other molecules present in the allograft. “Alloreactive Tlymphocytes are [in any case] requisite mediators of allograft rejection”[11].
At the level of detail used in the present work, CD4+ and CD8+ Tlymphocytes (“helper” and “cytotoxic”) are not distinguished, and we consider them together as Tcells mediating rejection. The reasons for this simplification are that both Tcell types increase in number during the immune response and that they can both be activated by exogenous antigens. It has recently been shown that both class I and II MHC molecules (respectively recognized by CD4+ and CD8+ Tcells) can bind extracellular antigens (“cross priming”), while it was previously thought that class I MHCs could only present intracellular antigen[14]. In the course of direct activation, APCs are either destroyed or will eventually undergo apoptosis. Since donor’s APCs do not reproduce in the allograft, direct Tcell activation is a timelimited process. The immune response, however, does not terminate since there is continuous supply of donorspecific molecules, produced from the constantly proliferating cells in the graft. Tlymphocyte activation occurs as long as the graft is present in the recipient.
We can thus appreciate two main phases of the immune response to the graft: an initial strong response, sustained by the activation of Tcells due to a large but fading quantity of donor APCs, necessitating aggressive immune suppression therapy; and a later, more or less constant indirect activation of Tcells, sustained by continuously produced graft epitopes, for which less aggressive therapy is sufficient. It has in fact been reported that the number of directly activated Tcells is larger than the number of indirectlyactivated Tcells, the latter constituting less than 10 percent of the total cellular alloimmune repertoire[15].
Two aspects of graft rejection have not been explicitly included in the model for simplicity. The first is the increased production of lymphocytes from lymphoid organs, which receive various signals from stimulating molecules (cytokines). This mechanism has been shown[16] to take place when an inflammatory process (i.e. rejection) is ongoing. The second aspect of rejection is the appearance of the graftversushost disease, where donor’s Tcells present in the allograft react towards recipient’s antigens. This last phenomenon depends on the type of the transplanted organ and may be negligible in most cases of solid organ transplantation.
MLR is an in vitro experiment used to study alloreactive Tcells response to exogenous MHC molecules: it is used in clinical practice as a prediction rejection test before performing organ transplantation. MLR is induced growing mononuclear cells of an individual with those of another individual, these cells being isolated from peripheral blood: the difference between MHC loci of the two individuals induces clonal expansion of alloreactive lymphocytes[17]. One of the model parameters, the one corresponding to the clonal expansion rate of Tcells, was determined by an MLR experiment.
In the present model we consider one of the most frequently used immunosuppression therapy protocols, based on calcineurine inhibitors (e.g. cyclosporine)[1], which block Tcell clonal expansion. Drugs of this class inhibit signal transduction when TCRs recognize the epitope, so that the cell does not proliferate even when activated. Other types of therapy, acting through different mechanisms, could also be considered, extending the present approach.
The model
The model is represented in block diagram form in Figure2, showing all compartments and the relative dynamics. Organ transplantation is assumed to occur at time t_{ τ } in the life of the patient. Before transplantation, the major components of the immune system involved in transplant rejection are assumed to be at equilibrium. The choice of the time t_{0}, at which simulations begin, is therefore irrelevant, as long as t_{0} < t_{ τ }. In order to represent organ damage, we assume that antigen is released into the blood stream proportionally with the viable mass of the corresponding tissue, and consequently that the viable graft mass is proportional to bloodstream antigen mass, so that a substantial decrease in antigen concentration will indicate organ failure.
With the transplantation of an allograft at time t = t_{ τ }, the state of the immune system is suddenly and dramatically altered. The entry of a large amount of already processed foreign antigen, and the continuous production of foreign antigen by the functional graft tissue, induce a strong response by the host’s immune system. If this response is not controlled by immunosuppression, rejection and loss of the organ follow. Depending on the time and mechanisms needed to activate the immune response and on the duration of antigen permanence in the body, we describe different antigen types and assume that a specific Tcell type corresponds to each antigen type. In our model we use nonlinear ordinary differential equations to describe antigen and Tcell dynamics, using appropriate coefficients to specify the effect of antigens on corresponding Tcells and viceversa. State variables are generally defined in terms of concentrations, and we suppose a single vast plasma/interstitial fluid volume space where, given time, all species distribute. State variables and parameters, with corresponding units of measurement, are listed in Tables1 and2, respectively.
E is a generic environmental antigen: we suppose that the individual is constantly exposed to bacteria, viruses etc. and that the immune system is continuously stimulated by the corresponding exogenous antigens. When the allograft is transplanted in the recipient’s body, however, different antigen types are introduced. S represents antigen produced by graft cells: while these cells reproduce, this antigen is continuously formed. This antigen necessitates processing by endogenous APCs and presentation to Tcells. S in fact represents the nonprocessed form, while L is the processed form of the same antigen. The distinction between the processed and the nonprocessed forms explains the delay in processed antigen effect with respect to the directly presented antigen. U indicates alloantigen directly presented by donor’s APCs, which does not form anew in the organ. This antigen is ready to be presented to Tcells, so its effect in activating Tlymphocytes is rapid and strong, but vanishes as APCs are progressively cleared from the organ. Alloantigens which are indirectly presented and which derive from cells not reproducing in the allograft, have to be processed by endogenous APCs. As it was discussed in the previous section, this type of antigen does not reform in the recipient once it is eliminated: Y represents the nonprocessed form, while Z represents its processed form, ready to activate Tcells. The processing mechanism is represented in the same way as for antigens S and L. The model is detailed in the following equations:
In equation 1, describing antigen E dynamics, the rate k_{ e } represents the entry of environmental antigens into the body, assumed to be constant throughout the considered period of time (before as well as after transplantation). The elimination terms k_{ xec } and k_{ xe } describe antigen neutralization due respectively to Tcell action and to Tcellindependent elimination of the antigen (as it happens e.g. through chemical and physical elimination mechanisms, such as lipases, mucus secreted by respiratory and gastrointestinal tracts etc.).
Equation 2 represents the dynamics of Tlymphocytes which react towards antigen E: k_{ c } indicates constant physiological Tlymphocyte production from lymphoid organs, k_{ cac }(F)E C_{ e } represents Tcell clonal expansion after antigen contact, which is inhibited by drug action, by setting the rate k_{ cacF } as a proper function of the drug concentration F. Tlymphocytes are also “consumed” by antigen, in the sense that upon Tcell interaction with antigen, the lymphocyte eventually undergoes apoptosis (programmed cell death): this is described by the elimination term k_{ xce }E C_{ e }. We introduced another elimination term, k_{ xc }C_{ e }: lymphocytes die for apoptosis even if they do not encounter any antigen after a certain period, and we assume this mechanism to be proportional to Tcell concentration.
The regenerating antigen is described by equation 3: at time t_{ τ } there is an impulsive entry, modeled by a Dirac delta (Eq. 12). Once the organ is transplanted, its cells regenerate: the growth rate is assumed to be logistic of parameter k_{ s }, limited by a carrying capacity S^{∗}, so that S concentration tends towards S^{∗}, whether it is above or below it. The elimination term k_{ xsc }C_{ l }S, depends on S concentration and on Tcells primed from the processed form of the antigen. The nonprocessed form is not ready to activate Tlymphocytes, but it is destroyed by Tcells activated from the processed form. In fact, Tcells activated from the processed antigen (C_{ l }) are primed to react against cells carrying the same epitopes (graft cells). Organ rejection is thus represented by the S antigen elimination, this being related to organ mass.
The processed form of the regenerating antigen is described by equation 4, where the only positive entry is represented by the term k_{ ls }S depending on the unprocessed antigen concentration; the two elimination terms are similar to those in the previous described equations, representing lymphocytedependent and independent elimination, respectively. Equation 5 describes the C_{ l } lymphocyte dynamics, with the constant entry, the clonal expansion depending on the processed antigen L, elimination upon encounter with both antigens S and L, and antigenindependent elimination. The directly presented antigen U is represented by equation 6: a Dirac’s delta describes impulsive entry at transplantation time, while elimination happens in two ways, dependent and independent from Tcells, respectively. Uspecific Tlymphocytes, C_{ u }, are represented in equation 7, including constant production rate, clonal expansion, antigendependent and independent elimination.
The last three equations represent the nonregenerating antigen which has to be processed, its nonprocessed form (Y ), the processed form (Z) and Zspecific Tcells. The dynamics are similar to the S, L and C_{ l } subsystem, the only difference being that this antigen does not regenerate. Equation 8 has an impulsive entry and elimination dependent and independent from Tlymphocytes; in equation 9 the entry depends on Y concentration, and the usual two elimination terms follow; equation 10 is similar to equation 5. Equations 12, 13 and 14 describe the impulsive antigen entries into the system. In each case an antigen concentration (respectively S_{ τ }, U_{ τ } and Y_{ τ }) is multiplied by a Dirac delta term acting at time t_{ τ }. The end result is the representation of the appropriate instantaneous change in S, U and Y compartments at the time of transplantation.
Finally, immunosuppressor pharmacokinetics is described by equation 11. For the purpose of the present model, given the long timescale considered, drug administration is assumed to be continuous, with average rate k_{ f } (equation 15) of delivery into the circulation. The drug is eliminated from the circulation following a linear, firstorder process with rate constant k_{ xf }. Since before transplantation no treatment is administered, k_{ f } is 0 before time t_{ τ }, while it is equal to κ_{ f } from the time at which therapy begins, which we assume to be t_{ τ }. The term\frac{{\kappa}_{f}}{{k}_{\mathrm{xf}}}\delta (t{t}_{\tau}) represents therefore the (impulsive) loading dose of the drug, necessary to bring it instantaneously to the equilibrium level, at which it is constant thereafter. After antigen contact, Tcells would spontaneously give rise to a (fast) clonal expansion. The immunosuppressive effect of the drug, leading to a slower increase of Tlymphocyte concentrations, is described by the exponential decrease in the clonal expansion coefficient k_{ cacF } produced by proportionally increasing drug concentrations F (with effect rate constant λ): this is represented by equation 16. The drug is in fact thought to block signal transduction after antigen contact, in a concentrationdependent fashion, thus disabling clonal expansion and reducing Tcell proliferation. Equation 11 represents the pharmacokinetics of the antirejection drug F with given constant entry (depending on the administration scheme), linear elimination and with impulsive entry assumed to be simultaneous with transplantation. Equations 1216 and Tables1 and2 define each symbol used. Tables1 and2 also report units of measurement for all variables and parameters.
Model parameters
The model has not been fitted to experimental data, although an MLR experiment has been performed to assess the maximal lymphocyte clonal expansion rate and compute the parameter k_{ cacF }. Indicative parameter values for the processes modelled are difficult to find in the literature. For this reason, the main criterion followed for parameter calibration (other than k_{ cacF }) was the production of relative timecourses of the relevant state variables, which appeared consistent with physiology and clinical experience to the medical doctors among the authors, while remaining within a broadly acceptable range of magnitude. For instance, k_{ xc } was set at 1/year, whereas values around 0.1/year can be computed indirectly from normal lymphocyte apoptosis results in vitro[18].
Mixed lymphocyte reaction
Blood samples were collected from four healthy donors and peripheral blood mononuclear cells (PBMCs), from three of the blood samples, were stained with the lipophilic fluorescent molecule Carboxyfluorescein Succinimidyl ester CFSE. The CellTrace CFSE Cell Proliferation Kit (Invitrogen) was used following the protocol provided by the manufacturer. Cells were stained with PBS/5%FCS/CFSE50 μ M for ten minutes at room temperature and then washed twice. Nonlabeled PBMCs obtained from the fourth healthy donor were used as allogenic stimuli and labeled PBMC were cultured in vitro in the presence or absence of nonlabeled PBMC with a ratio of 2:1 (2×107 labeled PBMC versus 107 nonlabeled PBMC). Three days later, cells were harvested and stained with PEconjugated antiCD3 mAb (Becton Dickinson). Cells were then analyzed at cytofluorimeter FACScalibur (Becton Dickinson) and data were acquired by the software CellQuest pro.
Parameter computation
Since lymphocytes were exposed to a large amount of allogenic cells, we assume that the experiment reflects the maximal clonal expansion that can be achieved at a “maximal” antigen concentration. From equations 2, 5, 7 and 10, k_{ cacF } dimensions are l/μ g/year. From the experiment we measured the percentage of replicating cells per day, which has to be divided by the maximal antigen concentration expressed as μ g/l. The concentration of MHC molecules in the experimental preparation was approximated as follows. On a cell surface there are approximately 10^{5} MHC molecules, and there were approximately 10^{6} cells/ml of blood. It thus follows that the concentration of MHC molecules was 10^{5}×10^{6}/ml, i.e. 10^{14}molecules per liter. Considering both MHC class I and class II molecules, the average molecular weight is approximately 60 kDa, meaning that one mole of MHC weights 60×10^{3}grams. Since in one mole there is one Number of Avogadro of molecules, 10^{14} molecules correspond to 10^{14}/(6×10^{23}) moles, or approximately (10^{18}×6)/(10^{23}×6)) g/l. The parameter k_{ cacF } has therefore been computed as the clonal expansion per day divided by the maximal antigen concentration (10^{−5}g/l or 10 μ g/l), afterwards multiplied by 365 days. No formal statistical parameter estimation (e.g. Maximum Likelihoodbased) was attempted. In keeping with the generally qualitative character of this model’s predictions, the desired outcome of the MLR experiment was an indicative, plausible value for the rate of clonal expansion, and, as is apparent, this plausible value itself is conditional on adhoc assumptions (e.g. that maximal stimulation is equivalent to an antigen concentration of 10 μ g/l).
Results
Clonal expansion rate from MLR
In Figure3 a picture from the cytofluorimetric analysis is shown. Cells labeled with antiCD3+ are Tlymphocytes, and those who have replicated have a lower content of CFSE, thus the upperleft panel shows Tcells that have reacted to the exogenous antigen. In the two samples that replicated we had respectively 23.9 and 23.1% of cells that had divided, that is an average of 23.5%, which means that 10023.5=76.5% of the cells did not replicate. Cells went through 3 cell cycles in 3 days, so it can be assumed that cells went through one cycle per day on average. At each cycle, cells which divide double, so the amount of cells that are alloreactive and start replicating at day 1 is 23.5%/2^{3}≃3% of the final number. So, for every 100 cells in the final count, on day 1 we have 76.5 (cells that did not replicate) + 3 (cells that replicated), that is 79.5% of the final 100%. The global replication rate at each rate was computed solving the following equation:
where y_{0} is the total number of cells on day 0 (approximately 79) and y_{3} is the number of total cells on day 3 (100), while k_{ max } is the replication rate. Solving (17) with t = 3 we obtain k_{ max } = 0.08 /day. k_{ cacF } = 0.08/10^{−5}l/g/day, which is 2.92 l/μg/year.
Model simulation
The model has been implemented in MatlabⒸ2010b, and simulations are presented showing the behavior of the several types of antigen and corresponding specific Tcell populations after transplantation of a solid organ. Three scenarios are depicted, corresponding respectively to the notherapy, moderate therapy and maximal therapy situations. A time range from 35 to 60 years is shown, hypothesizing that transplantation occurs at time 40 years.
In Figure4 all antigen types, the corresponding Tcell dynamics, and drug concentrations, are shown in the three therapy cases; subfigure4.1 shows drug dynamics. In all subfigures the solid line (–) represents the notherapy case, the dashed line ( ) refers to the moderate drug dose while the dotted line (.) refers to the high drug dose. Variable concentration ranges, as will be explained in the discussion, are in agreement with physiological limits.
Concentrations are assumed to be constant if no traumatic events happen during life. As the allograft is introduced (taking, e.g., t_{ τ } = 40 years), all antigen types described by the model (with the exception of the environmental antigen E) go from 0 to a high level. In response, specific Tcell concentration also grows from a normal, low level to a high one. The solid lines correspond to the notherapy case: at the beginning antigen concentrations are very high and so Tcells are quickly activated, causing antigen elimination, which only partially regenerates (case of antigen S and its processed form, L). After this first period, as Tcells react towards the graft to destroy it, antigen concentration (which is proportional to graft mass) decreases. Following antigen reduction, Tcell concentrations are also reduced, but remain consistently higher than baseline (i.e. before transplantation), as antigen from the organ is continuously produced and never vanishes (Tlymphocytes C_{ l }, Figure4.6). If no therapy is applied, Tcell levels remain high and will bring the organ to a minimal size (rejection and failure): as can be noticed, the continuous line in the regenerating antigen graph (concentration of antigens S and L, Figure4.4 and 4.5) is very low.
The dashed line represents the case in which administration of a moderate amount of immunosuppressive drug (such as cyclosporine) takes place. We assume that the drug is given simultaneously with the organ transplantation and that the patient is continuously and constantly treated (subfigure4.1). It should be noticed that, as the increase of Tcell concentration is much lower, antigen level remains higher than the notherapy case. This indicates that the graft is not totally destroyed by the immune system.
If therapy is much stronger (dotted line), the graft normally survives. The problem is that, with immunosuppression, Tcell levels fall well below normal. As Tcells are less aggressively attacking the allograft, its regeneration allows the attainment of a constant, sizable equilibrium level, but Tcells do not proliferate as much and are thereby less effective not only towards the graft’s antigen, but also towards environmental antigens, exposing the patient to opportunistic infections. In fact, while rejection is prevented, patients might not be able to defend themselves from severe general infections. Any chosen intensity of therapy represents a compromise between desirable graft tolerance (with attending functional organ size) and dangerous lowering of general immune defenses.
SubFigures4.2 to4.11 show the time course of specific antigens and their respective Tcell dynamics. In subFigures4.2 and4.3, generic antigen and generic Tlymphocyte concentrations are reported. If no therapy is administered, these dynamics are at equilibrium. With therapy, generic Tcell concentration decreases (depending on drug dose): this indicates that immunosuppression is not specific in lowering Tcell expansion towards the graft. Instead, it reduces proliferation of all Tcells and, as a consequence, environmental antigen permanence levels in the body (E) increase.
In subfigures4.4,4.5 and4.6, the concentrations of regenerating antigen in the nonprocessed (S) and processed (L) forms, as well as Tcells that respond to it (C_{ l }), are shown. As seen before, in the absence of therapy, low antigen and high Tcell levels are reached after transplantation, while with drug administration antigen and Tcell concentrations are respectively higher and lower. This is the only case in which antigen never goes to zero because it is produced from regenerating graft cells. Dynamics of directly presented, nonregenerating antigen (U) and corresponding Tcells (C_{ u }) are shown in subfigures4.7 and4.8: without therapy Tcells rapidly expand and consequently the antigen is rapidly eliminated, while, with therapy, the action of Tcells is less aggressive on this type of antigen, which will eventually be eliminated because it does not regenerate. The nonregenerating and indirectly presented antigen (subFigures4.9 and4.10) has to be processed: there is a delay in the increase of Z, which depends on nonprocessed antigen (Y ) dynamics. The unprocessed antigen rapidly grows and rapidly vanishes as it is processed to the indirectly presented form and is no longer produced (because it derives from APCs not reproducing in the graft).
In subfigure4.11 the dynamics of specific Tcells primed for antigen Y and Z is shown. In this case as well, as drug dosages increase, Tcell concentrations decrease and antigen levels rise.
From the graphs shown it is evident that is not easy to find the right drug dosage. In fact, it can be seen that when the dose is moderate Tcell levels are adequate but the graft tissue size is too small, while, as drug concentrations increase, Tlymphocyte levels are not sufficient to defend the patient from infections. The perfect situation would be to find a therapy level, which is effective in saving the graft from rejection, but does not unduly expose the individual to infections.
Effect of immunosuppression on Tcell reaction towards infection
In order to explore the predictive ability of the model, a simulation of an infection occurring two years after transplantation has been performed: the environmental antigen E and its specific Tcell population are shown in Figure5. It is clear that, if the patient is not under treatment with the immunopuppressor (solid line) the immune system reacts normally (big increase in Tlymphocyte concentration) and the antigen is quickly brought back to equilibrium levels. If the subject is immunosuppressed, Tcell population growth is inhibited in a drugdosagedependent fashion, and the patient is not able to fight the infection (elevated E antigen levels maintained for a long time).
Balance between organ survival and immunosuppression
As will be discussed in the Appendix, the subsystem S, L and C_{ l } has two equilibrium points. One is provided by the elementary triple S = 0, L = 0, C_{ l } = k_{ c } / k_{ xc }, corresponding to a situation of no antigen S, which may refer to the pretransplantation case, when the patient has no actual exogenous organ mass. According to a proper setting of the model parameters (see the Appendix for more details), this equilibrium point is shown to be stable with respect to perturbations of L and C_{ l }, and to be unstable with respect to perturbations of S. On the other hand, the other, asymptotically stable equilibrium point corresponds to a nonelementary solution for the three state variables and is eventually reached after transplantation, when antigen and Tcells dynamics, under a certain drug dose, are balanced. Let us denote by\stackrel{\u0304}{S} the value of S at the stable equilibrium that is eventually reached after transplantation, which represents the carrying capacity for antigen concentration (which we assume proportional to organ mass), as explained in subsection 3.2. A simulation is performed to show how\stackrel{\u0304}{S} changes in relation to varying immunosuppressant concentrations. The plot in Figure6 shows how the value of\stackrel{\u0304}{S} varies with increasing drug doses: at F equal 0,\stackrel{\u0304}{S} has a positive (non zero) value, which is physiologically plausible as the antigen would not be completely eliminated without drug, but its concentration would be low. As drug concentrations increase,\stackrel{\u0304}{S} increases following a saturation curve: at high drug levels, when lymphocytes are inhibited, the antigen equilibrium approximates its maximum level. Similar diagrams are reported in Figures7,8,9,10 (see Appendix).
Hypothetical cyclical therapy
An additional simulation was performed to test the effect of a hypothetical “intermittent” therapy. The aim was to investigate the physiological response in the case of a cyclical administration of immunosuppressant. The period of interruption equals the period of treatment, and the given dose is the same in each treatment period. A comparison between continuous and cyclical therapies is shown in Figure11. The time interval applied in the simulation shown is one year; several other time intervals were tested (1 week, 2, 3, 4, 6 months, 2 years) and results were similar (results not shown). As shown in Figure11.1, the immunosuppressant is given at one year intervals. The red line represents the continuous therapy (low dose), while the blue and green lines the intermittent therapies (low and high doses). In Figures11.2 and11.3 the effect of the three therapy schemes on the environmental antigen (E) and the respective Tcells (C_{ e }) are shown. Figures11.4,11.5 and11.6 show the dynamics of antigen S and its processed form L  representing the organ mass  and the corresponding Tlymphocytes (C_{ l }). The low intermittent dosage is less aggressive towards Tlymphocytes and, correspondingly, the subject is better protected against environmental antigens, but the organ is not protected from rejection as much as with the continuous dosage. The high intermittent dosage is less effective than the low continuous one, Tcells are not sufficiently inhibited and the graft antigen level is low.
Sensitivity analysis
Three variables have been identified as the most important ones to describe rejection and immunosuppression in a clinical setting: S and T, representing the organ mass (nonprocessed and processed form of the organ antigen) and C_{ e }, the Tcells specific for the environmental antigen. Sensitivity analysis of these variables at 1 year and at 10 years after transplantation (t = 41 and t = 50) has been performed (tornado plots are shown in Figure12).
C_{ e } at t = 41 and Ce at t = 50 plots are reported in Figures12.1 and12.2, respectively. The first observation is that the parameters affecting the two targets are the same, but in a different priority order. Also, at one year after transplantation, the variable C_{ e } is generally more sensitive to parameter variations compared to the (near) equilibrium state at 10 years posttransplantation. The k_{ xc } parameter is most influent on both variables: it represents physiological Tcell elimination, due to natural apoptosis. It is noticeable that the other elimination rate, k_{ xce }, which is the elimination due to Tcell interaction with the antigen, is much lower. This is understandable, since k_{ xc } represents a general mechanism that involves all cells while k_{ xce } only applies to cells in contact with their specific antigen. Another aspect worth noticing is that the rates which directly influence the environmental antigen (E), k_{ xe }, k_{ xec } and k_{ e }, are important at t = 41 but negligible at t = 50 (variation 0.05% or less).
Figures12.3 and12.4 show tornado plots for S at t = 41 and S at t = 50. The parameters, to which the two targets for the S antigen are most sensitive, are broadly the same. The major difference is in antigen regeneration. For both t = 41 and t = 50 variations in S^{∗} are relevant, whereas variations in k_{ s } only impact modelpredicted S levels at 10 years posttransplantation, likely due to an accumulated effect of the small rate change. Besides the above described parameters, which regulate the S regeneration term, the parameters to which the variable S is most sensitive are the ones which regulate drug action on Tcells (k_{ xf } inverse correlation with S increase, and lambda direct correlation with S increase) and the Tcell clonal expansion rate (k_{ cacF }). Also the transfer rate from S to L antigen forms (k_{ ls }) is a relevant parameter. It must be kept in mind that S is the nonprocessed form of the regenerating antigen, and it does not directly activate Tlymphocytes, while it is directly eliminated by them. In fact, the parameters governing the variation of L, the antigen processed form, are somewhat influent on both targets.
Tornado plots for L at t = 41 and L at t = 50 are reported in Figures12.5 and12.6, respectively. The processed form of the regenerating antigen is mostly affected by variations of the parameters lambda, k_{ cacF } and k_{ xf }, as described above for the nonprocessed form, while it is not much affected by those parameters which regulate the processed form of S (S^{∗} and k_{ s }), both at time 41 and 50. One interesting observation concerns the sensitivity to the parameter k_{ ls }, the transfer rate from the S to the L form: one would expect an increase in this parameter to result in an increase in the target, while at t = 41 there is 1% variation in L. This is probably due to the fact that, in the moment in which the Tcell expansion stimulation is most effective, an increase in processing rate leads to an increase in the form available to activate the cells, which results in a faster clonal expansion and thus lower antigen level. At time 50, instead, the effect is almost negligible, because we are at an equilibrium situation.
As a general comment to the sensitivity analysis it is evident that the parameter k_{ cacF } is surely relevant. Accurately assessing the value of the k_{ cacF } parameter (derived from the MLR experiment) seems important for a correct quantitative prediction of the time course of the lymphocyte populations, as could naturally be expected. This point should be kept in mind upon applying the model in a clinical context, possibly predicting an individual patient’s posttransplantation course.
Discussion
Mathematical models are increasingly used in biology and clinical medicine to express concise, mechanistic descriptions of ongoing phenomena. The possibility of representing a pathophysiological process by means of a mathematical model allows the investigator to formalize beliefs, compare interpretations and simulate hypothetical scenarios of interest. Mathematical models have in particular already been introduced into several areas of immunology[3–8]. So far, however, no mathematical model has yet been presented describing allograft rejection in order to support the evaluation of therapies.
The clinical problem, which characterizes the management of the transplanted patient, is the difficult adjustment of immunosuppressive therapy, walking the fine line between undersuppression, with ensuing organ rejection, and oversuppression, with the danger of potentially lethal opportunistic infections. While a wide spectrum of active pharmacological agents are now available to the transplantation specialist, their mechanism of action is often incompletely understood and their precise effect on the complex balance of immune system competence is not quantitatively determined. Therapy therefore follows ruleofthumb principles, intensive monitoring of potential damage indicators (like serum creatinine for kidney, hepatic enzymes for liver transplantation), meticulous monitoring of plasma drug levels. The relationship between the time courses of drug effect, Tcell cycle and organ damage is however a matter of guesswork, only partially mitigated by the relatively precise knowledge of the pharmacokinetics of the immunosuppressive drugs themselves. In fact, what pharmacological information is currently offered to clinicians consists largely of singledrug pharmacokinetics parameters. Recent experiences have indeed suggested that the possibility of studying drug pharmacokinetics and pharmacodynamics through modeling techniques (in silico) may greatly reduce the need for animal and cellular models[19], as well as the discomfort and risks associated with extensive human experimentation. In order to work in concrete, however, the modeling approach requires a tight interconnection of mathematical constructs and physiological knowledge.
The model presented here describes established physiology mechanisms, whose outcomes, however, are not directly detectable with clinical measurements. Specific Tlymphocyte clonal expansion after foreign antigen contact, which is known to be the first step of the cascade leading to organ rejection[11], is not directly measured in clinical practice[20].
Variables concentration ranges, reported in the plots, were obtained setting the parameters in order to have a physiological reproduction of rejection, with the exception of the experimentally derived parameter k_{ cacF } which was maintained constant. It is worth noticing that the concentration ranges reflect physiological values. Regarding Tcell concentrations, it has been reported that the physiological range for Tcells (CD4+ and CD8+) is 7503600 cells/μ l[17], which is of course variable depending on the subject and on the situation (immune status, period of the year, general physical conditions etc.). Experiments reported in the literature give ranges of approximately 12001700 cells/μ l[21] in physiological conditions. In the plots, if we sum the concentrations of all antigen types before transplantation, we have a value of about 2∗10^{3}cells/μ l. As regards drug concentrations, the reported values for cyclosporine treatment are 150350 ng/ml, depending on the period following transplantation[22] which is also in agreement with our simulation: in the plot, the hypothesized drug concentrations are approximately 100 (low dose) and 350 (high dose) ng/ml. Regarding antigen concentrations, the situation is more complicated, since values for it cannot be found (as far as we know) in the literature. However, while still remaining consistent with cellsurface molecule densities as they are generally known, we may simply assume arbitrary concentrations proportional to organ mass: a rapid decrease in this arbitrary antigen concentration would signal an ongoing rejection process.
The proposed model describes the pretransplantation equilibrium state, characterized by constant environmental antigen and Tcells level, which is dramatically perturbed by the entry of a large amount of alloantigens. Regenerating antigen determines the continuation over time of the immune system activation, leading to possible chronic rejection. The administration of therapy limits the immunological aggression towards the organ and the natural ability of the allograft to reproduce makes it so that an equilibrium is attained at a nonzero level of remaining allograft tissue. This is potentially the most useful area of application of future versions of the present model, which will incorporate, besides a general biological description of the immune response, also a precise quantification of the applicable pharmacokinetics (possibly depending on the patient or on patient subgroups).
Immunosuppressive therapy is very invasive and the substantial risks of potentially severe side effects have been widely discussed[23, 24]. Models for therapy improvement (e.g. drug dosage) have been proposed, so far only considering single aspects of the therapy (e.g. plasma drug concentrations) or focusing on the action of a single specific drug[25, 26]. The model presented here has instead the aim of framing drug kinetics and effects within a simplified representation of the relevant immune system biology. The somewhat empirical therapy adjustments in clinical practice, which at present are based on organ function damage indicators and drug level monitoring, may therefore be complemented, using a model similar to the one presented here, by a quantitative systemic assessment of the likely impact of posology alterations, considering therapy effects in the context of patient individual characteristics and immune system status.
The study of the present model prompts, in fact, some interesting considerations.
One aspect worth noticing, which in clinical practice is subject to iterative attempts, is therapy adjustment. The model explicitly shows that a constant dosage of one immunosuppressant is never satisfactory since no good compromise can be achieved in this way between organ survival and acceptable patient immune defenses. Even if this fact is widely appreciated among clinicians (and in fact therapy is adjusted testing it directly on the patient), there has never been, in our knowledge, a direct demonstration of it[22]. With the present model we have attempted to follow in detail the fate of several among the most meaningful cellular and chemical species involved in the immune response to organ transplantation. In so doing, we attempted a mechanistic description of those factors promoting species accumulation and decay, thereby falling naturally into the framework of mass action kinetics.
In the present work, a simulation where intermittent therapies were tested (see Figure11) indicates that this kind of intermittent treatment would not be effective. In fact, even if the intermittent dosage is almost three times the continuous one, therapy is more effective with the latter.
In the model, the dynamics of the several types of antigen and lymphocytes are coupled through the effects of immunosuppressive therapy, acting equally on all types of Tcells. Decoupling of the several state variables, linked only by the drug effect, is actually a relevant feature of the model: therapy acts on the whole immune system; if it could act only on those cells which specifically clear graft antigens, immunosuppression would not affect individual protection against infections. This aspect also underlies how important it would be to have more specific data from transplanted patients. Also, data availability would definitely be useful for parameter identification.
Once a robust biological model is in place, it becomes relatively easy to incorporate the effect of different drugs. It will therefore be possible to express the suppression of clonal expansion, a greater mortality of Tlymphocytes, or even a generalized action in suppressing the inflammatory response (as may happen when administering corticosteroids). The problem here will not be as much in introducing the specific actions of the array of available therapy schemes, commonly used in clinical practice, but rather in representing with some degree of accuracy those side effects, which make it undesirable to simply increase without bounds the dosage of immediately useful agents. It will become possible, in this way, to support the decisionmaking of the attending physician or surgeon, who has to choose a reasonable compromise between immediate therapeutic effect and longterm complications.
The model presented in this work has been developed with the aim of allowing the eventual representation of different mechanisms of action, hence of the effects, of different classes of immunosuppressive drugs. Mechanisms can differ either from the molecular or the cellular viewpoint. There are different steps, along the pathway of Tlymphocyte activation, at which drugs can act (resting state, early activation, late activation and proliferation). Polyclonal antilymphocyte antibodies act at the resting state. Calcineurine inhibitors (cyclosporine, tacrolimus) act in the early activation pathway, so they have the same inhibition mechanism from the cellular viewpoint. However, cyclosporine and tacrolimus have different chemical structure, and act with different mechanisms at the molecular level. It has in fact been reported that tacrolimus is more effective than cyclosporine, it is used in smaller concentrations, and there are differences in their side effects[27, 28]. Monoclonal antibodies and rapamycin (TOR) inhibitors act in the late activation step. Antiproliferative drugs (azathioprine and mycophenolate acid) act on the last step of the activation pathway. Corticosteroids have a very different mechanism of action in that they do not inhibit Tcell production, but they act nonspecifically on the inflammatory process, preventing organ failure without directly acting on Tcell dynamics[1]. The current model can also be modified by explicitly representing different steps of the activation pathway as well as focalizing on molecular aspects for a better description of different mechanisms of action.
In this representation, some simplifications have been deliberately introduced. Among these, no discrimination has been made concerning the different Tcell types: Tlymphocytes can be either naive or activated; once activated, they differentiate into cells with specific roles (mainly helper and cytotoxic); in the present model, however, the global class of Tcells is considered, representing the response to transplantation of the immune system as a whole. The consideration of different cellular types, besides TLymphocytes of the CD4 and CD8 classes, would in fact be helpful in refining the description of the chain of events involved in the inflammatory response. Other cells of the immune system (e.g. antigen presenting cells, Blymphocytes, macrophages etc.) as well as cytokines (responsible of cell proliferation, signaling and recruiting inflammation agents) are also involved in the rejection mechanisms. In the present work the need to limit model complexity has prompted the decision of representing only the cellular compartment (Tlymphocytes), most representative of solid organ rejection reaction, which directly increases in response to incoming antigen and triggers the rejection response. While the inflammatory response is not followed in detail, a measure of the inflammatory damage to the organ is however represented by the amount of circulating L antigen, assumed to be proportional to the tissue mass of living allograft.
Another simplification consists in not representing explicitly the increased overall Tcell production occurring in the presence of inflammation. When an inflammatory process is ongoing (e.g. during rejection) lymphoid organs are stimulated to nonspecifically increase cell production. These mechanisms are poorly understood and the actual increase in competent Tcells may not be so high as to substantially modify the response: for this reason, a constant Tcell production was assumed (k_{ c }).
The main limitations of the current model invest both the detail of therapeutic manipulations it describes and the plausibility of the represented biology. While the model, as discussed above, can be easily and naturally extended to account for more than the single pharmaceutical agent (F) incorporated so far, there are in fact important aspects of the immunological response to transplantation which have not yet been tackled. One such is the description of the Graft versus Host Response, which is of great importance in explaining the events following transplantation of lymphoid tissue (like bone marrow transplants), particularly after massive immunosuppression of the recipient before the operation: for this reason, the present model should be considered appropriate only for solid organ transplants (liver, kidney, pancreas, heart). Another area where greater biological detail would be useful is that of the description of the chain of events in the inflammatory process which underlie the clinical features of chronic rejection. While within the framework of the present model no distinction has been made between acute (or indeed hyperacute) and chronic rejection mechanisms, factors leading to the different types of rejection may be different, and may be the object of one type of model refinement. In particular, while it is well known that acute rejection is mediated by CD4+ and CD8+ Tlymphocytes stimulated from exogenous MHC, mechanisms leading to chronic rejection are still not completely clear, and the latter is now the most common reason of graft loss from the recipient.
While the model as reported does offer useful insights in the reciprocal variations of antigen and immune cell species during a generic, hypothetical solid organ transplantation, model parameter estimation has not been carried out and no quantitative prediction can be strictly constructed, not to mention the assessment of prediction uncertainty. Indicative parameter values for the processes modeled are difficult to find in the literature. For this reason, the main criterion followed for parameter calibration was the production of relative timecourses of relevant state variables, which appeared consistent with clinical experience to the medical doctors among the authors, while remaining within a broadly acceptable range of magnitude. A priori identifiability analysis of the model has not been performed, and no data fitting has allowed us to assess a posteriori regions of confidence on parameter values. As a consequence, the model identifiability issue remains completely open.
Conclusion
The graphs, reporting the time course of the different types of antigen and corresponding Tcell populations, show that the model captures well the clinically expected behavior of the transplanted organ mass and of the immune system reaction, under the three scenarios of no therapy, moderate therapy and aggressive therapy. From the graphs it is evident that is not easy to find the right drug dosage. It can be seen that when the drug dose is moderate Tcell levels are adequate to prevent opportunistic infections, but the Lantigen level, corresponding to the viable graft tissue, is rather low. Conversely, at a drug concentration effective in maintaining the entire transplanted tissue mass, the Tcell population is suppressed excessively and the risk of complications would appear to become substantial. The model therefore predicts that single drug therapy is likely to be inadequate to safely prevent graft rejection, in accord with the clinical experience so far accumulated. A perfect situation might not exist, but the theoretical exploration of drug combinations and of nonconstant therapy schemes could be one way to obtain useful indications for the biological experimentation of novel therapeutic protocols.
The present work proposes then a first mathematical model of the cellular immune response to solid organ transplantation, addressing both acute and chronic rejection. The model’s mathematical behavior is broadly consistent with known physiology and longterm variations in immune status and allograft survival. The model can be tailored to address specific organ transplantation situations, may be naturally adapted to the representation of different therapeutic regimens and may offer useful indications for the optimization of therapy protocols in the transplanted patient.
Appendix: Qualitative behavior of the solutions
From a mathematical point of view, the whole system (110) may be split into the four independent subsystems composing it, namely:

the environmental antigen subsystem:
\begin{array}{c}\frac{\mathrm{dE}}{\mathrm{dt}}={k}_{e}{k}_{\mathrm{xec}}{C}_{e}E{k}_{\mathrm{xe}}E\hfill \\ \frac{d{C}_{e}}{\mathrm{dt}}={k}_{c}+\left(\right.{k}_{\mathrm{cac}}\left(F\right){k}_{\mathrm{xce}}\left)\right.E{C}_{e}{k}_{\mathrm{xc}}{C}_{e}\end{array}(18) 
the allograft regenerating antigen indirectly presented subsystem:
\begin{array}{c}\frac{\mathrm{dS}}{\mathrm{dt}}={k}_{\mathrm{s\tau}}+{k}_{s}S{k}_{s}\frac{{S}^{2}}{{S}^{\ast}}{k}_{\mathrm{xsc}}{C}_{l}S\hfill \\ \frac{\mathrm{dL}}{\mathrm{dt}}={k}_{\mathrm{ls}}S{k}_{\mathrm{xlc}}{C}_{l}L{k}_{\mathrm{xl}}L\hfill \\ \frac{d{C}_{l}}{\mathrm{dt}}={k}_{c}+\left(\right.{k}_{\mathrm{cac}}\left(F\right){k}_{\mathrm{xcl}}\left)\right.L{C}_{l}{k}_{\mathrm{xcs}}S{C}_{l}{k}_{\mathrm{xc}}{C}_{l}\end{array}(19) 
the allograft non regenerating antigen directly presented subsystem:
\begin{array}{c}\frac{\mathrm{dU}}{\mathrm{dt}}={k}_{\mathrm{u\tau}}{k}_{\mathrm{xuc}}{C}_{u}U{k}_{\mathrm{xu}}U\hfill \\ \frac{d{C}_{u}}{\mathrm{dt}}={k}_{c}+\left(\right.{k}_{\mathrm{cac}}\left(F\right){k}_{\mathrm{xcu}}\left)\right.U{C}_{u}{k}_{\mathrm{xc}}{C}_{u}\end{array}(20) 
the allograft non regenerating antigen indirectly presented subsystem:
\begin{array}{c}\frac{\mathrm{dY}}{\mathrm{dt}}={k}_{\mathrm{y\tau}}{k}_{\mathrm{xyc}}{C}_{z}Y{k}_{\mathrm{xy}}Y\hfill \\ \frac{\mathrm{dZ}}{\mathrm{dt}}={k}_{\mathrm{xzc}}{C}_{z}Z{k}_{\mathrm{xz}}Z+{k}_{\mathrm{zy}}Y\hfill \\ \frac{d{C}_{z}}{\mathrm{dt}}={k}_{c}+\left(\right.{k}_{\mathrm{cac}}\left(F\right){k}_{\mathrm{xcz}}\left)\right.Z{C}_{z}{k}_{\mathrm{xc}}{C}_{z}{k}_{\mathrm{xcy}}Y{C}_{z}\end{array}(21)
Each subsystem is driven by the common input given by the drug concentration F, which evolves according to a stepwise trajectory:
Lemma
Each state component of the four subsystems endowed with a physiological initial condition (i.e. all positive components), admits nonnegative evolutions, ∀t ≥ 0.
Proof
Consider subsystem (18) and E(0)>0. Due to the continuity of both E(t) and dE/dt, the solution E(t) would become negative if there existed a time instant\stackrel{\u0304}{t}>0 such thatE\left(\stackrel{\u0304}{t}\right)=0 and{\left(\rightseparators="">\frac{\mathrm{dE}}{\mathrm{dt}}}_{}\n \n t\n =\n \n \n t\n \n \u0304\n \n \n, which is a contradiction because:
According to the same reasoning, it also readily appears that C_{ e }(t) never vanishes. The same approach can be repeated for the other three subsystems, since the deltaDirac functions simply model a positive instantaneous increase occurring at the time t_{ τ } of transplantation.•
Lemma
As far as the environmental antigen subsystem (18) is concerned, there exists a unique positive, locally asymptotically stable equilibrium point.
Proof
The equilibrium points of (18) satisfy the following algebraic equations:
from which it follows that the steady state of C_{ e }satisfies the following second order equation:
Since the second and zeroorder coefficients are positive and negative respectively, both solutions are real: one positive, the other negative, regardless to the sign of the firstorder coefficient (see, e.g.,[29]). Thus we have a unique positive solution for C_{ e }. As a matter of fact, by substituting the positive solution into the first equation of (24) we have a unique positive solution also for E:
As for the local stability analysis, we compute the Jacobian matrix:
from which the characteristic polynomial is:
Since all the coefficients are positive, the roots have negative real part[29], which means local asymptotic stability of the equilibrium point. •
Remark
If we interpret the positive solution of (25) as a function of the drug\stackrel{\u0304}{F} administered after transplantation it happens that, by increasing\stackrel{\u0304}{F} the corresponding equilibrium for C_{ e }reduces its value since the first order coefficient increases, keeping unchanged the second and zeroorder terms. And, as a matter of fact, the corresponding equilibrium of E increases its value. These equilibria are drawn in Figure13: increasing values of F go from the bottom right corner, corresponding to the pretransplantation case F=0, to the upper left.
Lemma
The allograft antigen subsystem (19) admits an elementary equilibrium point (S = 0,L = 0,C_{ l } = k_{ c }/k_{ xc }), which is locally asymptotically stable if:
Proof
The existence of the equilibrium point (S = 0, L = 0, C_{ l } = k_{ c }/k_{ xc }) comes out by ready computation. In order to investigate the local stability, compute the Jacobian matrix and evaluate it at the equilibrium point:
It clearly comes that the eigenvalues of J_{ l }are the elements of the diagonal: therefore, the equilibrium point is locally asymptotically stable if condition (29) is satisfied, since in that case all the eigenvalues of the Jacobian matrix are negative real. •
Remark
It is reasonable to assume, from a physiological point of view, that the stability condition (29) is not satisfied. In this case, according to the structure of the Jacobian matrix, it happens that perturbations of the type (S = 0,L = ε_{ L },C_{ l } = k_{ c }/k_{ xc } + ε_{ C l }) allow a trajectory definitely convergent to the equilibrium (S = 0,L = 0,C_{ l } = k_{ c }/k_{ xc }), if{\epsilon}_{L},{\epsilon}_{{C}_{l}} are small enough. On the other hand, for any arbitrarily small ε_{ S }>0, any perturbation of the type (S = ε_{ S },L = 0,C_{ l } = k_{ c }/k_{ xc }) will make the trajectory diverge from equilibrium. Indeed, simulations have been carried out by setting the model parameters in order not to have inequality (29) satisfied (see Table2).
The investigation for other equilibrium points of subsystem (29) requires the computation of the solutions of the following nonlinear algebraic system:
By making substitutions, the previous system becomes:
where:
A qualitative analysis would provide conditions too cumbersome to be easily treated, therefore a numerical bifurcation analysis has been carried out, according to the set of parameters reported in Table2. The bifurcation parameter is the drug amount F.
It turns out that there exists a unique triple of real positive solutions (S,L,C_{ l }) for system (32), whose values are depicted in Figure7 versus the drug concentration F. It has to be stressed that this equilibrium point is locally asymptotically stable whatever the value of F. Figures8,9,10 refer to the bifurcation diagrams for S, L and C_{ l } with respect to the drug concentration F
Lemma
There exists a unique nonnegative, locally asymptotically stable equilibrium point (U = 0, C_{ u } = k_{ c }/k_{ xc }) for the directly presented allograft antigen (20), and there exists a unique nonnegative, locally asymptotically stable equilibrium point (Y = 0, Z = 0, C_{ z } = k_{ c }/k_{ xc }) for the nondirectly presented allograft antigen (21).
Proof
The proof comes from direct computation, following the same lines of the previous Lemmas.•
Author’s contributions
ADG designed the research, wrote the model equations, assessed the numerical results, and participated in writing the paper. AM participated in writing the paper, performance of the research and conducted the MLR experiment. AA participated in the performance of the research. PP participated in writing the paper and conducted the mathematical analysis. FR participated in the performance of the research. SM participated in writing of paper, in the performance of the research and assessed the clinical significance of the results. All authors read and approved the final manuscript.
Abbreviations
 APC:

Antigen Presenting Cell
 CNI:

Calcineurine Inhibitor
 MHC:

Major Histocompatibility Complex
 ODE:

Ordinary Differential Equation
 TCR:

TCell Receptor.
References
Taylor AL, Watson CJ, Bradley JA: Immunosuppressive agents in solid organ transplantation: Mechanisms of action and therapeutic efficacy. Crit Rev Oncol Hematol. 2005, 56 (1): 2346. 10.1016/j.critrevonc.2005.03.012.
Kermack WO, McKendrick AG: A Contribution to the Mathematical Theory of an Epidemic. Proc R Soc Lond A. 1927, 115: 700721. 10.1098/rspa.1927.0118.
Nowak MA, May RM: Virus dynamics: Mathematical principles of immunology and virology. 2001, New York: Oxford University Press
Bajaria SH, Webb G, Cloyd M, Kirschner D: Dynamics of Naive and Memory CD4+ T Lymphocytes in HIV1 Disease Progression. J Acquir Immune Defic Syndr. 2002, 30 (1): 4158.
Marino S, Kirschner DE: The human immune response to Mycobacterium tuberculosis in lung and lymph node. J Theor Biol. 2004, 227 (4): 46386. 10.1016/j.jtbi.2003.11.023.
Matzavinos A, Chaplain MA, Kuznetsov VA: Mathematical modelling of the spatiotemporal response of cytotoxic Tlymphocytes to a solid tumour. Math Med Biol. 2004, 21 (1): 134. 10.1093/imammb/21.1.1.
Flower DR, Timmis J: In silico immunology. 2007, New York: Springer
Ciupe SM, Devlin BH, Louise Markert M, Kepler TB: The Dynamics of TCell Receptor Repertoire Diversity Following Thymus Transplantation for DiGeorge Anomaly. PLoS Computational Biology. 2009, 5: 113.
Forsdyke DR: Early evolution of MHC polymorphism. J Theor Biol. 1991, 150: 4516. 10.1016/S00225193(05)804394.
Cote I, Rogers NJ, Lechler RI: Allorecognition. Transfus Clin Biol. 2001, 8 (3): 31823. 10.1016/S12467820(01)001215.
Heeger PS: Tcell allorecognition and transplant rejection: a summary and update. Am J Transplant. 2003, 3 (5): 52533. 10.1034/j.16006143.2003.00123.x. Review
Gould DS, Auchincloss Jr H: Direct and Indirect Allorecognition: the role of MHC antigens in graft rejection. Immunol Today. 1999, 20 (2): 7782. 10.1016/S01675699(98)013942.
Benichou G, Takizawa PA, Olson CA, McMillan M, Sercarz EE: Donor Major Histocompatibility Complex (MHC) Peptides Are Presented by Recipient MHC Molecules during Graft Rejection. J Exp Med. 1992, 175 (1): 3058. 10.1084/jem.175.1.305.
Yewdell JW, Bennink JR: Immunodominance in Major Histocompatibility Complex classI restricted T lymphocyte responses. Annu Rev Immunol. 1999, 17: 5188. 10.1146/annurev.immunol.17.1.51. Review
Benichou G, Valujskikh A, Heeger PS: Contributions of Direct and Indirect T Cell Alloreactivity During Allograft Rejection in Mice. J Immunol. 1999, 162 (1): 3528.
Faderl S, Estrov Z: Hematopoietic Growth Factors and Cytokines. J Clin Ligand Assay. 2000, 23: 165176.
Abbas AK, Litchman AH: Immunologia cellulare e molecolare. 5th edition. Italian edition. Edited by: Sovazzani S, Masson Ed. 2006
Otton R, Soriano FG, Verlengia R, Curi R: Diabetes induces apoptosis in lymphocytes. J Endocrinol. 2004, 182: 145156. 10.1677/joe.0.1820145.
Mager DE, Jusko WJ: Development of Translational PharmacokineticPharmacodynamic Models. Clin Pharmacol Ther. 83 (6): 90912. Epub 2008 Mar 26. Review
Safinia N, Afzali B, Atalar K, Lombardi G, Lechler RI: Tcell alloimmunity and chronic allograft dysfunction. Kidney Int. 2010, 78 (Suppl 119): S2—S12
Torres AJL, Angelo ALD, Netto EM, Sampaio GP, Souza DFC, Inocêncio LA, Lemos JAR, Brites C: Reference Range for T Lymphocytes Populations in Blood Donors from Two Different Regions in Brazil. Brazilian Journal Infectious Diseases. 2009, 13 (3): 221225. 10.1590/S141386702009000300013.
Post DJ, Douglas DD, Mulligan DC: Immunosuppression in Liver Transplantation. Liver Transplantation. 2005, 11 (11): 13071314. 10.1002/lt.20614.
Humar A, Michael M: American Society of Transplantation Recommendations for Screening, Monitoring and Reporting of Infectious Complications in Immunosuppression Trials in Recipients of Organ Transplantation. Am J Transplantation. 2006, 6: 262274. 10.1111/j.16006143.2005.01207.x.
Morelon E, Touraine JL: Infectious complications due to immunosuppression in organ transplant patients. La Revue du Praticien. 2007, 57 (15): 16771686.
Ruggeri A, Martinelli M: A program for the optimization of cyclosporine therapy using population kinetics modeling. Comput Methods Progras Biomed. 2000, 61: 6169. 10.1016/S01692607(99)000280.
KeHua W, YiMin C, JinFeng G, Ying Z, SuoDi Z, FuDe C, Wei L: Population pharmacokinatics of cyclosporine in clinical renal transplant patients. Drug Metabolism and Disposition. 2005, 33: 12681275. 10.1124/dmd.105.004358.
Marchetti P: Newonset diabetes after transplantation. J Heart Lung Transplant. 2004, 23 (5 Suppl): S194201.
Gonçalves SC, DíazSerrano KV, de Queiroz AM, Palioto DB, Faria G: Gingival overgrowth in a renal transplant recipient using cyclosporine A. J Dent Child (Chic). 2008, 75 (3): 3137.
Franklin GF, Powell JD, EmamiNaeini A: Feedback control of dynamic systems. 1994, MA: Addison Wesley
Acknowledgements
We thank Chiara Niccolò who helped with the MLR experiment and Annabella Procoli who provided the cytofluorimeter.
Author information
Authors and Affiliations
Corresponding author
Additional information
Competing interests
The authors declare that they have no competing interests.
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
Rights and permissions
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
About this article
Cite this article
Gaetano, A., Matone, A., Agnes, A. et al. Modeling rejection immunity. Theor Biol Med Model 9, 18 (2012). https://doi.org/10.1186/17424682918
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/17424682918
Keywords
 Major Histocompatibility Complex
 Clonal Expansion
 Chronic Rejection
 Major Histocompatibility Complex Molecule
 Mixed Lymphocyte Reaction