## ORIGINAL RESEARCH

Transpl Int, 11 April 2022
https://doi.org/10.3389/ti.2022.10297

# Modeling the Potential of Treg-Based Therapies for Transplant Rejection: Effect of Dose, Timing, and Accumulation Site

Maya M. Lapp1, Guang Lin2, Alexander Komin3, Leah Andrews4, Mei Knudson5, Lauren Mossman4, Giorgio Raimondi3* and Julia C. Arciero6*
• 1Department of Mathematics, The College of Wooster, Wooster, OH, United States
• 2Department of Mathematics, Purdue University, West Lafayette, IN, United States
• 3Department of Plastic and Reconstructive Surgery, Johns Hopkins School of Medicine, Baltimore, MD, United States
• 4Department of Mathematics, St. Olaf College, Northfield, MN, United States
• 5Department of Mathematics, Carleton College, Northfield, MN, United States
• 6Department of Mathematical Sciences, Indiana University-Purdue University of Indianapolis, Indianapolis, IN, United States

Introduction: The adoptive transfer of regulatory T cells (Tregs) has emerged as a method to promote graft tolerance. Clinical trials have demonstrated the safety of adoptive transfer and are now assessing their therapeutic efficacy. Strategies that generate large numbers of antigen specific Tregs are even more efficacious. However, the combinations of factors that influence the outcome of adoptive transfer are too numerous to be tested experimentally. Here, mathematical modeling is used to predict the most impactful treatment scenarios.

Methods: We adapted our mathematical model of murine heart transplant rejection to simulate Treg adoptive transfer and to correlate therapeutic efficacy with Treg dose and timing, frequency of administration, and distribution of injected cells.

Results: The model predicts that Tregs directly accumulating to the graft are more protective than Tregs localizing to draining lymph nodes. Inhibiting antigen-presenting cell maturation and effector functions at the graft site was more effective at modulating rejection than inhibition of T cell activation in lymphoid tissues. These complex dynamics define non-intuitive relationships between graft survival and timing and frequency of adoptive transfer.

Conclusion: This work provides the framework for better understanding the impact of Treg adoptive transfer and will guide experimental design to improve interventions.

GRAPHICAL ABSTRACT

GRAPHICAL ABSTRACT

## Introduction

Following transplantation, lifelong immunosuppression is required to prevent allograft rejection (1). Unfortunately, because of drug-associated complications (13) and the unfeasibility of complete immunosuppression, 5 years survival rates for patients undergoing solid organ transplantation range anywhere from 40–70%, depending on the organ transplant type (46). Thus, there is an urgent need to develop alternative treatment strategies to promote graft tolerance and improve the quality of life for transplant recipients.

One such promising strategy is the adoptive transfer of regulatory T cells (Tregs). Tregs are lymphocytes that suppress the activity of other immune cells and are critical for maintaining peripheral tolerance and preventing autoimmune pathologies (7). The idea is that adoptive transfer of large numbers of Tregs can suppress transplant rejection and promote establishment of tolerance to the transplanted organ. Studies in pre-clinical animal models have shown feasibility and efficacy of polyclonal Treg infusion (813). Results from the completed clinical trials provide evidence that delivery of expanded polyclonal Treg is safe (1416) and possibly effective (17). Excitingly, preclinical studies suggest that the therapeutic effect can be improved by using alloantigen-specific Tregs (1820). To achieve this, Tregs can be modified to express Chimeric Antigen Receptors (CAR) or transgenic T cell receptors, which endow large-scale production of Tregs with the desired antigen specificity (10, 21, 22). In animal transplant models, CAR-Tregs have been shown to reach the grafts (19, 23, 24) and to control skin graft rejection to a greater extent than polyclonal Tregs (19, 24). CAR-Treg application and safety is currently being investigated in the first clinical trial in Europe (STEADFAST study, Sangamo Therapeutics, Inc.).

Despite these promising outcomes, much remains unknown about the consequences and efficacy of Treg adoptive transfer and what conditions maximize their therapeutic effect. It is difficult to compare experimental outcomes between studies due to differences in model organisms, transplanted organs, dose magnitude, dose timing, and Treg population quality (affected by in vitro expansion). Moreover, disparities in the results of some Phase I clinical trials highlight this lack of understanding of their optimal use (13, 25, 26). Further clinical trials will continue to improve our comprehension of adoptive transfer efficacy. However, it takes years to evaluate the long-term effects of this treatment, and development of methods to facilitate a more rapid understanding of the impact and optimization of treatment regimen would be invaluable.

Mathematical and computational models have been widely used in conjunction with experimental methods in cancer and virology to understand immune system dynamics and aid in the design of effective immunotherapies (2729), but their application in transplantation is lacking. Some theoretical models for solid organ transplant rejection have been proposed (3035), several of which focused on the impact of immunosuppression only (31, 3335) and used simplified representations of immune components. We established one of the first theoretical models to describe the immune system and transplant dynamics that give rise to transplant rejection, and ours is currently the only transplantation model using differential equations to track Tregs independently from other T cell populations (30). Both our work and that of De Gaetano et al. incorporated experimental methods in developing transplant rejection models (30, 31). However, no mathematical model to date provides a robust mechanism for analysis of Treg adoptive transfer on the immune response to transplantation.

In the present study, we adapted our model (30) to include new equations and terms to simulate the impact of alloantigen-specific Treg adoptive transfer on graft survival (Supplementary Figure S1). We updated the model equations for antigen-presenting cells (APCs) to match their behavior observed in vivo (36). We also performed parameter estimation on the updated model using a non-dominated sorting-based multi-objective evolutionary algorithm (37), called NSGA-II (Non-dominated Sorting Genetic Algorithm II). Although adoptive transfer is not sufficient to prevent graft rejection independently of immunosuppression (10), in this study we focus exclusively on adoptive transfer treatments to elucidate their direct effect on transplant rejection without the complicating and possibly confounding effects of concomitant immunosuppression. Thus, using the updated model, this study aimed at: 1) identifying optimal conditions for Treg delivery, specifically the activation status and tissue distribution, magnitude of dose, timing of delivery, and frequency of dosing, 2) analyzing immune dynamics to explain the effects of adoptive transfer on graft survival, and 3) suggesting future avenues for experimental studies into adoptive transfer treatment.

Our model simulations and analysis identified that timely inhibition of dendritic cell (DC) maturation in the graft and prolonged modulation of cytotoxic CD8 T cells and inflammatory macrophage activity in the graft by Treg were significantly more impactful than inhibiting the activation of alloreactive T cells in draining lymphoid tissues. Use of the model allowed us to identify a non-intuitive correlation between Treg dosing and administration frequency that delineates more effective interventions. Overall, our model enables the rapid simulation of a vast number of conditions that would be prohibitive to cover experimentally. It also provides the framework for future modeling efforts that will assess combinatorial treatments involving both Treg adoptive transfer and immunoregulatory agents.

## Materials and Methods

This study does not involve any new animal or human studies, and thus it was exempt from IRB approval.

### Model Description

We expanded our previously developed model describing the dynamics of murine heart transplant rejection (30) to examine the impact of adoptive transfer of Tregs on graft survival. The original model consists of 13 nonlinear first order ordinary differential equations (ODEs) tracking the following populations in a representative lymph node compartment (LN) and graft compartment (G): CD4 T cells ($THLN$ and $THG)$, CD8 T cells ($TELN$ and $TEG)$, regulatory T cells ($TRLN$ and $TRG)$, immature DCs $(Aimm)$, mature DCs ($AmatLN$ and $AmatG)$, inflammatory macrophages ($Ainf$), pro- and anti-inflammatory cytokines ($Cp$ and $Ca)$, and graft cells ($G$). An additional ODE is introduced to track naïve Tregs ($TRNLN$). All model equations and parameters appear in the Supplemental Digital Content, and all model variables and their initial values are provided in Supplementary Table S1. A complete list of model assumptions is provided in (30). While this model is parametrized for mouse heart transplantation, the predicted interactions between the immune response and transplanted organ would be similar for other transplants.

In this study, the model equations tracking the immature and mature DCs are adapted to include more realistic representations of experimental observations (36). Specifically, the decay rate of immature DCs is assumed to depend on the remaining graft mass (Supplementary Eq. S10, second term). In addition, DC maturation is assumed to occur in the presence of pro-inflammatory cytokines or CD4 T cells (Supplementary Eq. S10, third term; Supplementary Eq. S11, first term). For verification of the modified model and updated predictions for host immune dynamics without adoptive transfer treatment, see the Supplemental Digital Content (Supplementary Figures S2, S3).

A dosing function, $D(t)$, for the adoptive transfer of Tregs is defined in Eq. 1 (also in Supplementary Eq. S12). Parameters $fG$, $fLN$, and $fN$ correspond to the fraction of the Treg dose that enters the graft as activated Tregs (Eq. 2), the lymph node as activated Tregs (Eq. 3), and the lymph node as naïve Tregs (Eq. 4), respectively. For simulations involving adoptive transfer, $fG+fLN+fN=1$. These three parameters are set to zero when adoptive transfer is not simulated. The dosing function is composed of one or multiple exponentially decaying functions of time, where D0 indicates the dosing rate, n indicates the total number of doses delivered, $ti$ indicates the post-operative day (POD) of delivery of the ith Treg injection, and $β$ indicates the decay rate. A decay rate of $β=2$ day−1 is used to simulate the relatively fast absorption of Tregs during adoptive transfer. We assume around 50% of injected cells are lost or distribute in non-draining lymphoid tissues (which are excluded from the model), and we define the dose magnitude $C$ to be the total number of injected cells that localize to the modeled lymph node and graft compartments. Thus, for the majority of this study, we simulate adoptive transfer with Tregs. The dose magnitude is the area under the curve of our dosing function, which is given by cells. A single dose with cells/day or $n$ doses with cells/day both administer a total of Tregs. Figure 1 shows the dosing rate over time for a single dose administered on POD0 with cells/day (Figure 1A) and 5 doses administered on POD0, 1, 2, 3, and 4 with $D0=1065$ cells/day (Figure 1B). The dose magnitude $C$ is varied during model analysis.

$dTRGdt=keRTRLN−μRTRG+rRGTRG(TEG+THG)α6+TRG+fGD(t)(2)$
$dTRLNdt=aRTRNLNAmatLNγ2+AmatLN−μRTRLN+rRTRLN(TELN+THLN)α2+TRLN−eRTRLN+fLND(t)(3)$

FIGURE 1

FIGURE 1. Visual representations of Dosing function, $D(t)$, delivering a total of Tregs. (A) A single dose with dosing rate $D0=106$ cells/day. (B) Five doses delivered on POD0-4 with dosing rate $D0=1065$ cells/day.

### Model Simulations

The activation state, accumulation site, magnitude, timing, and frequency of adoptive transfer are varied in the model to simulate the impact of each of these factors on the immune response to the graft. Parameters $fG$, $fLN$, and $fN$ are varied to determine the impact of activation state (i.e., naïve or activated) and accumulation site of Tregs. The effects of dose magnitude are assessed by varying $C$. The timing and frequency of Treg doses are evaluated by varying the start day of the dose ($t0$) and the number of doses, n. If multiple doses are administered, it is assumed that the doses are given at 1-day intervals. Following our previous work, graft rejection is defined as a 75% reduction in the original number of graft cells (30). The time at which the model predicts graft rejection is used in this study to identify optimal dosing strategies for the adoptive transfer of Tregs.

### Parameter Estimation

To calibrate the model parameters in the updated model (see parameters shaded in Supplementary Table S2), we employed a non-dominated sorting-based multi-objective evolutionary algorithm (37), called NSGA-II (Non-dominated Sorting Genetic Algorithm II), which can greatly improve the efficiency in constrained multi-objective optimization tasks. NSGA-II is a popular non-domination based genetic algorithm for multi-objective optimization and parameter estimation. NSGA-II improves elitism and there is no need to choose sharing parameters a priori.

## Results

### Impact of Treg Accumulation Site

To analyze the effect of Treg adoptive transfer using our model, we first evaluated the impact of the activation status of the injected cells on graft survival. It is well recognized that pre-activated Treg (even after resting) are an order of magnitude more suppressive than naïve Tregs (38, 39). For this analysis, we simulated the equivalent of delivering Tregs on POD0. As expected, accumulation of naïve Tregs into the lymph node was not as effective as the adoptive transfer of pre-activated Tregs (Supplemental Digital Content, Supplementary Figure S4). This confirmed the benefit of using ex-vivo expanded Tregs, which is a necessary step in most clinical preparations of antigen-specific cells that renders homogeneous populations of pre-activated Treg. For the remainder of the study, we simulated the use of pre-activated Tregs.

Next, we focused on the question: does the tissue distribution of Tregs between the lymphoid compartment and the graft post-adoptive transfer impact transplant survival? Figure 2A depicts the change in the graft mass over time for the scenarios where the transferred Tregs exclusively accumulate in the lymphoid versus graft compartments, as well as for the case of no adoptive transfer (black curve). There is a substantial benefit with the accumulation of Tregs in the graft in comparison to lymphoid tissues, with estimated rejection on POD74 vs. POD25, respectively. The model can simulate a wide range of Tregs dose magnitudes; for almost all ranges, accumulation of Tregs directly to the graft (Figure 2B, blue curve) is more effective at extending survival than localization to the lymph node compartment (Figure 2B, red curve). For very small doses ($C<2.9×104$ cells), delivering activated Tregs to the lymph node is predicted to be more effective than delivering to the graft (Figure 2C); though, it improves graft survival by no more than 2.4 days.

FIGURE 2

FIGURE 2. Impact of distribution of adoptively transferred Tregs. (A) Number of graft cells shown over time for three different cases: activated Tregs administered to the graft (blue), activated Tregs administered to the lymph node (red), and no Tregs administered (black). Tregs are administered on POD0 with dose magnitude  cells. The horizontal dashed line indicates a 75% reduction in initial graft size. The intersections of the curves with this dashed line give the model predicted values of graft rejection time. (B) Model predicted rejection time as Treg dose magnitude ($C$) is varied. The different curves depict predictions for the delivery of activated Tregs to the graft (blue) and activated Tregs to the lymph node (red). The black line marks POD11, which is when rejection occurs without adoptive transfer. (C) Circled region in panel B is magnified to examine the effect of small dose magnitudes on graft rejection.

### Impact of Single Dose Timing

Having shown the benefit of graft localization, we studied the impact of the timing of a single injection of varying amounts of Tregs. Figure 3 shows a non-monotonic relationship between rejection time and day of Treg administration, observed for some Treg dose magnitudes (Figure 3A). Specifically, when using cells (blue curve), we observed that delaying Treg administration to POD1.5 yields optimal graft survival time. The non-monotonic behavior between dose timing and graft survival is highlighted for physiological dose magnitudes ($C<106$ cells) in Figure 3B. Once the dose magnitude is increased above cells, a monotonic relationship is re-established in which graft survival time decreases with the delay of Treg administration.

FIGURE 3

FIGURE 3. Impact of the timing of Treg administration on graft rejection. (A) Predicted graft rejection time as the day of dose administration (ti) is varied given four different dose magnitudes: (black), (blue), (red), and (green) cells. (B) Predicted graft rejection time as dose magnitude ($C$) is varied for Tregs administered on POD0 (blue), POD1.5 (red), and POD3 (green).

Although a theoretical model can assume that all Tregs administered to an individual accumulate exclusively in the graft or lymphoid tissues, physiological in vivo constraints and Treg properties dictate a variable partitioning between the two compartments. Thus, in Figure 4, we show the simulated impact of varying the fraction of Tregs that enter the graft (fg) between 0 and 1; the remaining fraction (1–fg) is assumed to enter the draining lymphoid tissue compartment. As indicated by the red curve, the model predicts that if more than 10% of the injected Tregs locate to the graft, delaying injection until POD1.5 is the most beneficial strategy to prolong the time to graft rejection.

FIGURE 4

FIGURE 4. Impact on rejection time of the fraction of activated Tregs that enters the graft. A dose of cells is administered on POD0 (blue), POD1.5 (red), and POD3 (green) and the impact on graft survival is projected in relation to the distribution of the injected Treg between graft and lymphoid compartments.

### Impact of Multiple Doses

The model allows us to compare the protective effect of a single injection of Tregs versus splitting that total number of cells among multiple daily injections (an approach that can also represent sustaining the presence of a smaller amount of Treg over time). As depicted in Figure 5, distributing a fixed dose prolongs graft survival. The relationship between the number of consecutive doses and graft survival was not monotonic, with graft survival peaking around 10–15 equivalent doses, depending on the first day of injection (Figure 5A). Significantly, when enough doses are administered, delaying the start of adoptive transfer treatment is no longer beneficial. The model predicts that administering cells divided evenly among 14 daily doses beginning on POD0 is the most effective treatment, extending graft survival until POD115. As shown in Figure 5B (graft mass over time), it is noteworthy that in addition to extending transplant survival, multiple doses of Treg provide better protection of the graft (i.e., slower rate of mass decline). Figure 5C shows that the impact of multiple doses compared with a single dose is apparent even if only a small fraction of Tregs (>5%) enters the graft.

FIGURE 5

FIGURE 5. Impact on graft rejection of splitting the same total number of Treg into multiple doses. (A) Plot of the simulated rejection time as a function of the indicated number of doses administered starting on POD0 (blue), POD1 (red), POD2 (black), and POD3 (green). The total dose magnitude for each simulation is cells. (B) Number of graft cells as a function of time for a single dose (blue curve) and for the optimal case of multiple doses (14 daily doses of Tregs beginning on POD0, red curve). (C) Comparison of rejection time for single (blue curve) and multiple (red curve) dose administration if the fraction of Tregs that reach the graft is varied.

### Impact of Treg Adoptive Transfer on Host Immune Dynamics

Differently from bioinformatic approaches (e.g., machine learning), mathematical modeling enables the analysis of the dynamics of host immune cells that determine the simulated outcome. This valuable property can be used to gain insight into why certain adoptive transfer treatments are more effective than others. We investigated in detail the difference in effects between using a single injection of graft-infiltrating Tregs versus splitting the total number over multiple injections. Figure 6 depicts the predicted dynamics of host immune cells in the graft with single dose adoptive transfer (blue) and 14 dose adoptive transfer (red). For a similar investigation into the impact of site of accumulation and timing of single injections, please see the Supplemental Digital Content (Supplementary Figures S5, S6).

FIGURE 6

FIGURE 6. Model predicted alterations of immune dynamics induced by Tregs delivered as a single dose (blue) or distributed across 14 doses (red). (A) Number of Tregs in the graft, $TRG$. (B) Number of mature DCs in the graft, $AmatG$. (C) Number of cytotoxic CD8 T cells in the graft, $TEG$. (D) Number of inflammatory macrophages in the graft, $AinfG$. (E) Rate of graft destruction caused by cytotoxic CD8 T cells. (F) Rate of destruction caused by inflammatory macrophages. Dashed boxes [in (C–F)] highlight timeframes when indicated destructive cell quantity is higher for the 14 dose treatment but the resulting rate of destruction is lower than for the single dose treatment.

Delivering a large single dose of Tregs minimizes the maturation of graft infiltrating DCs (Figure 6B) and therefore limits the activation and graft accumulation of cytotoxic CD8 T cells (Figure 6C) and delays the activation of macrophages (Figure 6D). Therefore, delivering a substantial portion of Tregs soon after transplantation is effective at limiting the early inflammatory response. Nevertheless, the regulation of the destructive capacity of those cells is not sustained with a single administration.

When the injected cells are split among 14 doses, elevated Treg levels are maintained during the surge of cytotoxic T cells and inflammatory macrophages (Figures 6A,C,D, red). In contrast, for single dose adoptive transfer, Tregs in the graft decay before the surge of these graft destructive subsets (Figures 6A,C,D, blue). The lack of Treg-mediated control in the last scenario allows each cytotoxic T cell and inflammatory macrophage to cause greater graft damage despite a lower accumulation level. This phenomenon is highlighted in the boxed regions in Figures 6C–F. We note that Figures 6E,F show the value of each term on the right-hand side of the differential equation for graft cells that contributes to graft destruction (Supplementary Eq. S9). Therefore, higher values of the curve in panel E or F indicate a greater contribution to graft destruction. For example, from POD7.2 to POD11.9, although more cytotoxic T cells are predicted to be present with the 14-dose regimen (Figure 6C; the 14-dose curve is above the 1-dose curve in the dotted boxed region), their effector functions are regulated and less graft destruction occurs (Figure 6E; the 1-dose curve is above the 14-dose curve). Similarly, from POD5.3 to POD7.7, although more inflammatory macrophages are predicted with the 14-dose regimen (Figure 6D), less graft destruction is mediated by this subset (Figure 6F).

Overall, analyzing the population dynamics for various Treg dose frequencies indicates that, to optimize adoptive transfer, cell delivery must 1) maintain elevated Treg levels during the surge of cytotoxic T cells and inflammatory macrophages in the graft and 2) deliver enough Tregs early after transplantation to inhibit DC maturation and therefore limit T cell activation. These same principles hold true when examining dose timing and can explain why administering a single dose on POD1.5 is more effective than delivering a single dose on POD0 (see Supplemental Digital Content, Supplementary Figure S6). When examining site of Treg accumulation, minimizing DC maturation is again critical to understand the benefit of delivering Tregs to the graft rather than the lymph node (see Supplemental Digital Content, Supplementary Figure S5).

## Discussion

In this study, we expanded our previous mathematical model of murine heart transplant rejection to incorporate the impact of Treg adoptive transfer. Importantly, and complementary to bioinformatic approaches (4043), we took advantage of the power of mathematical modeling to enable analysis of the immune dynamics underlying the results of simulations and thereby offer rationales for the therapeutic optimizations proposed. With the advent of feasible and reliable approaches of genetic engineering to generate high numbers of antigen-specific Treg for adoptive transfer (10), our data suggest important optimizations that would maximize the therapeutic efficacy of these cells.

### Impact of Treg Accumulation Site

Our model predicts that Tregs accumulating to the transplant are more therapeutically effective than Tregs distributing to the lymph node (Figure 2). Importantly, analysis of immune dynamics provides the rationale behind such a different outcome: inhibition of DC maturation in the graft is more effective than direct inhibition of T cell activation in the draining lymphoid tissue at minimizing the number of CD8 and CD4 T cells that reach the graft (Supplemental Digital Content, Supplementary Figure S5). Interestingly, experimental reports in rodent models have suggested that an immediate accumulation of Treg in the graft promotes longer survival (44, 45). Overall, our analysis suggests that identifying methods to increase the fraction of Tregs that translocate early on to the graft will maximize the impact of adoptive transfer. It is noteworthy that the conditions of ex-vivo Treg expansion impart specific migratory capacity to the cells and these differences impact the therapeutic result obtained (4549). Moreover, the breadth of genetic engineering that has become possible for clinical products (e.g., CRISPR/Cas9 based modifications) could provide useful tools to impart the ideal behavior in adoptively transferred Tregs.

### Impact of Single Dose Timing

When administering a single dose of Tregs, our simulations suggest that delaying cell delivery is more effective than administering cells on the day of transplantation (Figure 3B). The justification for this unexpected result is similar to that given for splitting Treg into multiple doses. Slightly delaying administration allows Tregs to directly inhibit the destructive functions of cytotoxic T cells and inflammatory macrophages while ensuring a timely reduction of DC activation and, therefore, T cell activation (Supplemental Digital Content, Supplementary Figure S5). Obviously, this scenario is specific to the unique conditions of administering Tregs without any additional manipulation of the recipient (like immunosuppression, see below), but it highlights the impact of examining cellular interactions at the system level.

### Impact of Multiple Doses

The model predicts that maintaining a prolonged influx of Tregs is the most effective method of lengthening graft survival. Delivering multiple, smaller doses of Tregs preserves the graft longer than delivering the same total number of Tregs in a single bolus. This regimen allows for elevated levels of Tregs to remain in the graft during the surge of graft-destructive cells (and control them), while also compromising on an “early enough” limitation of DC maturation. Although the clinical implementation of multiple daily dosing of Treg is unrealistic, this result highlights the important point of sustaining the survival and function of transferred Tregs. This is a contentious issue with the reported negative effect that many immunosuppressive drugs have on the homeostasis and function of Treg (50, 51). There is growing interest in devising approaches to sustain the persistence and function of Tregs. The use of IL-2/anti-IL2 complexes, IL-2 muteins, as well as the genetic engineering of Tregs to respond to “orthogonal” IL-2 are all examples of active investigations to promote Tregs persistence (10, 20, 52). In parallel, the utilization of biomaterials to promote the sustained accumulation of Tregs in proximity to the transplant represent a very promising strategy (53). Our model results provide the rationale to strongly support these ongoing efforts.

### Limitations, Parallelisms, and Future Work

While the processes of sensitivity analysis and parameter estimation performed in this study (based on experimental data) have improved the accuracy of several model parameters, some parameters remain uncertain. For example, the persistence and proliferation of Tregs in the graft are unquantified variables that have a profound impact on the dynamics of graft infiltrating immune cells. Similarly, the relationship between number of Tregs injected and the number of cells that reach the graft is not quantified; the factors that influence such a relationship are poorly defined, and thus, further experimentation is needed (54).

While ours is the only ODE transplantation model to date that tracks Tregs independently of other T cells, several immunological mathematical models have recently been developed to elucidate the role of Tregs on self-tolerance and to identify key Treg interactions with other immune populations (5558). The assumptions from these models may be used to improve our existing model to better replicate Treg behavior. In particular, given the importance of IL-2 to Treg survival, proliferation, and function, explicitly tracking IL-2 concentrations as in (57) would strengthen our current model and allow further investigation into IL-2 therapies as a method of extending Treg survival.

Although there is very limited quantitative information, other pre-clinical transplant models show important concordance with some of our model simulations. In a mouse model of pancreatic islet transplantation (45), Zhang et al. compared the i.v. infusion of Treg with the co-transplantation at the site of grafting. The significantly longer survival obtained in the latter case agrees with our model-suggested principle of a higher therapeutic impact when Treg can rapidly and directly modulate graft immune populations. Unfortunately, their report did not present different doses or timing of Treg administration. A similar scenario emerges from reports using Treg in mouse skin transplant models. The use of CAR-Treg (59), cell lines of Treg (60), or polyclonal Treg (61) injected on the day of the transplant promotes only a modest increase in transplant survival, a result that correlates with simulation of a fixed number of Treg that probably mostly home to the lymphoid compartment. In all cases, the combination with so called “adjunct therapies” (ranging from thymectomy and T cell depletion to irradiation and bone marrow co-transplantation) is demonstrated as necessary to achieve lasting therapeutic effects. Currently, we can only draw qualitative comparisons to these experimental models since theoretical model parameters would need to be adapted to each different transplant model.

Overall, our modified model is beneficial in identifying methods to maximize the benefits of Treg adoptive transfer and in generating hypotheses on the key immune dynamics that govern its outcome and that can be tested experimentally. However, as demonstrated in this study and by experimental evidence to date, Treg adoptive transfer alone is insufficient to prevent transplant rejection. These results highlight the need to understand what additional perturbations to the system would better support or enhance the protective function of Tregs. Our mathematical model provides the framework into which treatments like immunosuppression (existing or hypothetical) can be included and used to dissect their very complex effects. Combined with the promising technological advances in both the investigation and manipulation of cells, there is tangible optimism toward the ultimate goal of optimizing therapeutic strategies for transplantation.

## Data Availability Statement

The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding authors.

## Author Contributions

JA and GR conceived the idea of the study; all authors contributed to the rationalization of immunological concepts into equations; JA, ML, LA, LM, and MK created and optimized the model code; GL and LA performed model optimization; GR and JA defined immunological assumptions; JA, ML, GR, and AK wrote the manuscript.

## Funding

This work was supported by United States Army Medical Research Acquisition Activity (USMRAA) grant W81XWH-19-1-0352, and internal support from the Department of Plastic and Reconstructive Surgery (both to GR). This work was also supported by NSF DMS-1559745, NSF DMS-1654019, and NIH 1R01EY030851 (to JA).

## Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

## Acknowledgments

The authors thank Olga Dorabiala, Ross Walgren, Julia Shapiro, Hannah Schmidt, and Anirudh Arun for invaluable support with literature search and conceptualization of immune interactions.

## Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontierspartnerships.org/articles/10.3389/ti.2022.10297/full#supplementary-material

## Abbreviations

APC, antigen-presenting cell; AT, adoptive transfer; CAR, chimeric antigen receptors; DC, dendritic cell; IL-2, interleukin 2; LN, lymph node; ODE, ordinary differential equation; POD, post-operative day; Tregs, regulatory T cells.

## References

1. Allison, TL. Immunosuppressive Therapy in Transplantation. Nurs Clin North America (2016) 51(1):107–20. doi:10.1016/j.cnur.2015.10.008

2. Min, DI, and Monaco, AP. Complications Associated with Immunosuppressive Therapy and Their Management. Pharmacotherapy (1991) 11(5):119S–125S. doi:10.1002/j.1875-9114.1991.tb02641.x

3. de Fijter, JW, Holdaas, H, Øyen, O, Sanders, J-S, Sundar, S, Bemelman, FJ, et al. Early Conversion from Calcineurin Inhibitor- to Everolimus-Based Therapy Following Kidney Transplantation: Results of the Randomized ELEVATE Trial. Am J Transpl (2017) 17(7):1853–67. doi:10.1111/ajt.14186

4. Shyu, S, Dew, MA, Pilewski, JM, DeVito Dabbs, AJ, Zaldonis, DB, Studer, SM, et al. Five-year Outcomes with Alemtuzumab Induction after Lung Transplantation. J Heart Lung Transplant (2011) 30(7):743–54. doi:10.1016/j.healun.2011.01.714

5. Grant, D, Abu-Elmagd, K, Mazariegos, G, Vianna, R, Langnas, A, Mangus, R, et al. Intestinal Transplant Registry Report: Global Activity and Trends. Am J Transpl (2015) 15(1):210–9. doi:10.1111/ajt.12979

6. Hellemans, R, Hazzan, M, Durand, D, Mourad, G, Lang, P, Kessler, M, et al. Daclizumab versus Rabbit Antithymocyte Globulin in High-Risk Renal Transplants: Five-Year Follow-Up of a Randomized Study. Am J Transpl (2015) 15(7):1923–32. doi:10.1111/ajt.13191

7. Vignali, DAA, Collison, LW, and Workman, CJ. How Regulatory T Cells Work. Nat Rev Immunol (2008) 8(7):523–32. doi:10.1038/nri2343

8. Safinia, N, Leech, J, Hernandez-Fuentes, M, Lechler, R, and Lombardi, G. Promoting Transplantation Tolerance; Adoptive Regulatory T Cell Therapy. Clin Exp Immunol (2013) 172(2):158–68. doi:10.1111/cei.12052

9. Tang, Q, and Bluestone, JA. Regulatory T-Cell Therapy in Transplantation: Moving to the Clinic. Cold Spring Harb Perspect Med (2013) 3(11):a015552. doi:10.1101/cshperspect.a015552

10. Ferreira, LMR, Muller, YD, Bluestone, JA, and Tang, Q. Next-generation Regulatory T Cell Therapy. Nat Rev Drug Discov (2019) 18(10):749–69. doi:10.1038/s41573-019-0041-4

11. MacDonald, KN, Piret, JM, and Levings, MK. Methods to Manufacture Regulatory T Cells for Cell Therapy. Clin Exp Immunol (2019) 197(1):52–63. doi:10.1111/cei.13297

12. Issa, F, Strober, S, Leventhal, JR, Kawai, T, Kaufman, DB, Levitsky, J, et al. The Fourth International Workshop on Clinical Transplant Tolerance. Am J Transpl (2021) 21(1):21–31. doi:10.1111/ajt.16139

13. Sawitzki, B, Harden, PN, Reinke, P, Moreau, A, Hutchinson, JA, Game, DS, et al. Regulatory Cell Therapy in Kidney Transplantation (The ONE Study): a Harmonised Design and Analysis of Seven Non-randomised, Single-Arm, Phase 1/2A Trials. The Lancet (2020) 395(10237):1627–39. doi:10.1016/s0140-6736(20)30167-7

14. Trzonkowski, P, Bacchetta, R, Battaglia, M, Berglund, D, Bohnenkamp, HR, ten Brinke, A, et al. Hurdles in Therapy with Regulatory T Cells. Sci Transl Med (2015) 7(304):304ps18. doi:10.1126/scitranslmed.aaa7721

15. Brunstein, CG, Miller, JS, Cao, Q, McKenna, DH, Hippen, KL, Curtsinger, J, et al. Infusion of Ex Vivo Expanded T Regulatory Cells in Adults Transplanted with Umbilical Cord Blood: Safety Profile and Detection Kinetics. Blood (2011) 117(3):1061–70. doi:10.1182/blood-2010-07-293795

16. Bluestone, JA, Buckner, JH, Fitch, M, Gitelman, SE, Gupta, S, Hellerstein, MK, et al. Type 1 Diabetes Immunotherapy Using Polyclonal Regulatory T Cells. Sci Transl Med (2015) 7(315):315ra189. doi:10.1126/scitranslmed.aad4134

17. Brunstein, CG, Miller, JS, McKenna, DH, Hippen, KL, DeFor, TE, Sumstad, D, et al. Umbilical Cord Blood-Derived T Regulatory Cells to Prevent GVHD: Kinetics, Toxicity Profile, and Clinical Effect. Blood (2016) 127(8):1044–51. doi:10.1182/blood-2015-06-653667

18. Tang, Q, and Vincenti, F. Transplant Trials with Tregs: Perils and Promises. J Clin Invest (2017) 127(7):2505–12. doi:10.1172/jci90598

19. Boardman, DA, Philippeos, C, Fruhwirth, GO, Ibrahim, MAA, Hannen, RF, Cooper, D, et al. Expression of a Chimeric Antigen Receptor Specific for Donor HLA Class I Enhances the Potency of Human Regulatory T Cells in Preventing Human Skin Transplant Rejection. Am J Transpl (2017) 17(4):931–43. doi:10.1111/ajt.14185

20. Romano, M, Fanelli, G, Albany, CJ, Giganti, G, and Lombardi, G. Past, Present, and Future of Regulatory T Cell Therapy in Transplantation and Autoimmunity. Front Immunol (2019) 10:43. doi:10.3389/fimmu.2019.00043

21. Raffin, C, Vo, LT, and Bluestone, JA. Treg Cell-Based Therapies: Challenges and Perspectives. Nat Rev Immunol (2020) 20(3):158–72. doi:10.1038/s41577-019-0232-6

22. Boardman, DA, and Levings, MK. Cancer Immunotherapies Repurposed for Use in Autoimmunity. Nat Biomed Eng (2019) 3(4):259–63. doi:10.1038/s41551-019-0359-6

23. Dawson, NA, Lamarche, C, Hoeppli, RE, Bergqvist, P, Fung, VC, McIver, E, et al. Systematic Testing and Specificity Mapping of Alloantigen-specific Chimeric Antigen Receptors in Regulatory T Cells. JCI Insight (2019) 4(6):e123672. doi:10.1172/jci.insight.123672

24. MacDonald, KG, Hoeppli, RE, Huang, Q, Gillies, J, Luciani, DS, Orban, PC, et al. Alloantigen-specific Regulatory T Cells Generated with a Chimeric Antigen Receptor. J Clin Invest (2016) 126(4):1413–24. doi:10.1172/jci82771

25. Atif, M, Conti, F, Gorochov, G, Oo, YH, and Miyara, M. Regulatory T Cells in Solid Organ Transplantation. Clin Transl Immunol (2020) 9(2):e01099. doi:10.1002/cti2.1099

26. Koyama, I, Bashuda, H, Uchida, K, Seino, K-i., Habu, S, Nakajima, I, et al. A Clinical Trial with Adoptive Transfer of Ex Vivo-induced, Donor-specific Immune-Regulatory Cells in Kidney Transplantation-A Second Report. Transplantation (2020) 104(11):2415–23. doi:10.1097/tp.0000000000003149

27. Altrock, PM, Liu, LL, and Michor, F. The Mathematics of Cancer: Integrating Quantitative Models. Nat Rev Cancer (2015) 15(12):730–45. doi:10.1038/nrc4029

28. Culshaw, RV. Review of HIV Models: The Role of the Natural Immune Response and Implications for Treatment. J Biol Syst (2004) 12(02):123–35. doi:10.1142/s0218339004001099

29. Kirschner, D, Pienaar, E, Marino, S, and Linderman, JJ. A Review of Computational and Mathematical Modeling Contributions to Our Understanding of Mycobacterium tuberculosis Within-Host Infection and Treatment. Curr Opin Syst Biol (2017) 3:170–85. doi:10.1016/j.coisb.2017.05.014

30. Arciero, JC, Maturo, A, Arun, A, Oh, BC, Brandacher, G, and Raimondi, G. Combining Theoretical and Experimental Techniques to Study Murine Heart Transplant Rejection. Front Immunol (2016) 7:448. doi:10.3389/fimmu.2016.00448

31. Gaetano, A, Matone, A, Agnes, A, Palumbo, P, Ria, F, and Magalini, S. Modeling Rejection Immunity. Theor Biol Med Model (2012) 9:18. doi:10.1186/1742-4682-9-18

32. Day, JD, Metes, DM, and Vodovotz, Y. Mathematical Modeling of Early Cellular Innate and Adaptive Immune Responses to Ischemia/Reperfusion Injury and Solid Organ Allotransplantation. Front Immunol (2015) 6:484. doi:10.3389/fimmu.2015.00484

33. Banks, HT, Hu, S, Jang, T, and Kwon, H-D. Modelling and Optimal Control of Immune Response of Renal Transplant Recipients. J Biol Dyn (2012) 6:539–67. doi:10.1080/17513758.2012.655328

34. An, G. Introduction of a Framework for Dynamic Knowledge Representation of the Control Structure of Transplant Immunology: Employing the Power of Abstraction with a Solid Organ Transplant Agent-Based Model. Front Immunol (2015) 6:561. doi:10.3389/fimmu.2015.00561

35. Akinwumi, SM. Modeling the Kinetics of EBV in Primary Carriers and Transplant Recipients. Edmonton, AB: University of Alberta (2017). doi:10.7939/R3HH6CK80

36. Oberbarnscheidt, MH, Zeng, Q, Li, Q, Dai, H, Williams, AL, Shlomchik, WD, et al. Non-self Recognition by Monocytes Initiates Allograft Rejection. J Clin Invest (2014) 124(8):3579–89. doi:10.1172/jci74370

37. Deb, K, Pratap, A, Agarwal, S, and Meyarivan, T. A Fast and Elitist Multiobjective Genetic Algorithm: NSGA-II. IEEE Trans Evol Computat (2002) 6(2):182–97. doi:10.1109/4235.996017

38. Sakaguchi, S, Yamaguchi, T, Nomura, T, and Ono, M. Regulatory T Cells and Immune Tolerance. Cell (2008) 133(5):775–87. doi:10.1016/j.cell.2008.05.009

39. Wood, KJ, Bushell, A, and Hester, J. Regulatory Immune Cells in Transplantation. Nat Rev Immunol (2012) 12(6):417–30. doi:10.1038/nri3227

40. Bontha, SV, Maluf, DG, Mueller, TF, and Mas, VR. Systems Biology in Kidney Transplantation: The Application of Multi-Omics to a Complex Model. Am J Transpl (2017) 17(1):11–21. doi:10.1111/ajt.13881

41. Mastoridis, S, Martínez-Llordella, M, and Sanchez-Fueyo, A. Emergent Transcriptomic Technologies and Their Role in the Discovery of Biomarkers of Liver Transplant Tolerance. Front Immunol (2015) 6:304. doi:10.3389/fimmu.2015.00304

42. Stegall, MD, and Borrows, R. Computational Biology: Modeling Chronic Renal Allograft Injury. Front Immunol (2015) 6:385. doi:10.3389/fimmu.2015.00385

43. Yang, JYC, and Sarwal, MM. Transplant Genetics and Genomics. Nat Rev Genet (2017) 18(5):309–26. doi:10.1038/nrg.2017.12

44. Issa, F, Hester, J, Milward, K, and Wood, KJ. Homing of Regulatory T Cells to Human Skin Is Important for the Prevention of Alloimmune-Mediated Pathology in an In Vivo Cellular Therapy Model. PLoS One (2012) 7(12):e53331. doi:10.1371/journal.pone.0053331

45. Zhang, N, Schröppel, B, Lal, G, Jakubzick, C, Mao, X, Chen, D, et al. Regulatory T Cells Sequentially Migrate from Inflamed Tissues to Draining Lymph Nodes to Suppress the Alloimmune Response. Immunity (2009) 30(3):458–69. doi:10.1016/j.immuni.2008.12.022

46. Campbell, DJ. Control of Regulatory T Cell Migration, Function, and Homeostasis. J.Immunol. (2015) 195(6):2507–13. doi:10.4049/jimmunol.1500801

47. Issa, F, Robb, RJ, and Wood, KJ. The where and when of T Cell Regulation in Transplantation. Trends Immunol (2013) 34(3):107–13. doi:10.1016/j.it.2012.11.003

48. Jhunjhunwala, S, Chen, LC, Nichols, EE, Thomson, AW, Raimondi, G, and Little, SR. All-trans Retinoic Acid and Rapamycin Synergize with Transforming Growth Factor-Β1 to Induce Regulatory T Cells but Confer Different Migratory Capacities. J Leukoc Biol (2013) 94(5):981–9. doi:10.1189/jlb.0312167

49. MacDonald, KG, Orban, PC, and Levings, MK. T Regulatory Cell Therapy in Transplantation. Curr Opin Organ Transplant (2012) 17(4):343–8. doi:10.1097/mot.0b013e328355aaaf

50. Camirand, G, and Riella, LV. Treg‐Centric View of Immunosuppressive Drugs in Transplantation: A Balancing Act. Am J Transpl (2017) 17(3):601–10. doi:10.1111/ajt.14029

51. Furukawa, A, Wisel, SA, and Tang, Q. Impact of Immune-Modulatory Drugs on Regulatory T Cell. Transplantation (2016) 100(11):2288–300. doi:10.1097/tp.0000000000001379

52. Eskandari, SK, Sulkaj, I, Melo, MB, Li, N, Allos, H, Alhaddad, JB, et al. Regulatory T Cells Engineered with TCR Signaling-Responsive IL-2 Nanogels Suppress Alloimmunity in Sites of Antigen Encounter. Sci Transl Med (2020) 12(569):eaaw4744. doi:10.1126/scitranslmed.aaw4744

53. Fisher, JD, Zhang, W, Balmert, SC, Aral, AM, Acharya, AP, Kulahci, Y, et al. In Situ recruitment of Regulatory T Cells Promotes Donor-specific Tolerance in Vascularized Composite Allotransplantation. Sci Adv (2020) 6(11):eaax8429. doi:10.1126/sciadv.aax8429

54. Brown, LV, Gaffney, EA, Ager, A, Wagg, J, and Coles, MC. Quantifying the Limits of CAR T-Cell Delivery in Mice and Men. J R Soc Interf (2021) 18(176):20201013. doi:10.1098/rsif.2020.1013

55. Alexander, HK, and Wahl, LM. Self-tolerance and Autoimmunity in a Regulatory T Cell Model. Bull Math Biol (2011) 73(1):33–71. doi:10.1007/s11538-010-9519-2

56. Fouchet, D, and Regoes, R. A Population Dynamics Analysis of the Interaction between Adaptive Regulatory T Cells and Antigen Presenting Cells. Plos One (2008) 3(5):e2306. doi:10.1371/journal.pone.0002306

57. Garcia-Martinez, K, and Leon, K. Modeling the Role of IL2 in the Interplay between CD4+ Helper and Regulatory T Cells: Studying the Impact of IL2 Modulation Therapies. Int Immunol (2012) 24(7):427–46. doi:10.1093/intimm/dxr120

58. León, K, Lage, A, and Carneiro, J. Tolerance and Immunity in a Mathematical Model of T-Cell Mediated Suppression. J Theor Biol (2003) 225(1):107–26. doi:10.1016/s0022-5193(03)00226-1

59. Sicard, A, Lamarche, C, Speck, M, Wong, M, Rosado‐Sánchez, I, Blois, M, et al. Donor‐specific Chimeric Antigen Receptor Tregs Limit Rejection in Naive but Not Sensitized Allograft Recipients. Am J Transpl (2020) 20(6):1562–73. doi:10.1111/ajt.15787

60. Golshayan, D, Jiang, S, Tsang, J, Garin, MI, Mottet, C, and Lechler, RI. In Vitro-expanded Donor Alloantigen-specific CD4+CD25+ Regulatory T Cells Promote Experimental Transplantation Tolerance. Blood (2007) 109(2):827–35. doi:10.1182/blood-2006-05-025460

61. Joffre, O, Santolaria, T, Calise, D, Saati, TA, Hudrisier, D, Romagnoli, P, et al. Prevention of Acute and Chronic Allograft Rejection with CD4+CD25+Foxp3+ Regulatory T Lymphocytes. Nat Med (2008) 14(1):88–92. doi:10.1038/nm1688

Keywords: regulatory T cells, rejection, heart transplant, mathematical model, adoptive transfer, immune response

Citation: Lapp MM, Lin G, Komin A, Andrews L, Knudson M, Mossman L, Raimondi G and Arciero JC (2022) Modeling the Potential of Treg-Based Therapies for Transplant Rejection: Effect of Dose, Timing, and Accumulation Site. Transpl Int 35:10297. doi: 10.3389/ti.2022.10297

Received: 10 December 2021; Accepted: 17 March 2022;
Published: 11 April 2022.

Copyright © 2022 Lapp, Lin, Komin, Andrews, Knudson, Mossman, Raimondi and Arciero. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Giorgio Raimondi, g.raimondi@jhmi.edu; Julia C. Arciero, jarciero@iupui.edu