Front. Immunol.Frontiers in ImmunologyFront. Immunol.1664-3224Frontiers Media S.A.10.3389/fimmu.2024.1463915ImmunologyOriginal ResearchMechanistic PKPD modeling to describe cytokine release associated with CD3 T-cell engager therapiesLefèvreApolline12*Parra-GuillenZinnia P.23TrocónizIñaki F.234BoetschChristophe1FrancesNicolas11Roche Pharma Research and Early Development (pRED), Pharmaceutical Sciences PS, Roche Innovation Center Basel, Basel, Switzerland2Pharmacometrics & Systems Pharmacology, Department of Pharmaceutical Science, School of Pharmacy and Nutrition, University of Navarra, Pamplona, Spain3Navarra Institute for Health Research (IdiSNA), Pamplona, Spain4Institute of Data Science and Artificial Intelligence (DATAI), University of Navarra, Pamplona, Spain
Edited by: Jose A. Garcia-Sanz, Spanish National Research Council (CSIC), Spain
Reviewed by: Yanguang Cao, University of North Carolina at Chapel Hill, United States
Luis Alberto Vallejo-Castillo, National Polytechnic Institute (IPN), Mexico
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.
Introduction
T-cell engagers (TCE), a therapeutic class of cancer immunotherapy (CIT), offer a novel approach to cancer treatment by harnessing and reactivating the patient’s immune system to eradicate tumor cells. However, the use of TCE in the clinic can lead to severe side effects, including cytokine release syndrome (CRS). Therefore, innovative dosing strategies need to be implemented to mitigate the risk of developing CRS.
Method
In the presented work, a mechanistic pharmacokinetics/pharmacodynamics (PKPD) model describing cytokine release following TCE therapy has been developed combining literature knowledge and preclinical data. The model was developed to explore and test hypotheses regarding the mechanisms behind the decrease of cytokine release following two repeated TCE administrations.
Results
The model is able to successfully reproduce the observed dynamics of cytokine levels associated with the initial and subsequent TCE doses, accounting for different dosing intervals. In addition, the model suggests a mechanism of action that uncouples cytokine release from tumor cell killing.
Discussion
This model provides an initial mechanistic framework to support the design of experiments and paves the way for the application of mathematical modeling to support clinical dosing regimen selection of any TCE.
cancer immunotherapy (CIT)T-cell engager (TCE)cytokine releasepharmacokinetics/pharmacodynamics (PK/PD) modelingeffector T-cell dynamicssection-in-acceptanceCancer Immunity and ImmunotherapyIntroduction
Cancer immunotherapy (CIT) has recently drawn remarkable attention by showing its ability to improve the overall survival in the clinic (1) and has led to the approval of more than 40 drugs by health authorities (2). Various classes of therapeutic agents exist in this area and share the concept of enhancing the existing anti-tumor immunity to eradicate malignant cells (3). One approach is to use synthetic immunity by redirecting the T-cells against the tumor cells (1). T-cell engagers (TCE) are bispecific antibodies (bsAbs) (4) that bind simultaneously to T-cells via the CD3 epsilon chain (CD3) within the T-cell Receptor (TCR) complex and to tumor cells via tumor-associated antigen (TAA) (5). The crosslinking of a tumor cell with a T lymphocyte enables the formation of an immunological synapse, similar to that of a natural TCR and peptide–major histocompatibility complex (MHC) complex (6). The synapse formation redirects any engaged T-cells in an inflammatory process (cytokine release, triggering further immune activity) and tumor killing.
TCEs are thus promising anti-tumor therapies but have remaining challenges associated with their clinical safety profiles. The primary toxicities following their administration are immune effector cell-associated neurotoxicity syndrome (ICANS) and cytokine release syndrome (CRS) (7). The latter is generally associated with cytokine release upon T-cell activation (8). Patients experiencing CRS may have a variety of symptoms (9) for which the grading has been defined by the Common Terminology Criteria for Adverse Events (CTC-AE) v5 (10). CRS incidence (grade I or more) in clinical trials of FDA approved bispecific TCEs ranged from 15% to 89% (11).
Strategies are implemented in the clinic to mitigate CRS, such as step-up dosing (SUD) (12), which involves gradually increasing the dose within the same patient until reaching the target dose; the use of a cytokine blockage agent (e.g. Tocilizumab (13), an anti-IL6 antibody) (14) or the use of corticosteroids (e.g. Dexamethasone) (15). In addition, in vitro data suggest that SUD can be used as a prophylactic method to prevent CRS in the clinical setting (16). Selection of the best dosing regimen, particularly in the context of SUD, is currently based on empirical methods. Mathematical models, including the proposed mechanistic PKPD model, can be used to inform dosing regimen selection and have the potential to reduce both the length of dose escalation trials and the empirical explorations of dose and schedule.
Some mathematical models predicting cytokine release after TCE administration have been previously published (17–19). Those models describe the decrease in cytokine release observed upon repeated TCE administration as a consequence of previous target cell depletion. The model from Chen and colleagues (19) additionally assumes that cytokines, once released, inhibit the release of new cytokines via a negative feedback loop. The hypothesis that the loss of cytokine release is a consequence of the T-cell killing of target cells was refuted by the work of Li et al. (16). They demonstrated in vitro that, when a second TCE stimulation was done seven days after the first one, lower levels of cytokines were associated with the second dose versus the first one while cytotoxic activity was conserved. This was observed despite the in vitro system being re-supplemented with fresh tumor cells prior to the second stimulation. Alternatively, the pool model proposed by Movin-Osswald and Hammarlund-Udenaes (20) could provide another hypothesis for the reduction in cytokine release following two TCE administrations. The original pool model assumes the existence of an initial reservoir (a pool) of available prolactin (PRL) that is depleted upon remoxipride administration, thus preventing subsequent releases upon multiple dosing until the pool has been replenished. Although this approach has not been previously applied within the context of CRS, it could be adapted by considering the cytokines as a similar pool.
Hereby, a hypothesis-driven mechanistic PKPD model is presented with the aim of exploring potential mechanistic insights behind cytokine release upon the two first TCE administrations. The proposed PKPD model expands upon previously published models (17–20) by incorporating evidence from newly available published preclinical data (16). The presented model focuses on predicting cytokine release upon the first two consecutive TCE administrations with varying dosing intervals and uncouples cytokine levels from tumor burden.
Materials and methodsPreclinical literature data
The PKPD model was developed based on literature preclinical data from Li et al. (16) where mechanisms of cytokine release were investigated. The group performed an in vivo study in which mammary tumor-bearing mice were treated with two subsequent doses of anti-HER2 T-cell Bispecific (0.5mg/kg) at different dosing intervals (administrations at day 0 and day 1; day 0 and day 7; day 0 and day 14; day 0 and day 21 or day 0 and day 28). Systemic cytokine release (interleukin-6, IL6) was monitored two hours after each treatment. Published data from Li et al. (16) were digitized using WebPlotDigitizer (21) and are available in the appendix section (Supplementary Table S1).
PKPD model development
A mechanistic PKPD model was built using a bottom-up approach including the main biological players involved in cytokine release: TCE, tumor cells, T-cells and interleukins. The model (Figure 1) assumed that upon TCE exposure, T-cells are activated and transition through different states in which they can trigger cytokine release and kill tumor cells, as supported by (22). The dynamics of the different entities were modeled using ordinary differential equations (ODE) describing the mass balance between the synthesis, transition and degradation or apoptosis processes and implemented using zero-, first- or second-order rate constants. Michaelis-Menten or Hill kinetic functions were also implemented to account for saturation processes. Given the limited available pre-clinical data to build the model, parameters could not be estimated. Model parameters characterizing the above mentioned processes were taken from literature or fine tuned (calibrated) to describe the available data of in vivo IL6 (16). Some parameters had to be borrowed from human references to simulate the PKPD model in mice.
Graphical representation of the mechanistic pharmacokinetic/pharmacodynamic (PK/PD) model. The model describes the T-cell engager (TCE) plasma concentration (Cpla) and its impact on the transition from naive T-cells (Tnaive) to activated T-cells (Tact), leading to cytokine release (IL6) and tumor cell killing. Biological entities and parameters are described in the text.
Model evaluation and exploration
Model performance was graphically assessed by confronting model predictions to the in vivo digitized data from (16). The unique IL6 sample (two hours following TCE administration) was assumed to be reflective of the maximum IL6 levels (23). In addition, the cytokine release and killing functions were explored under different relevant dosing scenarios to understand the behavior of the system (same dose levels, 0.5mg/kg, and different dosing intervals). Predictions from the PKPD model were also compared to simulations derived from a simpler model, the pool model from (20), which was repurposed to describe cytokine release (Appendix). Lastly, a parameter sensitivity analysis was performed to evaluate the effects of varying each parameter individually (+/- 20%) (24) on the ratio of maximum cytokine release at second TCE administration versus first administration, and identify the most influential parameters. Parameters involved in tumor cell killing were not included in the sensitivity analysis as they have no impact on cytokine release. Model simulations (ODE computations) were conducted in Rstudio version 4.1.3 (25) using the RxODE2 package (26).
ResultsModel description
The structure of the developed mechanistic PKPD model is presented in Figure 1.
Pharmacokinetics
The TCE plasma concentration (Cpla) was assumed to follow a 1-compartment PK model characterized by a linear clearance (CL) and distribution volume (V) parameters from a similar anti-HER2 TCE (27) (Equation 1). Initial drug concentration prior to dosing, Cpla(t=0), is equal to zero.
dCpladt=−CLV×CplaCpla(t=0)=0T-cell dynamics
Three T-cell states were included in the model: baseline naive T-cells (Tnaive), activated T-cells (Tact) and desensitized T-cells (Tdesens). Baseline naive T-cells (Equation 2) are activated upon TCE exposure as follow:
where kin,naive and kapop are the zero and first-order rate constants for baseline naive T-cells input and apoptosis, respectively. At baseline, the amount of naive T-cells Tnaive(t=0) is given by kin,naive/kapop. The activation signal, fKact, described in Equation 3 follows an Emax -like relationship with respect to the TCE plasma concentration, characterized by a Hill coefficient hact, a maximum activation signal Emax,act and an EC50act, which is the TCE concentration required to reach half of maximal T-cell activation. As only one dose level (0.5mg/kg) was tested in vivo, the model assumes a complete transfer of the naive T-cells to activated T-cells at the tested dose. The parameter values should be adjusted with more evidence from new experiments.
fKact=Emax,act×(CplahactCplahact+EC50acthact)
The transition from baseline naive T-cells to activated T-cells state (Tactdelay) was modeled with 5 transit compartments (Tactdelay1 to Tactdelay5) to account for the lag time between maximum TCE exposure and maximum cytokine level, and is governed by Equations 4, 5. At baseline, no activated T-cells are present and Tactdelay(t=0) = 0. The parameter value for the first-order transit rate constant kdelay between the transit compartments was selected to mimic a delay of a few hours between drug administration and cytokine release (16). Activated T-cells (Tact) undergo apoptosis at the same rate parameter as naive T-cells, kapop (see below).
The proliferation of T-cells was assumed to be associated with an IL2 cytokine release upon TCE activation, as IL2 secreted by T-cells binds to T-cell receptors and induces their clonal proliferation (28). An IL2effect variable was therefore introduced (Equations 4, 5) mimicking an autocrine loop for T-cell proliferation in presence of IL2 as follow:
where kIL2 is the maximum proliferation rate parameter of T-cells induced by IL2. In turn, this process can be inhibited when effector cells (T-cells) outnumber target cells (29). Hereby, T-cell proliferation by IL2 (Equation 6) is inhibited by a maximum inhibition rate parameter (Imax), an EC50E:T,IL2 parameter which correspond to the Effector to Target (E:T) ratio parameter needed to achieve half of maximum inhibition IL2 effect and a shape parameter (hIL2). All T-cells (naive, activated and desensitized) as well as all tumor cells are accounted for in the E:T ratio. For this analysis and in absence of clinical information, the initial values of T-cells and tumor cells were selected to match an optimal in vitro E:T ratio of 10 (16).
Cytokine release
Only activated T-cells from the fifth T-cells activated transit compartment (Tactdelay5) release cytokine (IL6), Cyt, (Equation 7) according to the following equation:
in which kin,cyt and kdeg,cyt are zero and first-order rate constants representing the physiological production and degradation of cytokines, respectively. Baseline cytokine levels, Cyt(t=0), are set to 10pg/mL (30) and defined by kin,cyt/kdeg,cyt. Function fKCyt (Equation 8) describes the increase of cytokine production upon T-cell activation as follow:
Where Emax,cyt is the maximum increase in cytokine production, hcyt is the Hill coefficient representing the steepness of the exposure-cytokine release relationship, and EC50cyt is the TCE concentration required to reach half of maximal cytokine production increase. In Equation 8, cytokine release also depends on the amount of Tactdelay5 and the shape parameter α. Those parameters were calibrated to describe the in vivo data from Li et al. (16).
Subsequently, the activated T-cells transit into a desensitized T-cell status, Tdesens, which represent T-cells that retained their ability to kill tumor cells without releasing cytokines (Equation 9). This desensitization process might be driven by the release of cytokine itself. As an example, the release of IL10, a regulatory cytokine, is thought to inhibit the release of more pro-inflammatory cytokines like IL6 and TNFα (31). Baseline levels of desensitized T-cells, Tdesens(t=0), are set to zero. The same parameter value (kdelay) has been assumed for the desensitization rate, in absence of data to differentiate them.
It has been hypothesized that desensitized T-cells are, as activated T-cells, affected by the T-cell proliferation effect (IL2effect) (Equation 6).
Tumor cell killing
Li et al. (16) demonstrated in vitro that reduced cytokine release upon the second dose was independent of tumor cell depletion and did not compromise cytotoxic efficacy. Tumor cell growth and killing were included in the model to determine if the observed in vitro data could be explained by the hereby proposed mechanisms. Tumor cells (Tumor) were assumed to follow an exponential growth (Equation 10) with a growth rate kg. They are killed by activated (Tactdelay1, …) and desensitized (Tdesens) T-cells as determined by the fKkill killing function (Equation 14) following an Emax relationship versus drug concentration (Cpla, EC50kill) and on the number of T-cells available for killing (Tkill, Equation 15). It was assumed that the TCE was necessary to modulate the activation of the T-cells and thus, no baseline killing was described when no drug was administered. Tumor cell killing was modeled with 3-transit compartments (Equations 11-13, Tumor1 to Tumor3) with a delay of kdelay,tumor (32).
As no tumor data was available from the in vivo experiments, growth rate (kg) and transit delay (kdelay,tumor) were taken from the literature (32) while plausible values were selected for the Emax,kill and EC50kill.
The model parameters of the mechanistic PKPD model are displayed in Table 1.
Mechanistic PKPD model parameters (in mince) and their corresponding values, units and source (literature or calibrated with the available in vivo data).
Parameter
Value
Units
Source/method
Description
PK
V
2.3
mL
Adapted fromYu et al. (27)
Volume of the central compartment
CL
0.025
mL/h
Elimination rate constant
T-cell kinetics
kin,naive
16
cells/uL/h
Lythe et al. (38)(theoretical computations)
Synthesis rate of naive T-cells *
kapop
0.01
1/h
Adapted fromFriberg et al. (39)
Apoptosis rate of naive T-cells *
T-cell activation
Emax,act
1.5
1/h
Calibrated
Maximum of T-cell activation signal
hact
2
–
Calibrated
Hill coefficient for activation signal
EC50act
1
µg/mL
Calibrated
Exposure of TCE to achieve half of maximal T-cell activation
kdelay
2.5
1/h
Adapted from Li et al. (16), Friberg et al. (39)
Transit rate between compartments
Cytokine release
kin,cyt
10* kdeg,cyt
1/h
Frances et al.(30)
Synthesis rate of cytokines **
kdeg,cyt
0.41
1/h
Degradation rate of cytokines **
Emax,cyt
0.00086
pg/mL/h
Calibrated
Maximum cytokine release from T-cells
hcyt
2.5
–
Calibrated
Hill coefficient for cytokine release signal
EC50cyt
5
µg/mL
Calibrated
Exposure of TCE to achieve half of maximal cytokine release from T-cells
α
2
–
Calibrated
Exponent
T-cell proliferation by IL2
kIL2
0.0008
1/h
Macallan et al.(40)
Maximum proliferation rate of T-cells (activated and desensitized) induced by IL-2 *
Imax
1
–
Fixed
Maximum inhibition of IL2 effect
hIL2
5
–
Fixed
Hill coefficient for IL2 effect
EC50E:T,IL2
7
–
Adapted from Bacac et. al(5)
E:T to achieve half of maximal IL2 effecta
Tumor kinetics
kg
0.006
1/h
Simeoni et al.(32)
Tumor growth rate
kdelay,tumor
0.02
1/h
Transit rate between tumor compartments
Emax,kill
0.000025
1/cells/uL/h
Calibrated
Maximum killing of tumor cells
EC50kill
1
µg/mL
Adapted from Staflin et al. (23)
TCE plasma concentration needed to achieve half of maximal tumor cell killing for the TCE concentration
In vitro conditions for maximal cytotoxic effect, *From humans, **From cynomolgus monkeys.
Model evaluation
The developed model was used to simulate the time course of IL6 cytokine as in the experimental design of Li et al., i.e. following two consecutive TCE administrations given at different dosing intervals (administrations at day 0 and day 1; day 0 and day 7; day 0 and day 14; day 0 and day 21 or day 0 and day 28). Observed IL6 levels versus PKPD model predictions are shown in Figure 2A. The mechanistic PKPD model was capable of reproducing experimental data (16) with (1) a lower cytokine release after the second compared to the TCE first administration, (2) an increase in cytokine release at the 2nd TCE administration when increasing the dosing interval between the two administrations (3) and a similar cytokine levels between the first and second TCE administration when dosing interval is of 28 days. The observed versus predicted plot (Figure 2B) did not reveal any systematic bias. The model was also able to reproduce the dynamic of other cytokines such as TNFα (Supplementary Figure S1, Appendix), following repeated administration of TCE.
Mechanistic PKPD model performance (A) Model predicted IL6 time course in mice (solid lines) following two consecutive TCE administrations at different dosing intervals (administrations at day 0 and day 1; day 0 and day 7; day 0 and day 14; day 0 and day 21 or day 0 and day 28). Dots represent the digitized experimental data from Li et al. (2019) along with their standard error of mean (SEM). (B) Observed versus model predicted IL6 cytokine levels.
The model was then used to explore the behavior of the different model entities under different scenarios (Figure 3). In the proposed PKPD model, it is assumed that IL6 is released only in presence of naive T-cells, and that the amount released is proportional to the quantity of naive T-cells available at the time of the second TCE administration. Accordingly, the model posits that as almost no naive T-cells are available 7 days after a first TCE administration, if a second dose were to be administered, no additional increments in activated T-cells would be observed, neither in IL6 release. Tumor cell killing would be maintained by the presence of desensitized T-cells (Figure 3B). On the other hand, 21 days after the first TCE administration, the PKPD model assumes that naive T-cells levels are close to their baseline values, which leads to IL6 levels similar to those observed following the initial TCE administration if a second TCE stimulation occurs on that day (Figure 3B).
Mechanistic PKPD model simulations of (A) T-cell dynamics following one TCE administration on day 0. (B) Evolution of the killing function (fKkill) and cytokine release function (fKCyt) over time with two TCE administrations 7 (upper panel) or 21 (lower panel) days apart. The gray arrows represent the TCE dosing events (same dose).
The results of the parameter sensitivity analysis for a dosing interval of 21 days are presented in Figure 4. The analysis revealed that the parameters with highest influence on IL6 release predictions are the rate constants associated with naive T-cells (kin,naive and kapo,naive), as well as the parameters associated with the cytokine release function in relation to plasma concentration (α, EmaxCyt and EC50Cyt). The same parameters were highlighted by the sensitivity analysis that was performed for the other tested dosing intervals. The alpha (α) parameter influences the dynamics of the cytokine release (eg, return to baseline after Cmax). Its value could be informed with a dataset capturing the dynamics of cytokines after Cmax.
Parameter sensitivity analysis: Impact of a sequential 20% deviation on each parameter on the simulated ratio of IL6 maximum concentration at first administration versus second administration 21 days later. Parameters influencing tumor cell killing were omitted (not affecting cytokine release).
Discussion
TCEs have shown promising potential in cancer therapy, but their clinical development can be risky for patients due to their propensity of inducing cytokine release, a potentially life-threatening condition. Cytokine blockage agents like Tocilizumab are often used as a rescue method for high grades CRS. Conservative approaches like the in vitro minimum anticipated biological effect level (MABEL) approach (8) are also recommended to prevent high grades of CRS for the first-in-human (FIH) dose selection. Nevertheless, the in vitro MABEL derived FIH dose can be substantially lower than the therapeutic dose, leading to long dose escalation phases (12). Consequently, with this approach, many patients are treated at sub-therapeutic doses while some others are put at risk when approaching the Maximum Tolerated Dose (MTD). The final recommended dose and schedule for the drug label will be empirically decided based on all evidence gathered during the clinical development phase. The occurrence of severe CRS in the clinic can also be mitigated by using a SUD approach. Nevertheless, the determination of the optimal SUD dosing regimen remains a major challenge and is often the result of an empirical approach involving many patients. Hereby, a mathematical model capable of predicting cytokine release during the first two TCE administrations is presented. It brings quantitative insights on the dosing interval for the step-up dosing procedure that could be used for any CD3 TCE. The mechanistic PKPD model can be eligible for a Model-Informed Drug Development (MIDD) approach to prospectively design clinical studies (33). It can serve to inform dosing regimen to minimize cytokine release, potentially reducing the CRS risk for the patients and hopefully shortening the early phases of clinical development. It can also be used in a sequential approach to refine the next cohort dosing regimen based on emerging clinical data, similar to the methodology used for the Escalation With Overdose Control (EWOC) clinical trial design (34).
The presented mechanistic PKPD model was expanded from previously published models (17–19) to integrate newly available data (16). In the absence of PK data from this study, the PKPD model uses a simple one-compartment PK model to capture TCE concentration. This choice provides a straightforward framework that can be adjusted based on the specific properties of the studied molecule such as Target-Mediated-Drug-Disposition (TMDD) or FcRn recycling (35). The mechanistic PKPD model describes the main findings from Li et al. (decreased cytokine release upon second TCE administration) as a shift in the T-cell phenotype, and as opposed to the loss of tumor cells suggested in previous models (17, 18). The proposed mechanistic PKPD model incorporates three T-cell states as in Hosseini and colleagues (18), with the difference that upon dosing, the T-cells move to a desensitized state in which they can no longer release cytokines but retain their cytotoxic potential (as suggested by (16, 36)). Here, it is assumed that there is an initial limited pool of naive T-cells available for activation (and therefore available for cytokine release) which takes time to fully replenish and prevent cytokine release upon following TCE administrations. Upon first TCE stimulation, the pool of baseline naive T-cells transition to an activated T-cell status able to release cytokines. Those T-cells transit then to a status referred to as desensitized, by losing their ability to release cytokines. After this first TCE stimulation, it takes about 28 days to fully replenish the pool of baseline T-cells (in absence of any new stimulation), and retrieve the cytokine levels following the first TCE administration. In the meantime, any new TCE stimulation will lead to less cytokine release compared to the first TCE stimulation, as the number of naive T-cells is also lower. A shorter dosing interval, as opposed to a longer one, will result in reduced systemic release of cytokines. The model simulation behaves similarly to the clinical observations where a step-up dosing approach (a priming low dose is given followed by a higher dose within a short dosing interval) is taken. Cytokine levels are reduced after the second dose, despite it being higher than the first one (12, 37). With one dose level tested in vivo (0.5mg/kg), it was not possible to establish a relationship between drug exposure, T-cell activation and consequently cytokine release. Thus, the model assumes a complete transfer of the naive T-cells to activated T-cells at the tested dose. To allow for an increase in cytokine release, the release function (fKcyt) depends thus on both the amount of available T-cells and the drug concentration. The principle of an initial reservoir (pool) of biological entities to be emptied by drug exposure was first proposed by Movin-Osswald and Hammarlund-Udenaes (20) in a model referred to as the “pool model” (Appendix). In the pool model, the group assumed that an initial pool of prolactin (PRL) is depleted following drug exposure and takes time to replenish. The original pool model was repurposed to describe cytokine release (Supplementary Figure S2, Appendix) and its performance was compared to the presented PKPD model. Both models performed equally well (Supplementary Figure S3, Appendix). Nevertheless, the newly proposed PKPD model provides a new interpretation of the data and integrates new components considered essential for further development and use of the model, such as the tumor response that can be uncoupled from cytokine release upon second TCE administration, as suggested in vitro (16). T-cell proliferation is also included in the mechanistic PKPD model as it is an essential element to describe the individual patient’s response in the clinic (extent of tumor killing, delay in response…). An empirical approach driven by the E:T ratio was used to describe the proliferation effect, which is assumed to be driven by IL2 (IL2effect) but can in reality integrate indirectly other mechanisms. At high ratios (above 10), the proliferation effect is inhibited to avoid an infinite T-cell proliferation. Various processes like IL2 downregulation or regulatory T-cells may explain this inhibition. The choice of the E:T ratio was driven by in vitro knowledge since the in vivo E:T ratio was unknown. It must be noted that this may not reflect the in vivo and/or clinical conditions. With the proposed mechanistic PKPD model, T-cells count and status could be used as biomarkers to inform on individual dosing regimen, allowing for personalized dosing.
The presented model was developed with the goal of providing a theoretical framework for understanding cytokine release dynamics upon the first two TCE administrations and generating hypotheses, ultimately contributing to the optimization of TCE dosing regimens. The PKPD model suggests biological mechanisms that can be tested in preclinical or clinical settings. In particular, in vivo longitudinal studies of cytokine release, T-cell state and evolution, and tumor cell killing for different dosing regimens (dose levels, number of doses, dosing intervals) would better characterize the relationship between T-cell number/state and cytokine release. Consequently, this would allow to better inform the activation function of the T-cells (fKact), potentially removing the dependency on drug plasma concentration for cytokine release. The model was developed based on in vitro and in vivo data, and clinical data would be needed to confirm the model predictions. In particular, testing different dose levels would also enable the refinement of the relationship between cytokine release and TCE exposure after the first dose. Furthermore, due to the true complexity of the immune system and of the mechanisms involved in cytokine release, as well as the lack of available data, some simplifications were made in the model. For instance, the cellular source of cytokines can be diverse and the presented model considers the cytokine release from T-cells only. As suggested in (15, 22), monocytes or macrophages may become activated following T-cell activation triggered by TCE, potentially contributing to additional cytokine release. In addition, the parameters for T-cell kinetics may differ from reality. Specifically, the assumption of identical rates for activation and desensitization (kdelay) may not accurately represent the underlying biological processes, and some parameters (kin,naive, kapop) were derived from human data while simulating cytokine release in mice. Moreover, the model could be updated with more data to link the cytokine release to the trimeric complex formation (TCE, T-cell, tumor cell), as in the model from Jiang et al. (17). This would better explain the variability observed in the patient population potentially due to the differences in tumor burden and T-cell levels.
Conclusion
A new mechanistic PKPD model has been developed to describe IL6 cytokine levels associated with the first and second TCE doses with varying dosing intervals between the two administrations. In the proposed model, the decrease of cytokine release upon second TCE administration observed in Li et al. (16) is hypothesized to result from a transition of T-cells, and not from the loss of tumor cells. The shift takes the T-cells from a state capable of cytokine release to another state where they are not capable of such release. In particular, the model is able to link the second administration cytokine levels to the interval versus the first TCE administration. The hypotheses stated in the model are thought to be valid for any CD3-TCE and could thus be used for the drug development of multiple molecules. The model suggests refined pre-clinical experiments to further support drug development by providing more mechanistic insights on the TCE associated cytokine release. It also opens opportunities to leverage PKPD modeling during the clinical development, potentially increasing the likelihood of success of developing new TCE drugs, by moving quicker in optimized drug dosing schemes, maximizing TCE associated patient’s benefits while minimizing the risk of CRS.
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 author.
The author(s) declare that no financial support was received for the research, authorship, and/or publication of this article. ZP-G and IFT received funding as part of the collaboration for PhD supervision.
Acknowledgments
Figure 1 was created with BioRender.com.
Conflict of interest
CB, NF, and AL are employed by F. Hoffmann-La Roche Ltd, Switzerland. CB and NF own Roche stocks.
The remaining 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.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fimmu.2024.1463915/full#supplementary-material
ReferencesMaHWangHSoveRJJafarnejadMTsaiCHWangJ. A quantitative systems pharmacology model of T cell engager applied to solid tumor. AAPS J. (2020) 22:85. doi: 10.1208/s12248-020-00450-3Sancho-AraizAMangas-SanjuanVTroconizIF. The role of mathematical models in immuno-oncology: challenges and future perspectives. Pharmaceutics. (2021) 13(7):1016. doi: 10.3390/pharmaceutics13071016ShimuASWeiHXLiQZhengXLiB. The new progress in cancer immunotherapy. Clin Exp Med. (2023) 23:553–67. doi: 10.1007/s10238-022-00887-0ShanshalMCaimiPFAdjeiAAMaWW. T-cell engagers in solid cancers-current landscape and future directions. Cancers (Basel). (2023) 15(10):2824. doi: 10.3390/cancers15102824BacacMFautiTSamJColombettiSWeinzierlTOuaretD. A novel carcinoembryonic antigen T-cell bispecific antibody (CEA TCB) for the treatment of solid tumors. Clin Cancer Res. (2016) 22:3286–97. doi: 10.1158/1078-0432.CCR-15-1696MiddelburgJKemperKEngelbertsPLabrijnAFSchuurmanJvan HallT. Overcoming challenges for CD3-bispecific antibody therapy in solid tumors. Cancers (Basel). (2021) 13(2):287. doi: 10.3390/cancers13020287Godbersen-PalmerCCoupetTAGradaZZhangSCSentmanCL. Toxicity induced by a bispecific T cell-redirecting protein is mediated by both T cells and myeloid cells in immunocompetent mice. J Immunol. (2020) 204:2973–83. doi: 10.4049/jimmunol.1901401SaberHDel VallePRicksTKLeightonJK. An FDA oncology analysis of CD3 bispecific constructs and first-in-human dose selection. Regul Toxicol Pharmacol. (2017) 90:144–52. doi: 10.1016/j.yrtph.2017.09.001Shimabukuro-VornhagenAGodelPSubkleweMStemmlerHJSchlosserHASchlaakM. Cytokine release syndrome. J Immunother Cancer. (2018) 6:56. doi: 10.1186/s40425-018-0343-9U.S. Department of Health and Human ServicesNational Institutes of HealthNational Cancer Institute. Common Terminology Criteria for Adverse Events (CTCAE) Version 5.0. (2017).ShahDSoperBShoplandL. Cytokine release syndrome and cancer immunotherapies - historical challenges and promising futures. Front Immunol. (2023) 14:1190379. doi: 10.3389/fimmu.2023.1190379BallKDovediSJVajjahPPhippsA. Strategies for clinical dose optimization of T cell-engaging therapies in oncology. MAbs. (2023) 15:2181016. doi: 10.1080/19420862.2023.2181016MiharaMKasutaniKOkazakiMNakamuraAKawaiSSugimotoM. Tocilizumab inhibits signal transduction mediated by both mIL-6R and sIL-6R, but not by the receptors of other members of IL-6 cytokine family. Int Immunopharmacol. (2005) 5:1731–40. doi: 10.1016/j.intimp.2005.05.010PharmCJTurnerDCGibianskyLLiFFreyNMenuetJ-F. New tocilizumab (TCZ) dosing guidance for T-cell engaging bispecific antibody-related cytokine release syndrome (CRS) in patients with relapsed/refractory (R/R) B-cell non-hodgkin lymphoma (B-NHL): insights from pooled clinical trial safety experience and quantitative clinical pharmacology (qCP) analyses. Poster presented at American Society of Hematology (ASH). (2023). doi: 10.1182/blood-2023-177679Van De VyverAJMarrer-BergerEWangKLehrTWalzAC. Cytokine release syndrome by T-cell-redirecting therapies: can we predict and modulate patient risk? Clin Cancer Res. (2021) 27:6083–94. doi: 10.1158/1078-0432.CCR-21-0470LiJPiskolRYbarraRChenYJLiJSlagaD. CD3 bispecific antibody-induced cytokine release is dispensable for cytotoxic T cell activity. Sci Transl Med. (2019) 11:eaax8861. doi: 10.1126/scitranslmed.aax8861JiangXChenXJaiprasartPCarpenterTJZhouRWangW. Development of a minimal physiologically-based pharmacokinetic/pharmacodynamic model to characterize target cell depletion and cytokine release for T cell-redirecting bispecific agents in humans. Eur J Pharm Sci. (2020) 146:105260. doi: 10.1016/j.ejps.2020.105260HosseiniIGadkarKStefanichELiCCSunLLChuYW. Mitigating the risk of cytokine release syndrome in a Phase I trial of CD20/CD3 bispecific antibody mosunetuzumab in NHL: impact of translational system modeling. NPJ Syst Biol Appl. (2020) 6:28. doi: 10.1038/s41540-020-00145-7ChenXKamperschroerCWongGXuanD. A modeling framework to characterize cytokine release upon T-cell-engaging bispecific antibody treatment: methodology and opportunities. Clin Transl Sci. (2019) 12:600–8. doi: 10.1111/cts.12662Movin-OsswaldGHammarlund-UdenaesM. Prolactin release after remoxipride by an integrated pharmacokinetic-pharmacodynamic model with intra and interindividual aspects. J Pharmacol Exp Ther. (1995) 274:921–7.RohatgiA. Webplotdigitizer: Version 4.6 (2022). Available online at: https://automeris.io/WebPlotDigitizer (Accessed July 19, 2024).LeclercqGServeraLADanilinSChallierJSteinhoffNBossenC. Dissecting the mechanism of cytokine release induced by T-cell engagers highlights the contribution of neutrophils. Oncoimmunology. (2022) 11:2039432. doi: 10.1080/2162402X.2022.2039432StaflinKZuch de ZafraCLSchuttLKClarkVZhongFHristopoulosM. Target arm affinities determine preclinical efficacy and safety of anti-HER2/CD3 bispecific antibody. JCI Insight. (2020) 5(7):e133757. doi: 10.1172/jci.insight.133757Parra-GuillenZPSchmidUJandaAFreiwaldMTroconizIF. Model-informed dose selection for xentuzumab, a dual insulin-like growth factor-I/II-neutralizing antibody. Clin Pharmacol Ther. (2020) 107:597–606. doi: 10.1002/cpt.1648Team RC. R: A language and environment for statistical computing. Vienna, Austria: R Foundation for Statistical Computing. (2022)FidlerMHallowMWilkinsJWangW. {RxODE}: Facilities for Simulating from ODE-Based Models. (2024). Available at: https://CRAN.R-project.org/package=rxode2.YuSZhangJYanYYaoXFangLXiongH. A novel asymmetrical anti-HER2/CD3 bispecific antibody exhibits potent cytotoxicity for HER2-positive tumor cells. J Exp Clin Cancer Res. (2019) 38:355. doi: 10.1186/s13046-019-1354-1BachmannMFOxeniusA. Interleukin 2: from immunostimulation to immunoregulation and back again. EMBO Rep. (2007) 8:1142–8. doi: 10.1038/sj.embor.7401099Negative feedback expansion of tregs caused by endogenous IL-2 limits the activity of IL-2-based therapies. J Cancer Immunol. (2023) 5(1):29–39. doi: 10.33696/cancerimmunol.5.074. R AILNGMB-MMWTB.FrancesNBacacMBray-FrenchKChristenFHintonHHusarE. Novel in vivo and in vitro pharmacokinetic/pharmacodynamic-based human starting dose selection for glofitamab. J Pharm Sci. (2022) 111:1208–18. doi: 10.1016/j.xphs.2021.12.019YiJSCoxMAZajacAJ. T-cell exhaustion: characteristics, causes and conversion. Immunology. (2010) 129:474–81. doi: 10.1111/j.1365-2567.2010.03255.xSimeoniMMagniPCammiaCDe NicolaoGCrociVPesentiE. Predictive pharmacokinetic-pharmacodynamic modeling of tumor growth kinetics in xenograft models after administration of anticancer agents. Cancer Res. (2004) 64:1094–101. doi: 10.1158/0008-5472.can-03-2524KuemmelCYangYZhangXFlorianJZhuHTegengeM. Consideration of a credibility assessment framework in model-informed drug development: potential application to physiologically-based pharmacokinetic modeling and simulation. CPT Pharmacometrics Syst Pharmacol. (2020) 9:21–8. doi: 10.1002/psp4.12479Tighiouart MouradRA. Dose finding with escalation with overdose control (EWOC) in cancer clinical trials. Stat Sci. (2010) 25(2):217–226. doi: 10.1214/10-STS333HarayaKTsutsuiHKomoriYTachibanaT. Recent advances in translational pharmacokinetics and pharmacodynamics prediction of therapeutic antibodies using modeling and simulation. Pharm (Basel). (2022) 15(5):508. doi: 10.3390/ph15050508DolinaJSVan Braeckel-BudimirNThomasGDSalek-ArdakaniS. CD8(+) T cell exhaustion in cancer. Front Immunol. (2021) 12:715234. doi: 10.3389/fimmu.2021.715234GaballaSNairRJacobsRWDevataSChoSGStevensDA. Double step-up dosing (2SUD) regimen mitigates severe icans and CRS
while maintaining high efficacy in subjects with relapsed/refractory (R/R) B-cell non-hodgkin
lymphoma (NHL) treated with AZD0486, a novel CD19xCD3 T-cell engager (TCE): updated safety and
efficacy data from the ongoing first-in-human ( FIH) phase 1 trial.
Blood. (2023) 142(Supplement
1):1662. doi: 10.1182/blood-2023-174668LytheGCallardREHoareRLMolina-ParisC. How many TCR clonotypes does a body maintain? J Theor Biol. (2016) 389:214–24. doi: 10.1016/j.jtbi.2015.10.016FribergLEHenningssonAMaasHNguyenLKarlssonMO. Model of chemotherapy-induced myelosuppression with parameter consistency across drugs. J Clin Oncol. (2002) 20:4713–21. doi: 10.1200/JCO.2002.02.140MacallanDCBuschRAsquithB. Current estimates of T cell kinetics in humans. Curr Opin Syst Biol. (2019) 18:77–86. doi: 10.1016/j.coisb.2019.10.002