Translational PK/PD and the first-in-human dose selection of a PD1/IL15: an engineered recombinant targeted cytokine for cancer immunotherapy.

Introduction: Interleukin 15 (IL-15) is a potential anticancer agent and numerous engineered IL-15 agonists are currently under clinical investigation. Selective targeting of IL-15 to specific lymphocytes may enhance therapeutic effects while helping to minimize toxicities. Methods: We designed and built a heterodimeric targeted cytokine (TaCk) that consists of an anti-programmed cell death 1 receptor antibody (anti-PD-1) and an engineered IL-15. This "PD1/IL15" selectively delivers IL-15 signaling to lymphocytes expressing PD-1. We then investigated the pharmacokinetic (PK) and pharmacodynamic (PD) effects of PD1/IL15 TaCk on immune cell subsets in cynomolgus monkeys after single and repeat intravenous dose administrations. We used these results to determine the first-in-human (FIH) dose and dosing frequency for early clinical trials. Results: The PD1/IL15 TaCk exhibited a nonlinear multiphasic PK profile, while the untargeted isotype control TaCk, containing an anti-respiratory syncytial virus antibody (RSV/IL15), showed linear and dose proportional PK. The PD1/IL15 TaCk also displayed a considerably prolonged PK (half-life range ∼1.0-4.1 days) compared to wild-type IL-15 (half-life ∼1.1 h), which led to an enhanced cell expansion PD response. The PD was dose-dependent, durable, and selective for PD-1+ lymphocytes. Notably, the dose- and time-dependent PK was attributed to dynamic TMDD resulting from test article-induced lymphocyte expansion upon repeat administration. The recommended first-in-human (FIH) dose of PD1/IL15 TaCk is 0.003 mg/kg, determined based on a minimum anticipated biological effect level (MABEL) approach utilizing a combination of in vitro and preclinical in vivo data. Conclusion: This work provides insight into the complex PK/PD relationship of PD1/IL15 TaCk in monkeys and informs the recommended starting dose and dosing frequency selection to support clinical evaluation of this novel targeted cytokine.

A Minimal PBPK/PD Model with Expansion-Enhanced Target-Mediated Drug Disposition to Support a First-in-Human Clinical Study Design for a FLT3L-Fc Molecule.

FLT3L-Fc is a half-life extended, effectorless Fc-fusion of the native human FLT3-ligand. In cynomolgus monkeys, treatment with FLT3L-Fc leads to a complex pharmacokinetic/pharmacodynamic (PK/PD) relationship, with observed nonlinear PK and expansion of different immune cell types across different dose levels. A minimal physiologically based PK/PD model with expansion-enhanced target-mediated drug disposition (TMDD) was developed to integrate the molecule's mechanism of action, as well as the complex preclinical and clinical PK/PD data, to support the preclinical-to-clinical translation of FLT3L-Fc. In addition to the preclinical PK data of FLT3L-Fc in cynomolgus monkeys, clinical PK and PD data from other FLT3-agonist molecules (GS-3583 and CDX-301) were used to inform the model and project the expansion profiles of conventional DC1s (cDC1s) and total DCs in peripheral blood. This work constitutes an essential part of our model-informed drug development (MIDD) strategy for clinical development of FLT3L-Fc by projecting PK/PD in healthy volunteers, determining the first-in-human (FIH) dose, and informing the efficacious dose in clinical settings. Model-generated results were incorporated in regulatory filings to support the rationale for the FIH dose selection.

Preclinical and translational pharmacology of afucosylated anti-CCR8 antibody for depletion of tumour-infiltrating regulatory T cells.

BACKGROUND AND PURPOSE: RO7502175 is an afucosylated antibody designed to eliminate C-C motif chemokine receptor 8 (CCR8)+ Treg cells in the tumour microenvironment through enhanced antibody-dependent cellular cytotoxicity (ADCC). EXPERIMENTAL APPROACH: We report findings from preclinical studies characterizing pharmacology, pharmacokinetics (PK)/pharmacodynamics (PD) and safety profile of RO7502175 and discuss the translational PK/PD approach used to inform first-in-human (FiH) dosing strategy and clinical development in solid tumour indications. KEY RESULTS: RO7502175 demonstrated selective ADCC against human CCR8+ Treg cells from dissociated tumours in vitro. In cynomolgus monkeys, RO7502175 exhibited a biphasic concentration-time profile consistent with immunoglobulin G1 (IgG1) antibodies, reduced CCR8+ Treg cells in the blood, induced minimal and transient cytokine secretion, and was well tolerated with a no-observed-adverse-effect level (NOAEL) of 100 mg·kg-1. Moreover, RO7502175 caused minimal cytokine release from peripheral blood mononuclear cells (PBMCs) in vitro. A quantitative model was developed to capture surrogate anti-murine CCR8 antibody PK/PD and tumour dynamics in mice and RO7502175 PK/PD in cynomolgus monkeys. Subsequently, the model was used to project RO7502175 human PK and receptor occupancy (RO) in patients. Because traditional approaches resulted in a low FiH dose for this molecule, even with its superior preclinical safety profile, an integrated approach based on the totality of preclinical data and modelling insights was used for starting dose selection. CONCLUSION AND IMPLICATIONS: This work demonstrates a translational research strategy for collecting and utilizing relevant nonclinical data, developing a mechanistic PK/PD model and using a comprehensive approach to inform clinical study design for RO7502175.

Nonclinical pharmacokinetics, pharmacodynamics and safety assessment of a FLT3L-Fc molecule for cancer immunotherapy.

FLT3L-Fc is a cytokine-Fc fusion agonizing receptor-type tyrosine-protein kinase FLT3 (fms-related tyrosine kinase 3; CD135). FLT3 is expressed on dendritic cells (DCs) as well as myeloid and lymphoid progenitors. Nonclinical pharmacokinetics, pharmacodynamics and safety of FLT3L-Fc were investigated in rats and cynomolgus monkeys. FLT3L-Fc induced robust pharmacodynamic responses, evidenced by marked expansion of peripheral blood cDC1s, cDC2s, and pDCs (up to 301-fold in rats and 378-fold in monkeys), peaking at 8-10 days after the first dose. FLT3L-Fc was well tolerated with no adverse findings at doses up to 10 mg/kg administered intravenously twice three weeks apart. In both species, major clinical pathology findings consisted of expansion of white blood cell (WBC) populations including lymphocytes, monocytes, neutrophils, basophils, and large unstained cells, which were pronounced after the first dose. The WBC findings were associated microscopically with histiocytic and mononuclear cell infiltrates in multiple organs. Tissue immunohistochemistry in monkeys showed that the leukocyte infiltrates consisted of hematopoietic progenitor cells and histiocytes with a reactive morphology and were associated with a slight stimulation of regional T and B cell populations. Additional FLT3L-Fc-associated changes included decreases in red blood cell (RBC) mass, increases in RBC distribution width, variable changes in reticulocytes, and transient alterations in platelet counts (rats only). The RBC and WBC findings were associated microscopically with increased hematopoietic cellularity of the bone marrow in both species and increased splenic megakaryocytic extramedullary hematopoiesis in rats. The totality of nonclinical safety data support the clinical development of FLT3L-Fc.

"Digital twins elucidate critical role of Tscm in clinical persistence of TCR-engineered cell therapy".

Despite recent progress in adoptive T cell therapy for cancer, understanding and predicting the kinetics of infused T cells remains a challenge. Multiple factors can impact the distribution, expansion, and decay or persistence of infused T cells in patients. We have developed a novel quantitative systems pharmacology (QSP) model of TCR-transgenic T cell therapy in patients with solid tumors to describe the kinetics of endogenous T cells and multiple memory subsets of engineered T cells after infusion. These T cells undergo lymphodepletion, proliferation, trafficking, differentiation, and apoptosis in blood, lymph nodes, tumor site, and other peripheral tissues. Using the model, we generated patient-matched digital twins that recapitulate the circulating T cell kinetics reported from a clinical trial of TCR-engineered T cells targeting E7 in patients with metastatic HPV-associated epithelial cancers. Analyses of key parameters influencing cell kinetics and differences among digital twins identify stem cell-like memory T cells (Tscm) cells as an important determinant of both expansion and persistence and suggest that Tscm-related differences contribute significantly to the observed variability in cellular kinetics among patients. We simulated in silico clinical trials using digital twins and predict that Tscm enrichment in the infused product improves persistence of the engineered T cells and could enable administration of a lower dose. Finally, we verified the broader relevance of the QSP model, the digital twins, and findings on the importance of Tscm enrichment by predicting kinetics for two patients with pancreatic cancer treated with KRAS G12D targeting T cell therapy. This work offers insight into the key role of Tscm biology on T cell kinetics and provides a quantitative framework to evaluate cellular kinetics for future efforts in the development and clinical application of TCR-engineered T cell therapies.

Complex PK-PD of an engineered IL-15/IL-15Rα-Fc fusion protein in cynomolgus monkeys: QSP modeling of lymphocyte dynamics.

XmAb24306 is a lymphoproliferative interleukin (IL)-15/IL-15 receptor α (IL-15Rα) Fc-fusion protein currently under clinical investigation as an immunotherapeutic agent for cancer treatment. XmAb24306 contains mutations in IL-15 that attenuate its affinity to the heterodimeric IL-15 receptor ßγ (IL-15R). We observe substantially prolonged pharmacokinetics (PK) (half-life ∼ 2.5 to 4.5 days) in single- and repeat-dose cynomolgus monkey (cyno) studies compared to wild-type IL-15 (half-life ∼ 1 hour), leading to increased exposure and enhanced and durable expansion of NK cells, CD8+ T cells and CD4-CD8- (double negative [DN]) T cells. Drug clearance varied with dose level and time post-dose, and PK exposure decreased upon repeated dosing, which we attribute to increased target-mediated drug disposition (TMDD) resulting from drug-induced lymphocyte expansion (i.e., pharmacodynamic (PD)-enhanced TMDD). We developed a quantitative systems pharmacology (QSP) model to quantify the complex PKPD behaviors due to the interactions of XmAb24306 with multiple cell types (CD8+, CD4+, DN T cells, and NK cells) in the peripheral blood (PB) and lymphoid tissues. The model, which includes nonspecific drug clearance, binding to and TMDD by IL15R differentially expressed on lymphocyte subsets, and resultant lymphocyte margination/migration out of PB, expansion in lymphoid tissues, and redistribution to the blood, successfully describes the systemic PK and lymphocyte kinetics observed in the cyno studies. Results suggest that after 3 doses of every-two-week (Q2W) doses up to 70 days, the relative contributions of each elimination pathway to XmAb24306 clearance are: DN T cells > NK cells > CD8+ T cells > nonspecific clearance > CD4+ T cells. Modeling suggests that observed cellular expansion in blood results from the influx of cells expanded by the drug in lymphoid tissues. The model is used to predict lymphoid tissue expansion and to simulate PK-PD for different dose regimens. Thus, the model provides insight into the mechanisms underlying the observed PK-PD behavior of an engineered cytokine and can serve as a framework for the rapid integration and analysis of data that emerges from ongoing clinical studies in cancer patients as single-agent or given in combination.

Systems-based digital twins to help characterize clinical dose-response and propose predictive biomarkers in a Phase I study of bispecific antibody, mosunetuzumab, in NHL.

Phase I oncology clinical trials often comprise a limited number of patients representing different disease subtypes who are divided into cohorts receiving treatment(s) at different dosing levels and schedules. Here, we leverage a previously developed quantitative systems pharmacology model of the anti-CD20/CD3 T-cell engaging bispecific antibody, mosunetuzumab, to account for different dosing regimens and patient heterogeneity in the phase I study to inform clinical dose/exposure-response relationships and to identify biological determinants of clinical response. We developed a novel workflow to generate digital twins for each patient, which together form a virtual population (VPOP) that represented variability in biological, pharmacological, and tumor-related parameters from the phase I trial. Simulations based on the VPOP predict that an increase in mosunetuzumab exposure increases the proportion of digital twins with at least a 50% reduction in tumor size by day 42. Simulations also predict a left-shift of the exposure-response in patients diagnosed with indolent compared to aggressive non-Hodgkin's lymphoma (NHL) subtype; this increased sensitivity in indolent NHL was attributed to the lower inferred values of tumor proliferation rate and baseline T-cell infiltration in the corresponding digital twins. Notably, the inferred digital twin parameters from clinical responders and nonresponders show that the potential biological difference that can influence response include tumor parameters (tumor size, proliferation rate, and baseline T-cell infiltration) and parameters defining the effect of mosunetuzumab on T-cell activation and B-cell killing. Finally, the model simulations suggest intratumor expansion of pre-existing T-cells, rather than an influx of systemically expanded T-cells, underlies the antitumor activity of mosunetuzumab.

Current practices for QSP model assessment: an IQ consortium survey.

Quantitative Systems Pharmacology (QSP) modeling is increasingly applied in the pharmaceutical industry to influence decision making across a wide range of stages from early discovery to clinical development to post-marketing activities. Development of standards for how these models are constructed, assessed, and communicated is of active interest to the modeling community and regulators but is complicated by the wide variability in the structures and intended uses of the underlying models and the diverse expertise of QSP modelers. With this in mind, the IQ Consortium conducted a survey across the pharmaceutical/biotech industry to understand current practices for QSP modeling. This article presents the survey results and provides insights into current practices and methods used by QSP practitioners based on model type and the intended use at various stages of drug development. The survey also highlights key areas for future development including better integration with statistical methods, standardization of approaches towards virtual populations, and increased use of QSP models for late-stage clinical development and regulatory submissions.

Quantitative Clinical Pharmacology of T-Cell Engaging Bispecifics: Current Perspectives and Opportunities.

T-cell directing/engaging bispecifics (TDBs) enable a powerful mode of action by activating T-cells through the creation of artificial immune synapses. Their pharmacological response involves the dynamic inter-relationships among T-cells, tumor cells, and TDBs. This results in complex and challenging issues in understanding pharmacokinetics, tissue distribution, target engagement, and exposure-response relationship. Dosing strategy plays a crucial role in determining the therapeutic window of TDBs because of the desire to maximize therapeutic efficacy in the context of known mechanism-related adverse events, such as cytokine release syndrome and neurological adverse events. Such adverse events are commonly reported as the most prominent events during the initial treatment cycles and dissipate over time. Therefore, the kinetic characterization of the inter-relationships between exposure/target engagement and safety/efficacy outcomes is crucial in designing the optimal dosing regimen to maximize the benefit/risk of TDB agents. In this review, we discuss the key clinical pharmacological considerations in drug discovery and development for TDBs and provide a summary of TDBs currently in clinical development. We also propose forward-looking perspectives and opportunities to derive insights through quantitative clinical pharmacology approaches.

Personalized Cancer Vaccines: Clinical Landscape, Challenges, and Opportunities.

Tremendous innovation is underway among a rapidly expanding repertoire of promising personalized immune-based treatments. Therapeutic cancer vaccines (TCVs) are attractive systemic immunotherapies that activate and expand antigen-specific CD8+ and CD4+ T cells to enhance anti-tumor immunity. Our review highlights key issues impacting TCVs in clinical practice and reports on progress in development. We review the mechanism of action, immune-monitoring, dosing strategies, combinations, obstacles, and regulation of cancer vaccines. Most trials of personalized TCVs are ongoing and represent diverse platforms with predominantly early investigations of mRNA, DNA, or peptide-based targeting strategies against neoantigens in solid tumors, with many in combination immunotherapies. Multiple delivery systems, routes of administration, and dosing strategies are used. Intravenous or intramuscular administration is common, including delivery by lipid nanoparticles. Absorption and biodistribution impact antigen uptake, expression, and presentation, affecting the strength, speed, and duration of immune response. The emerging trials illustrate the complexity of developing this class of innovative immunotherapies. Methodical testing of the multiple potential factors influencing immune responses, as well as refined quantitative methodologies to facilitate optimal dosing strategies, could help resolve uncertainty of therapeutic approaches. To increase the likelihood of success in bringing these medicines to patients, several unique development challenges must be overcome.

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.

Mosunetuzumab, a T-cell dependent bispecific antibody that binds CD3 and CD20 to drive T-cell mediated B-cell killing, is currently being tested in non-Hodgkin lymphoma. However, potent immune stimulation with T-cell directed therapies poses the risk of cytokine release syndrome, potentially limiting dose and utility. To understand mechanisms behind safety and efficacy and explore safety mitigation strategies, we developed a novel mechanistic model of immune and antitumor responses to the T-cell bispecifics (mosunetuzumab and blinatumomab), including the dynamics of B- and T-lymphocytes in circulation, lymphoid tissues, and tumor. The model was developed and validated using mosunetuzumab nonclinical and blinatumomab clinical data. Simulations delineated mechanisms contributing to observed cell and cytokine (IL6) dynamics and predicted that initial step-fractionated dosing limits systemic T-cell activation and cytokine release without compromising tumor response. These results supported a change to a step-fractionated treatment schedule of mosunetuzumab in the ongoing Phase I clinical trial, enabling safer administration of higher doses.

gQSPSim: A SimBiology-Based GUI for Standardized QSP Model Development and Application.

Quantitative systems pharmacology (QSP) models are often implemented using a wide variety of technical workflows and methodologies. To facilitate reproducibility, transparency, portability, and reuse for QSP models, we have developed gQSPSim, a graphical user interface-based MATLAB application that performs key steps in QSP model development and analyses. The capabilities of gQSPSim include (i) model calibration using global and local optimization methods, (ii) development of virtual subjects to explore variability and uncertainty in the represented biology, and (iii) simulations of virtual populations for different interventions. gQSPSim works with SimBiology-built models using components such as species, doses, variants, and rules. All functionalities are equipped with an interactive visualization interface and the ability to generate presentation-ready figures. In addition, standardized gQSPSim sessions can be shared and saved for future extension and reuse. In this work, we demonstrate gQSPSim's capabilities with a standard target-mediated drug disposition model and a published model of anti-proprotein convertase subtilisin/kexin type 9 (PCSK9) treatment of hypercholesterolemia.

gPKPDSim: a SimBiology®-based GUI application for PKPD modeling in drug development.

Modeling and simulation (M&S) is increasingly used in drug development to characterize pharmacokinetic-pharmacodynamic (PKPD) relationships and support various efforts such as target feasibility assessment, molecule selection, human PK projection, and preclinical and clinical dose and schedule determination. While model development typically require mathematical modeling expertise, model exploration and simulations could in many cases be performed by scientists in various disciplines to support the design, analysis and interpretation of experimental studies. To this end, we have developed a versatile graphical user interface (GUI) application to enable easy use of any model constructed in SimBiology® to execute various common PKPD analyses. The MATLAB®-based GUI application, called gPKPDSim, has a single screen interface and provides functionalities including simulation, data fitting (parameter estimation), population simulation (exploring the impact of parameter variability on the outputs of interest), and non-compartmental PK analysis. Further, gPKPDSim is a user-friendly tool with capabilities including interactive visualization, exporting of results and generation of presentation-ready figures. gPKPDSim was designed primarily for use in preclinical and translational drug development, although broader applications exist. gPKPDSim is a MATLAB®-based open-source application and is publicly available to download from MATLAB® Central™. We illustrate the use and features of gPKPDSim using multiple PKPD models to demonstrate the wide applications of this tool in pharmaceutical sciences. Overall, gPKPDSim provides an integrated, multi-purpose user-friendly GUI application to enable efficient use of PKPD models by scientists from various disciplines, regardless of their modeling expertise.

Multiplexed Component Analysis to Identify Genes Contributing to the Immune Response during Acute SIV Infection.

Immune response genes play an important role during acute HIV and SIV infection. Using an SIV macaque model of AIDS and CNS disease, our overall goal was to assess how the expression of genes associated with immune and inflammatory responses are longitudinally changed in different organs or cells during SIV infection. To compare RNA expression of a panel of 88 immune-related genes across time points and among three tissues - spleen, mesenteric lymph nodes (MLN) and peripheral blood mononuclear cells (PBMC) - we designed a set of Nanostring probes. To identify significant genes during acute SIV infection and to investigate whether these genes are tissue-specific or have global roles, we introduce a novel multiplexed component analysis (MCA) method. This combines multivariate analysis methods with multiple preprocessing methods to create a set of 12 "judges"; each judge emphasizes particular types of change in gene expression to which cells could respond, for example, the absolute or relative size of expression change from baseline. Compared to bivariate analysis methods, our MCA method improved classification rates. This analysis allows us to identify three categories of genes: (a) consensus genes likely to contribute highly to the immune response; (b) genes that would contribute highly to the immune response only if certain assumptions are met - e.g. that the cell responds to relative expression change rather than absolute expression change; and (c) genes whose contribution to immune response appears to be modest. We then compared the results across the three tissues of interest; some genes are consistently highly-contributing in all tissues, while others are specific for certain tissues. Our analysis identified CCL8, CXCL10, CXCL11, MxA, OAS2, and OAS1 as top contributing genes, all of which are stimulated by type I interferon. This suggests that the cytokine storm during acute SIV infection is a systemic innate immune response against viral replication. Furthermore, these genes have approximately equal contributions to all tissues, making them possible candidates to be used as non-invasive biomarkers in studying PBMCs instead of MLN and spleen during acute SIV infection experiments. We identified clusters of genes that co-vary together and studied their correlation with regard to other gene clusters. We also developed novel methods to faithfully visualize multi-gene correlations on two-dimensional polar plots, and to visualize tissue specificity of gene expression responses.

APOBEC3G-Augmented Stem Cell Therapy to Modulate HIV Replication: A Computational Study.

The interplay between the innate immune system restriction factor APOBEC3G and the HIV protein Vif is a key host-retrovirus interaction. APOBEC3G can counteract HIV infection in at least two ways: by inducing lethal mutations on the viral cDNA; and by blocking steps in reverse transcription and viral integration into the host genome. HIV-Vif blocks these antiviral functions of APOBEC3G by impeding its encapsulation. Nonetheless, it has been shown that overexpression of APOBEC3G, or interfering with APOBEC3G-Vif binding, can efficiently block in vitro HIV replication. Some clinical studies have also suggested that high levels of APOBEC3G expression in HIV patients are correlated with increased CD4+ T cell count and low levels of viral load; however, other studies have reported contradictory results and challenged this observation. Stem cell therapy to replace a patient's immune cells with cells that are more HIV-resistant is a promising approach. Pre-implantation gene transfection of these stem cells can augment the HIV-resistance of progeny CD4+ T cells. As a protein, APOBEC3G has the advantage that it can be genetically encoded, while small molecules cannot. We have developed a mathematical model to quantitatively study the effects on in vivo HIV replication of therapeutic delivery of CD34+ stem cells transfected to overexpress APOBEC3G. Our model suggests that stem cell therapy resulting in a high fraction of APOBEC3G-overexpressing CD4+ T cells can effectively inhibit in vivo HIV replication. We extended our model to simulate the combination of APOBEC3G therapy with other biological activities, to estimate the likelihood of improved outcomes.

Multi-scale modeling of HIV infection in vitro and APOBEC3G-based anti-retroviral therapy.

The human APOBEC3G is an innate restriction factor that, in the absence of Vif, restricts HIV-1 replication by inducing excessive deamination of cytidine residues in nascent reverse transcripts and inhibiting reverse transcription and integration. To shed light on impact of A3G-Vif interactions on HIV replication, we developed a multi-scale computational system consisting of intracellular (single-cell), cellular and extracellular (multicellular) events by using ordinary differential equations. The single-cell model describes molecular-level events within individual cells (such as production and degradation of host and viral proteins, and assembly and release of new virions), whereas the multicellular model describes the viral dynamics and multiple cycles of infection within a population of cells. We estimated the model parameters either directly from previously published experimental data or by running simulations to find the optimum values. We validated our integrated model by reproducing the results of in vitro T cell culture experiments. Crucially, both downstream effects of A3G (hypermutation and reduction of viral burst size) were necessary to replicate the experimental results in silico. We also used the model to study anti-HIV capability of several possible therapeutic strategies including: an antibody to Vif; upregulation of A3G; and mutated forms of A3G. According to our simulations, A3G with a mutated Vif binding site is predicted to be significantly more effective than other molecules at the same dose. Ultimately, we performed sensitivity analysis to identify important model parameters. The results showed that the timing of particle formation and virus release had the highest impacts on HIV replication. The model also predicted that the degradation of A3G by Vif is not a crucial step in HIV pathogenesis.