Next Article in Journal
Sweeteners Show a Plasticizing Effect on PVP K30—A Solution for the Hot-Melt Extrusion of Fixed-Dose Amorphous Curcumin-Hesperetin Solid Dispersions
Next Article in Special Issue
Analysis of Lipophilicity and Pharmacokinetic Parameters of Dipyridothiazine Dimers with Anticancer Potency
Previous Article in Journal
Critical View on the Qualification of Electronic Tongues Regarding Their Performance in the Development of Peroral Drug Formulations with Bitter Ingredients
Previous Article in Special Issue
Enhancing Paclitaxel Efficacy with Piperine-Paclitaxel Albumin Nanoparticles in Multidrug-Resistant Triple-Negative Breast Cancer by Inhibiting P-Glycoprotein
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

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

Genentech, Inc., South San Francisco, CA 94080, USA
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Pharmaceutics 2024, 16(5), 660; https://doi.org/10.3390/pharmaceutics16050660
Submission received: 28 March 2024 / Revised: 27 April 2024 / Accepted: 3 May 2024 / Published: 15 May 2024
(This article belongs to the Special Issue Role of Pharmacokinetics in Drug Development and Evaluation)

Abstract

:
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.

1. Introduction

Cancer immunotherapy (CIT) boosts the immune system’s ability to identify and eliminate cancer cells. Checkpoint inhibitors including anti-programmed death-1 (PD-1) or anti-programmed death-ligand 1 (PD-L1) therapies such as atezolizumab, nivolumab, and pembrolizumab, have been approved for several oncology indications. However, the majority of patients do not respond to treatment with checkpoint inhibitors, or they only experience transient disease stabilization. Thus, a significant unmet need persists for patients with metastatic solid tumors [1,2,3] necessitating the application of innovative and informed strategies such as model-informed drug development (MIDD) for effective dose-finding and the efficacy prediction of novel molecules. It has been hypothesized that low dendritic cell (DC) counts within the tumor microenvironment could hinder the effectiveness of CIT therapy in cancer patients. Among DC subsets, conventional DC1s (cDC1s) reportedly have high intrinsic capacity to cross-present internalized antigens via MHC class I, enabling CD8+ T cell activation and inducing anti-tumor T cell-mediated immunity [4,5,6,7]. cDC1s drive the initiation of the cancer immunity cycle by priming T cells and maintaining antigen-specific T cell recall responses [8,9,10]. cDC1s are scarce within tumors, due in part to the immunosuppressive nature of the tumor microenvironment which leads to suboptimal recruitment, differentiation, or viability of cDC1s [11,12]. Conversely, an abundance of cDC1s in the tumor microenvironment has been associated with increased T cell infiltration, the response to checkpoint inhibitors, and overall survival rates in patients with cancer [13,14,15].
Fms-related tyrosine kinase 3 ligand (FLT3L) is a soluble growth factor (or cytokine) that plays a critical role in regulating hematopoiesis. FLT3L functions by agonizing the fms-related tyrosine kinase 3 receptor (FLT3; also known as CD135), which is expressed on DCs as well as myeloid and lymphoid progenitor cells [16,17,18,19,20]. Agonism of this pathway initiates signaling cascades that promote DC proliferation, differentiation, mobilization, and survival [21,22,23,24,25,26,27]. FLT3L treatment has been shown to expand not only cDC1s, but also conventional type 2 DCs (cDC2s), plasmacytoid DCs (pDCs), and other immune cells of myeloid origin in blood, lymphoid organs, and tumors [28,29,30,31]. Preclinical proof-of-concept studies have demonstrated that FLT3L treatment improves anti-tumor efficacy by enhancing T cell priming [32,33,34,35,36,37,38], providing a rationale for combining FLT3L therapy with T cell-directed CIT including checkpoint inhibition.
Clinical studies have demonstrated that targeting the FLT3 signaling cascade is a promising approach to safely expanding DCs and enhancing anti-cancer immunity in patients. CDX-301, a recombinant human FLT3L molecule (Celldex Therapeutics, Inc., Hampton, NJ, USA), was tested as a single agent in healthy volunteers, as well as in combination with immune checkpoint inhibitors, radiation therapy, and other innate immune-activating agents in phase I/II clinical trials [23,34,36,39]. In healthy volunteers, CDX-301 was well-tolerated at 25 μg/kg/day for 7 or 10 days and 75 μg/kg/day for 5 days [23]. In another study, recombinant human FLT3L was well-tolerated when dosed up to 100 μg/kg/day for 14 days [22]. GS-3583, another FLT3L agonist with a crystallizable fragment (Fc) fusion protein molecule (FLT3L-Fc) (Gilead Sciences, Inc., Foster City, CA, USA), was recently tested in healthy volunteers as a single dose and found to be well-tolerated up to 2000 μg (40). Both recombinant and Fc fusion FLT3L formats demonstrated significant expansion of DCs in healthy volunteers [22,23,40] and, in combination with adjuvants, can potentially promote T cell-mediated anti-tumor activity in cancer patients. The PK properties of recombinant human FLT3L molecules, e.g., CDX-301, necessitate daily dosing, presenting logistical considerations that potentially impede the broader application of FLT3L-targeting treatments in clinical settings.
FLT3L-Fc (Genentech, Inc., South San Francisco, CA, USA) is a fusion of the truncated extracellular domain of human FLT3L and the Fc portion of human IgG1. The Fc portion has been further engineered to impair binding to human Fcγ receptors, thereby attenuating Fc-mediated effector functions [41]. The FLT3L-Fc molecule has been designed to extend the FLT3L half-life and decrease the clinical dosing frequency required to achieve therapeutically significant cDC1 expansion. It thus has the potential to enhance anti-tumor immunity, especially when used in combination with DC maturation agents and checkpoint inhibitor immunotherapy.
A minimal physiologically based pharmacokinetic/pharmacodynamic (PBPK/PD) model [42,43] with expansion-enhanced target-mediated drug disposition (TMDD) was developed to (1) integrate the mechanism of action of FLT3L-Fc, (2) explain and provide insight into the complex PK/PD relationship of FLT3L-Fc and CDX-301 observed in preclinical and clinical settings, and (3) provide a rationale for the selection of the first-in-human (FIH) dose in healthy volunteers. The model was first used to capture PK of FLT3L-Fc in cynomolgus monkeys and then used to project clinical PK by translating model parameter values, including clearance, physiological volumes, and flow rates [42,43,44], and accounting for clinical target-related parameters. Given that FLT3L-Fc and CDX-301 are expected to exhibit comparable PK/PD relationships due to their similar affinities to FLT3 receptors and mechanisms of action [22,23], the clinical PK/PD data from healthy volunteers treated with CDX-301 were incorporated into the model and used to project cDC1 or total DC expansion in peripheral blood after administration of FLT3L-Fc to healthy volunteers. This work describes the model and the MIDD strategy used to project a FIH dose aimed at achieving a minimum pharmacologically active dose (mPAD) level with an adequate safety margin for a Phase IA study (clinical trial number: 92655801, www.isrctn.com, accessed on 26 April 2024), in order to support clinical study design and regulatory submissions.

2. Materials and Methods

2.1. Ethics Statement

The authors confirm that they have obtained appropriate institutional review board approval for all animal experimental investigations at Charles River Laboratories (CRL, Reno, NV, USA and Sherbrooke, QC, Canada). All procedures were approved by the Institutional Animal Care and Use Committee (IACUC) at CRL (protocol code 20239581, 30 March 2020; protocol code 5003546, 7 December 2020) and were performed in compliance with the Animal Welfare Act, the Guide for the Care and Use of Laboratory Animals, and the Office of Laboratory Animal Welfare.

2.2. Datasets Used in the Workflow of Minimal PBPK/PD Model Development

Preclinical and clinical data were obtained from the published literature or the studies conducted by Genentech. A summary of the data used for model calibration or validation, including study titles, dose levels/regimens, measurements, and the source of the data is provided in Table S1. The list of studies includes two internal preclinical studies in cynomolgus monkeys using FLT3L-Fc (Genentech Studies) and three clinical studies from other FLT3L-based therapies [22,23,40]. The PK of FLT3L-Fc in preclinical studies 1 and 2 (Genentech Studies) were measured using ELISA methods (anti-FLT3L capture and anti-human Fc detection) with a lower limit of quantification (LLOQ) of 31.3 ng/mL for preclinical study 1 and 7.8 ng/mL for preclinical study 2. Detection of anti-drug antibodies (ADAs) to FLT3L-Fc for preclinical study 1 utilized a direct ELISA where FLT3L-Fc was directly coated on the ELISA plate, and cynomolgus monkey IgG directed against FLT3L-Fc was detected with an anti-cynomolgus monkey IgG antibody. Preclinical study 2 utilized a bridging ELISA with biotin- and digoxigenin-conjugated RO7497987 to capture and detect ADAs directed against RO7497987. The PK/PD of CDX-301 and recombinant human FLT3L in clinical studies 1 and 2, and the PK of GS-3583 in clinical study 3 were digitized using “GetData Graph Digitizer” version 2.26 (http://getdata-graph-digitizer.com, accessed on 26 April 2024). Since other recombinant human FLT3Ls are expected to behave similarly to CDX-301, we will henceforth refer to them collectively as CDX-301 for the purpose of our modeling. Figure 1 (created with Biorender) provides a summary of the workflow and the use of data for model calibration and validation.

2.3. Model Structure and Key Assumptions

2.3.1. PK Model Structure of FLT3L-Fc and GS-3583

A minimal PBPK model with expansion-enhanced TMDD was developed to describe the available FLT3L-Fc PK data from cynomolgus monkeys and project FLT3L-Fc PK profiles (Figure 1A). The minimal PBPK model structure and modeling assumptions were previously described by Jusko et al. [42,45]. Briefly, the minimal PBPK model integrates physiological factors and divides tissue spaces into two groups based on their vascular endothelial structure, leaky and tight, which subsequently provides two general sites of distribution. The model also includes convection as the primary mechanism of antibody trafficking from plasma into tissues and return to plasma [42,45]. In this work, free FLT3L-Fc or GS-3583 in the central compartment binds to free FLT3 receptors, forming single-bound or double-bound complexes with the receptors. Subsequently, homodimer (double-bound) complexes induce signaling cascades that promote the expansion of the available target. The differential equations describing the minimal PBPK model with expansion-enhanced TMDD are shown below:
d C p d t = D o s e V p k o n 1 · C p · C R + k o f f 1 · C R p + ( C l y m p h · L C p · L t i g h t · 1 σ t i g h t C p · L l e a k y · 1 σ l e a k y C p · C L p ) / V p
d C R d t = k s y n k d e g · C R k o n 1 · C p · C R + k o f f 1 · C R p k o n 2 · C R · C R p + k o f f 2 · C R p R + C R T O T · v m p r o l i f · R O D B 100 α R O D B 100 α + k m p r o l i f α
d C R p d t = k o n 1 · C p · C R k o f f 1 · C R p k o n 2 · C R · C R p + k o f f 2 · C R p R k d e g · C R p
d C R p R d t = k o n 2 · C p · C R p k o f f 2 · C R p R k d e g · C R p R
where Cp is the drug concentration in the central compartment, and Ctight, Cleaky, and Clymph are drug concentrations in the tight, leaky, and lymph compartments, respectively. Vp denotes the species-specific volume of the central compartment, Vtight represents the interstitial fluid volumes of tissues with non-fenestrated capillaries (e.g., muscle, skin, adipose, and brain), and Vleaky denotes interstitial fluid volumes of tissues with fenestrated capillaries (e.g., heart, liver, kidney, etc.) [46]. Ltight and Lleaky are the flows to collective tight and leaky tissues, and L represents the total lymph flow equal to the sum of Ltight and Lleaky. Parameters σtight and σleaky are vascular reflection coefficients for tight and leaky tissues, and σlymph is the lymphatic reflection coefficient [42]. CLp denotes the non-specific clearance from the plasma. CR, CRp, and CRpR represent concentrations of the free, single-bound, and double-bound forms of the FLT3 target. CR is assumed to be at homeostasis in the absence of the drug and is synthesized and degraded at zero-order and first-order rate constants ksyn and kdeg, respectively. Cp binds to CR and dissociates from CRp via rate constants kon1 and koff1, respectively. CRp binds to CR and dissociates from CRpR via rate constants kon2 and koff2, respectively. Note that it is assumed that the degradation rate constant for the complexes is the same as the degradation rate constant for the free receptor. Complex formation between FLT3L-Fc and two FLT3 receptors leads to downstream intracellular signaling and subsequent expansion of the total pool of FLT3 receptors [23,47]. In the model, the percentage of FLT3 receptors bound in a double-bound complex (RODB or double-bound receptor occupancy) expands the pool of free FLT3 receptors. Parameter vmprolif represents the maximum rate of FLT3 receptor expansion, kmprolif denotes the percentage of FLT3 RODB to achieve 50% of maximum expansion, and α represents the Hill coefficient. RODB is calculated as follows:
R O D B = 2 · C R p R C R T O T × 100
C R T O T = C R + C R p + 2 · C R p R
where CRTOT represents the total pool of FLT3 receptors. ROSB and RO represent the percentage of FLT3 receptors bound in a single-bound complex and the percentage of the total FLT3 receptors bound, respectively, and are calculated as follows:
R O S B = C R p C R T O T × 100
R O = 1 C R C R T O T × 100

2.3.2. PK Model Structure of CDX-301

A standard single-compartment PK model with subcutaneous absorption and nonlinear elimination was used to describe the measured concentration of CDX-301 from healthy volunteers (Figure 1B). The CDX-301 PK model is shown below:
d A S C , C D X d t = D o s e C D X k a b s · A S C ,   C D X
d A C ,   C D X d t = k a b s · A S C , C D X V m ,   C D X · C C , C D X K m ,   C D X + C C , C D X
C C , C D X = A C , C D X V d , C D X
where ASC,CDX and AC,CDX represent the amount of CDX-301 in the subcutaneous and central compartments, respectively, and CC,CDX represents the concentration of CDX-301 in the central compartment. kabs denotes the absorption rate of CDX-301 from the subcutaneous compartment into the central compartment. Vd,CDX denotes the volume of the central compartment. Vm,CDX represents the maximum nonlinear elimination rate of CDX-301, and Km,CDX denotes the concentration of CDX-301 needed to achieve 50% of Vm,CDX.

2.3.3. PD Model Structure for DC Expansion

FLT3L-Fc and CDX-301 are assumed to elicit the same exposure response on cDC1 and cDC2 expansion [22,23], and, hence, concentrations of FLT3L-Fc and CDX-301 in nM can be used interchangeably to drive cell expansion (Figure 1C). The expansion of immune cells increases the total FLT3 receptors, thereby enhancing the TMDD. The expansion-enhanced TMDD was previously modeled using an approach that decouples PK and PD [48,49]. The current model adopts this decoupled approach wherein unbound drug concentration drives a delayed expansion of cDCs following the initiation of dosing via a series of transit compartments. cDCx (x takes values of 1 and 2 to represent cDC1 or cDC2, respectively) proliferates in a drug-dependent and target-dependent manner [23], as represented by the following equation:
d c D C x d t = k d e g x · c D C x , 0 k d e g x · c D C x + c D C x · k d e g x · v m x · C 3 x n x k m x n x + C 3 x n x f D C 1 · k d e g 2 · m a x c D C x c D C x , 0 , 0 2 c D C x , 0
where cDCx,0 represents cDCx counts at baseline and C3x represents the drug-concentration in the third transit compartment to induce expansion of cDCx cells. vmx represents the maximum rate of drug-induced expansion of cDCx cells, kmx denotes the drug concentration required to achieve 50% of vmx, and nx represents the Hill coefficient. kdeg2 is a second-order rate constant for apoptosis of cDCx cells. The term fDC1 only appears in the equation for cDC1 and denotes the fractional term that relates the second-order rate constant for apoptosis cDC1 to that of cDC2. In the figures, we show cDC1 and total DC, which is the sum of cDC1 and cDC2.

2.4. Overview of Model Calibration, Validation, and Simulation

As previously described, a mPBPK/PD model with expansion-enhanced TMDD in the central compartment was used to capture FLT3L-Fc preclinical data and project clinical PK/PD in healthy volunteers. The model-building strategy consists of eight steps (Figure 1D):
  • The minPBPK model was calibrated to FLT3L-Fc PK data from cynomolgus monkeys (0.1 to 10 mg/kg; preclinical study 1).
  • Simulations using the fitted model parameters were validated against a second cynomolgus monkey PK dataset (1 and 3 mg/kg; preclinical study 2).
  • Cynomolgus monkey PK parameters were translated to human PK parameters using the appropriate human non-specific clearance, physiological volumes, and flows to provide clinical PK predictions based solely on interspecies translation (the cyno-derived scenario).
  • In parallel, a separate nonlinear single-compartment PK model was built to capture FLT3L (CDX-301) clinical PK data across the range of doses evaluated in healthy volunteers (3 μg/kg to 75 μg/kg once daily subcutaneous (SC) injections for 5, 7 or 10 consecutive days; clinical study 1).
  • Using the CDX-301 PK model, parameters of human PD were calibrated to the published time-course data on cDC1 (clinical study 1) [23] and total DC counts (clinical study 2, part 1) [22].
  • Human PD simulations were validated using total DC expansion data in healthy volunteers following a range of FLT3L doses (10 μg/kg to 100 μg/kg SC injections once daily for 14 consecutive days; clinical study 2, part 2) [22].
  • Alternative clinical PK predictions of FLT3L-Fc were projected using clinical PK data from subjects treated with a single-dose FLT3 agonist Fc fusion protein molecule, GS-3583 (225 μg and 675 μg; clinical study 3) [40] (human-derived scenario).
  • Projected clinical PK from both cynomolgus monkey- and human-derived scenarios were used to predict cDC1 and total DC expansion for different doses of FLT3L-Fc.
The model was developed using the Simbiology® toolbox in MATLAB® (version R2018a, Mathworks Inc., Natick, MA, USA), and model parameterization/simulations were performed using the gQSPSim toolbox [50]. Particle swarm methodology [48] was used to calibrate model parameters to capture preclinical and clinical PK or PD measurements. The list of model parameters, rate constants, and variables for preclinical and clinical simulations is provided in Table S2.

2.5. Data Availability

The datasets generated and/or analyzed during the current study are available from the corresponding author upon reasonable request.

2.6. Code Availability

The code for execution of the model, along with the corresponding ordinary differential equations (ODEs), are available in the Supplementary Materials.

3. Results

3.1. Brief Description of the PK/PD Model of the FLT3L-Based Treatments

The minimal PBPK/PD model with expansion-enhanced TMDD was designed to capture FLT3L-Fc PK data in cynomolgus monkeys and Gilead’s GS-3583 PK data in healthy volunteers, including the effect of the drug-induced expansion of the target. The binding of FLT3L-Fc to one or two monomeric FLT3 receptors in circulation accounted for the observed TMDD. Homodimer binding drives the proliferation of immune cells, which leads to the expansion of the total FLT3 receptors on DCs and other bone marrow common progenitors, all of which were lumped into one total target capacity. The PK of CDX-301, which is administered as daily SC injections, was captured empirically using a single-compartment model with nonlinear clearance; the simple PK model was used for CDX-301 because its short half-life hinders the estimation and identification of TMDD-related parameters and non-specific clearance. After capturing the PK for both modalities, the CDX-301 PK was linked to an indirect PD response model via a series of transit compartments to capture the observed expansion of cDC1s and total DC cells. Exposure to the drug results in the first-order proliferation of the cells, whereas both first- and second-order apoptosis terms are used to capture the observed plateau in PD and the subsequent rapid decay once the drug clears. The model structure is depicted in Figure 1.

3.2. Calibration and Validation in Nonclinical Species to Establish Confidence in the Model

Surface Plasmon Resonance (SPR) binding analysis of FLT3L-Fc demonstrated comparable affinities to the recombinant human and cynomolgus monkey FLT3 receptors, which supports cynomolgus monkeys as a relevant test species for evaluation of FLT3L-Fc PK and its translation to humans. PK data in cynomolgus monkeys and the corresponding model fits are provided in Figure 2A. Plasma concentrations of FLT3L-Fc exhibited nonlinear PK within the tested dose range of 0.1 to 10 mg/kg. Anti-drug antibodies (ADAs) to FLT3L-Fc were detected by day 14 in all nine animals in preclinical study 1 (Figure 2A), and in all 12 animals in preclinical study 2 (Figure 2B). The dose-normalized area under the curve from days 0 to 21 (AUC0–21) for 0.1 to 10 mg/kg ranged from 7 ± 0.497 to 1660 ± 336 µg/mL·day, and corresponding total clearance (CL) ranged from 13.7 ± 0.901 to 4.92 ± 1.96 mL/day/kg. The dose-normalized AUC0–21 and CL were observed to be dose-dependent, suggesting nonlinear PK [50]. As shown in Figure 2A, for 10 mg/kg, the model successfully captured the Cmax and the linear slow elimination phase for ADA-negative PK samples but did not capture ADA-positive samples, particularly near the end of the study, as expected. The model also captured the nonlinear portion of the PK (for ADA-negative time points) well for lower doses at 0.1 and 1 mg/kg. The predictive capability of the model was verified by comparing simulations at 1 and 3 mg/kg doses with PK data from preclinical study 2 (Figure 2B). Measured concentrations of ADA-positive samples were again lower than predicted at later time points, likely reflecting the impact of ADA on exposure. However, the model was able to capture serum plasma levels in ADA-negative samples, which provided confidence in the predictive capabilities of the PK model for nonclinical species. The model was then used to simulate the dynamics of FLT3 receptor concentrations and receptor occupancy following FLT3L-Fc treatment in cynomolgus monkeys (Figure 3). Simulations suggest that increasing the dose from 0.01 to 10 mg/kg leads to increased total occupancy of FLT3 receptors in single-bound or double-bound complexes. At 1 mg/kg and higher, the model predicts near-complete receptor occupancy. Unlike single-bound receptor occupancy, which exhibited a dose-dependent increase, the percentage of double-bound receptor occupancy reached its peak at 0.1 mg/kg (83%) and declined following doses ranging from 1 to 10 mg/kg. This reflects a shift from double-bound complexes towards single-bound ones at higher drug levels.

3.3. The Model Reproduces Clinical PK and PD from Healthy Volunteers Treated with FLT3L (CDX-301)

A single-compartment PK model with nonlinear clearance was used to capture the observed PK of CDX-301 in healthy volunteers treated with daily SC injections at 3, 10, 25, and 75 μg/kg/day in clinical study 1 [23] (Figure 4A). CDX-301 PK exhibited a dose-linear increase in the maximum observed concentration (Cmax) and a dose-dependent increase in the area under the curve (AUC) [23]. The PD model parameters were fitted to the cDC1 time-course data from Anandasabapathy et al. (clinical study 1) [23] and total DC counts (equivalent to sum of cDC1s and cDC2s in the model) from Maraskovsky et al. (clinical study 2, part 1) [22]. A subset of total DC data from Maraskovsky et al. [22] (clinical study 2, part 2) was set aside for validation. The model simulations successfully reproduced clinical CDX-301 PK and captured the dynamics and dose-dependent saturable expansion of cDC1 and total DC counts reasonably well across both studies (Figure 4). Given that FLT3L-based therapies lead to a significant fold expansion of DCs, the group-specific baseline levels for cDC1 were incorporated into the model to improve the goodness of fit. In clinical study 1, CDX-301 daily injections from 3 to 75 μg/kg produced a dose-dependent expansion of cDC1s (Figure 4B). The model captured the dynamics of expansion for all groups. Although the model slightly overpredicted the peak of cDC1 expansion in the 25 μg/kg × 5 cohort and underpredicted the peak in the 25 μg/kg × 10 cohort, it successfully captured peak expansion for the upper (75 μg/kg × 5) and lower dose regimens (3 and 10 μg/kg × 5) (Figure 4B). Conversely, in clinical study 2, total DC expansion was comparable between the lowest dose (10 μg/kg/day) and the highest dose (100 μg/kg/day) following 14 days of daily treatments. Further, the total DC profile showed signs of saturation at around 10 days post initiation of FLT3L treatment (Figure 2D and Figure 4C). The model captured the rate of total DC expansion and the magnitude and timing of saturation following treatment with 100 μg/kg/day for 14 days (Figure 4C) and successfully predicted the saturation from 10–100 μg/kg (Figure 4D). These calibration and validation steps established confidence in the predictive capabilities for clinical PD projection.

3.4. The Model Projected Clinical PK and PD for FLT3L-Fc Treatment across Dose Levels

Cynomolgus monkey PK of FLT3L-Fc was used to project human PK by translating cyno-derived physiologically based PK parameters including lymph flow rate, interstitial fluid volume, and volume of distribution to human parameters and allometrically scaling CL, as previously published [42,43]. Two PK scenarios were simulated using different target-related parameters (i.e., FLT3 target capacity and proliferation-related parameters): (1) in the cyno-derived scenario, parameters were fixed based on the estimated parameters from cynomolgus monkey PK data or (2) in the human-derived scenario, parameters were adjusted to better match the PK of limited PK data of GS-3583 in healthy volunteers at two dose levels (Figure 5). The model simulated plasma exposures following multiple doses at levels ranging from 0.01–1 mg/kg for both scenarios. These PK profiles were then used to project the dynamics and level of cDC1 and total DC expansion using the PK/PD relationship established from the CDX-301 clinical data. Nonlinear PK was apparent in both scenarios at ≤0.1 mg/kg, and the loss of exposure due to TMDD at low dose levels was more pronounced in the cyno-derived scenario than the human-derived scenario because of the greater target concentration in the cyno-derived scenario (Figure 6A,D). As a result, the cyno-derived scenario required about a three-fold greater dose than the human-derived scenario to maximize the DC expansion effect (0.3 mg/kg vs. 0.1 mg/kg), causing cDC1 and total DC to plateau at approximately 70,000 and 1,000,000 cells/mL, respectively (Figure 6). Given that the human-derived scenario required a lower dose to induce a similar PD effect than the cyno-derived scenario, the human-derived scenario was chosen for FIH dose selection as the more conservative approach to achieve mPAD in healthy volunteers.

3.5. Rationale for FIH Dose Selection in Healthy Volunteers

A model-based FIH dose for FLT3L-Fc was selected based on an mPAD approach and supported with safety margins determined from FLT3L-Fc preclinical toxicity data. The criteria for the FIH dose were set to result in less than 20% of maximal potential peripheral blood cDC1 expansion. Assuming that the typical patient weighs 70 kg, a starting dose of 700 μg (0.01 mg/kg) was estimated to induce a peak expansion for peripheral blood cDC1 at 9500 cells/mL (about 8-fold expansion), which is approximately 13% of the maximum expansion level for cDC1 (Figure 6E). For the proposed FIH dose, the maximum total RO in the central compartment is estimated at 67%, with an average RO over 21 days of about 16% (Figure S1). In addition, the projected clinical exposure at the FIH dose level was estimated to be within the range of measured exposures following CDX-301 administration to healthy volunteers [23]. Furthermore, the exposures at the 3 mg/kg no-observed-adverse-effect level (NOAEL) from the cynomolgus monkey GLP-compliant toxicity study provided safety margins of more than 220-fold and more than 1500-fold over the clinical Cmax and AUC0–21, respectively, at the FIH dose [50]. The selected FIH dose also enables a full characterization of the PK/PD of FLT3L-Fc in healthy volunteers across a reasonable number of cohorts.
The model-based FIH dose selection and justification were communicated with the FDA in the regulatory submission, and the molecule is currently being tested in healthy volunteers (clinical trial number: 92655801, www.isrctn.com, accessed on 26 April 2024). The clinical simulations and adequate preclinical safety margins collectively supported the mPAD-based starting dose of 700 μg FLT3L-Fc IV for phase Ia clinical trials.

4. Discussion

Therapeutics that expand target-expressing cell populations are likely to exhibit expansion-enhanced TMDD, whereby higher doses drive greater expansion and thus increased total drug clearance [48,49]. The mechanistic minimal PBPK/PD model presented here integrates complex PK/PD relationships observed in preclinical and clinical studies of multiple FLT3L-based therapies. This model enabled the projection of clinical PK/PD profiles for FLT3L-Fc to inform clinical trial design, and was effectively utilized in our regulatory documents to propose the FIH dose selection for the clinical trial (clinical trial number: 92655801, www.isrctn.com, accessed on 26 April 2024). PK/PD was simulated based on the relationship observed from an extensive range of FLT3L (CDX-301) doses (3 μg/kg for 5 days to 100 μg/kg for 14 days) (Figure 4), encompassing a range of projected clinical exposures of FLT3L-Fc following FIH dosing. These clinical simulations provided quantitative guidance for the selection of 700 μg (0.01 mg/kg) as a safe starting dose in a phase Ia clinical trial (supported by adequate safety margins), while inducing an approximately 8-fold expansion of cDC1s in human subjects and avoiding exposure of human subjects to sub-pharmacological dose levels (clinical trial number: 92655801, www.isrctn.com, accessed on 26 April 2024). This approach and model-informed FIH dose selection were communicated with the FDA in the Investigator’s Brochure (IB) and the phase 1 protocol.
Nonlinear TMDD-driven PK behavior for FLT3L-Fc is expected in the phase Ia clinical trial based on the nonlinear PK observed in cynomolgus monkeys and the PK nonlinearity will likely be influenced by the level of target expansion observed in healthy volunteers. Cynomolgus monkeys are an appropriate animal model for clinical translation because human non-specific clearance can generally be well-predicted from observed PK in this species [44]. However, PK simulations translated from cynomolgus monkeys underpredicted GS-3583 PK data in healthy volunteers (Figure 5, dashed lines), suggesting that the target expression and other target-related parameters between humans and cynomolgus monkeys may not be equivalent. Reducing FLT3 receptor concentrations by 2.8-fold and adjusting the proliferation dynamics of FLT3 improved the human GS-3583 PK predictions, better capturing the initial and terminal time points in healthy volunteers (Table S2, Figure 5). In addition, these revised simulations led to a higher receptor occupancy and cDC1 and total DC expansion for a given FLT3L-Fc dose (Figures S1 and S2), providing a more conservative mPAD-based starting dose (0.01 mg/kg) compared to a starting dose selected based on the cyno-derived scenario (0.03 mg/kg).
The FIH starting dose of 700 μg (0.01 mg/kg) was selected to minimize the risk of potential adverse events in trial participants while ensuring that healthy volunteers are not subjected to sub-pharmacological exposures. The potential adverse effects of FLT3L-Fc are informed by findings from preclinical toxicity studies and the clinical safety profiles of previously studied FLT3L-based therapies. FLT3L-Fc was well-tolerated in cynomolgus monkeys up to 10 mg/kg q3w (with a NOAEL of 3 mg/kg identified in the GLP-compliant toxicity study), and the findings were generally consistent with the pharmacodynamic profile of FLT3L-Fc [50]. The clinical FLT3L-based molecules, CDX-301 and GS-3583, were also well-tolerated in healthy volunteers at comparable exposures to the proposed FIH starting dose [23,40]. A fraction of healthy volunteers (29%) receiving CDX-301 at a dose ≥ 25 μg/kg exhibited transient lymphadenopathy [23], suggesting the possibility of infiltration and expansion of histiocytes and leukocytes in lymphoid tissues. In FLT3L-Fc clinical studies, all study participants will be monitored for lymphadenopathy, splenomegaly, and organ enlargement. GS-3583 was well-tolerated as a single dose in healthy volunteers up to 2000 μg (or about 0.029 mg/kg assuming a 70 kg patient) with no reported serious grade 3 adverse effects [40]. The FIH dose of FLT3L-Fc is predicted to produce 8-fold cDC1 expansion, which is less than 20% of the maximum expansion observed in subjects treated with CDX-301. Furthermore, simulations suggest that DC counts are expected to return to baseline within the trial duration (Figure S2), minimizing the risk of extended immune cell expansion in trial participants. Collectively, preclinical data with FLT3L-Fc from the repeat-dose toxicity studies, available clinical data from FLT3L-based therapies, and model simulations support the proposed FIH dose for entry into the clinic.
The nonlinear PK observed in cynomolgus monkeys was reproduced using a minimal PBPK model that captured expansion-enhanced TMDD, but the model did not account for the impact of ADA on exposure due to lack of preclinical-to-clinical translation. ADAs were detected in all animals on or after day 14 and generally impacted FLT3L-Fc exposures in cynomolgus monkeys with high titers. However, the ADA response to biological therapeutics in nonclinical species has a generally low predictive power for estimating immunogenicity in humans [51,52,53]. Furthermore, no ADAs were detected on day 21 in healthy volunteers receiving daily treatments of CDX-301 [23]. Considering FLT3L-Fc is a fully human recombinant FLT3L directly fused to a native human Fc, it is not anticipated to result in significant immunogenicity in the clinical study.
Interspecies PK translation of non-specific clearance for large molecules is relatively well-studied in predicting exposure in humans [42,43,44], but, because the system-specific factors can vary between species, translating TMDD, let alone expansion-enhanced TMDD, and predicting PD response are more challenging and not as robust [54]. Interspecies differences in PD were apparent in the magnitude of response and expansion profiles of cDC1 and total DCs, whereas mice, cynomolgus monkeys, and humans treated with FLT3L-based therapies exhibited substantially different PK/PD relationships. In mice, a dose-dependent expansion of cDC1s and total DCs was observed in animals treated with a mouse surrogate FLT3L-Fc at doses of 0.1 to 30 mg/kg with no signs of PD saturation or bell-shaped response (see Figure 3 by Decalf et al. [41]). In cynomolgus monkeys treated with 0.1, 1, and 10 mg/kg of FLT3L-Fc, observations show characteristics of PD saturation and a bell-shaped response (see Figure 4 in Wu et al., 2024 [50]), presumably because higher doses favor formation of binary (single-bound) complexes rather than the active ternary (double-bound) complexes. However, the collective clinical PD data in healthy volunteers treated with FLT3L-based therapies suggests a dose-dependent saturable PD response with no apparent bell-shaped behavior. Due to the uncertainty in the translatability of PK/PD between species, we primarily relied on the PK/PD relationship observed in humans treated with CDX-301 and projected clinical PK profiles under two scenarios: (1) directly translating PK from cynomolgus monkeys treated with FLT3L-Fc using cyno-based target-related parameters, and (2) translating PK-related parameters from cynomolgus monkeys in addition to incorporating clinical PK data from GS-3583 to estimate target-related parameters. The projected PK profiles were then used to predict cDC1 and total DC expansion under the two scenarios. In this way, PK/PD projections for healthy volunteers integrated preclinical and clinical PK data from FLT3L-Fc and GS-3583, respectively, with the clinical PK/PD relationship from CDX-301.
Like other engineered cytokines, the mechanism of action of FLT3L-Fc and the resulting PK/PD is complex. Upon binding of the molecule to two FLT3 receptors and forming a ternary complex, the complex induces several downstream pathways that lead to cell differentiation and proliferation in different tissues. The complex is also internalized, which leads to target-mediated drug disposition. Expansion of immune cells due to FLT3L-Fc treatment creates feedback that increases the total FLT3 target capacity and hence further augments TMDD for the drug. In this work, a decoupled approach [48,49] was used in model development to capture PK/PD data [48]. When employing the decoupled approach, it is crucial to utilize relevant PK/PD information, given that translation across species might be limited. Both coupled [49] and decoupled approaches have been used in modeling engineered cytokines. Coupling the PK and PD via FLT3 receptors could potentially provide more insight into mechanisms of drug activity. However, it requires extensive information on key system parameters (e.g., cell numbers, tissues of origin, receptor expression per cell type, and internalization rates of free and bound receptors, as well as the required threshold for the number of ligand-receptor complexes per cell type to induce the downstream pathways), complicating the model development by requiring elaborate, biological systems-based model structures used in quantitative system pharmacology (QSP) models. Further, the model predictions stand the risk of being inaccurate if the underlying data are not robust. Note that species differences across these parameters and their combined effect could lead to significant PK/PD differences and likely explains the complex PK and PD relationships observed in mice, cynomolgus monkeys, and humans. A decoupled approach is useful as it aims to capture PK and PD separately; in this modeling approach, an assumption is made that a total target capacity in the circulation can mimic the effects of FLT3 receptor expression on different cells in different tissues. Additionally, this total FLT3 target capacity can proliferate to mimic the effect of cell expansion in response to FLT3L-Fc treatment. Once the model captures the PK across different dose levels, the exposure can be used as the key driver of expansion for any cell type of interest, e.g., cDC1s and total DCs. This approach simplifies the biology, which involves the expansion of precursors in the bone marrow and differentiation to several immune cell populations. This model has several key practical advantages: the process of model calibration is more straightforward, and the model is more amenable to including other cell types without a significant need to re-calibrate the other aspects of the model. However, future QSP modeling could be utilized to properly account for key system parameters and improve the understanding of the role of FLT3L-based treatments in biological systems.
This work presents a model constructed using fixed physiological parameters with well-characterized values, effectively reducing concerns about parameter uncertainty. The key model parameters influencing PK and PD include target expression and PK parameters which were suitably fitted using data from cynomolgus monkeys. Validation was successfully performed using an independent dataset, which reinforced confidence in the generalizability and robustness to variations in parameter values. Additional model fitting and validation were performed using sparse and digitized clinical PD data from CDX-301. However, the model parameters can easily be refined to capture additional clinical data from FLT3L-Fc once they become available. FLT3L-Fc has the potential to improve and enhance anti-tumor immunity in combination with other treatments such as DC maturation agents and checkpoint inhibitors. Therefore, future uses of the model could include informing combination strategies to evaluate DC activation and anti-tumor efficacy in patients.

Supplementary Materials

The following supporting information can be downloaded at https://www.mdpi.com/article/10.3390/pharmaceutics16050660/s1, Supplemental ODEs and Repeated Assignments; Supplemental Instructions of gQSPsim, Table S1: Overview of clinical and preclinical studies used for model calibration/validation, Table S2: Model Parameters, Figure S1: Projected time profiles of single-bound, double-bound and total receptor occupancy in healthy volunteers treated with a single dose of FLT3L-Fc, Figure S2: Projected time profiles of clinical PK, cDC1, and total DC cells in healthy volunteers treated with a single dose of FLT3L-Fc. References [55,56] are cited in the Supplementary Materials.

Author Contributions

Conceptualization, methodology, software, validation, formal analysis, investigation, writing—original draft, writing—review and editing, visualization, supervision, and project administration, I.H. Methodology, software, validation, formal analysis, investigation, data curation, writing—original draft, writing—review and editing, and visualization, B.F. Investigation, formal analysis, resources, data curation, and writing—review and editing, J.G. Investigation, resources, data curation, and formal analysis, J.D. Investigation, resources, data curation, and formal analysis, M.K. Conceptualization, methodology, and data curation, M.O. Investigation, resources, data curation, and formal analysis, T.W.B. Investigation, resources, data curation, formal analysis, writing—reviewing and editing, C.M. Investigation, resources, writing—review and editing, project administration, G.K.R.; Writing—review and editing, and supervision, K.G. Resources, writing—review and editing, supervision, project administration, and funding acquisition, A.V.K. Conceptualization, methodology, resources, investigation, writing—original draft, writing—review and editing, and supervision, S.R. All authors have read and agreed to the published version of the manuscript.

Funding

The work was funded by Genentech, Inc.

Institutional Review Board Statement

The animal study protocol was approved by the Institutional Review Board of Charles River Laboratories (protocol code 20239581, 30 March 2020; protocol code 5003546, 7 December 2020).

Informed Consent Statement

Not applicable.

Data Availability Statement

Data is available upon request from the authors.

Acknowledgments

The authors would like to thank Kai Connie Wu, Tien Nguyen, Tanja Zabka, and Adeyemi Adedeji for their support in assessing the preclinical safety data. We thank Xiangdan Wang, and Laeticia Comps-Agrar for quantifying the binding affinity of FLT3L-Fc against human and cynomolgus monkey FLT3 receptors, as well as Pamela Chan, Jenny Li, Chris Lowe, and Weili Yan for preclinical PK and ADA support.

Conflicts of Interest

I.H., B.F., J.G., J.D., M.K., M.O., T.W.B., C.M., G.K.R., K.G. and A.V.K. are all employees of Genentech, Inc. S.R. was an employee of Genentech, Inc. at the time the studies were conducted and is currently an employee at Generate Biomedicines. All authors may be shareholders in F. Hoffmann-La Roche AG.

References

  1. Reck, M.; Rodríguez-Abreu, D.; Robinson, A.G.; Hui, R.; Csőszi, T.; Fülöp, A.; Gottfried, M.; Peled, N.; Tafreshi, A.; Cuffe, S.; et al. Pembrolizumab versus Chemotherapy for PD-L1-Positive Non-Small-Cell Lung Cancer. N. Engl. J. Med. 2016, 375, 1823–1833. [Google Scholar] [CrossRef]
  2. Robert, C.; Long, G.V.; Brady, B.; Dutriaux, C.; Maio, M.; Mortier, L.; Hassel, J.C.; Rutkowski, P.; McNeil, C.; Kalinka-Warzocha, E.; et al. Nivolumab in previously untreated melanoma without BRAF mutation. N. Engl. J. Med. 2015, 372, 320–330. [Google Scholar] [CrossRef]
  3. Rosenberg, A.S.; Worobec, A. A Risk-Based Approach to Immunogenicity Concerns of Therapeutic Protein Products, Part 1: Considering Consequences of the Immune Response to a Protein. BioPharm Int. 2004, 17, 22–26. [Google Scholar]
  4. Bachem, A.; Güttler, S.; Hartung, E.; Ebstein, F.; Schaefer, M.; Tannert, A.; Salama, A.; Movassaghi, K.; Opitz, C.; Mages, H.W.; et al. Superior antigen cross-presentation and XCR1 expression define human CD11c+CD141+ cells as homologues of mouse CD8+ dendritic cells. J. Exp. Med. 2010, 207, 1273–1281. [Google Scholar] [CrossRef]
  5. Jongbloed, S.L.; Kassianos, A.J.; McDonald, K.J.; Clark, G.J.; Ju, X.; Angel, C.E.; Chen, C.J.; Dunbar, P.R.; Wadley, R.B.; Jeet, V.; et al. Human CD141+ (BDCA-3)+ dendritic cells (DCs) represent a unique myeloid DC subset that cross-presents necrotic cell antigens. J. Exp. Med. 2010, 207, 1247–1260. [Google Scholar] [CrossRef] [PubMed]
  6. Poulin, L.F.; Salio, M.; Griessinger, E.; Anjos-Afonso, F.; Craciun, L.; Chen, J.L.; Keller, A.M.; Joffre, O.; Zelenay, S.; Nye, E.; et al. Characterization of human DNGR-1+ BDCA3+ leukocytes as putative equivalents of mouse CD8alpha+ dendritic cells. J. Exp. Med. 2010, 207, 1261–1271. [Google Scholar] [CrossRef]
  7. Haniffa, M.; Shin, A.; Bigley, V.; McGovern, N.; Teo, P.; See, P.; Wasan, P.S.; Wang, X.N.; Malinarich, F.; Malleret, B.; et al. Human tissues contain CD141hi cross-presenting dendritic cells with functional homology to mouse CD103+ nonlymphoid dendritic cells. Immunity 2012, 37, 60–73. [Google Scholar] [CrossRef] [PubMed]
  8. Salmon, H.; Idoyaga, J.; Rahman, A.; Leboeuf, M.; Remark, R.; Jordan, S.; Casanova-Acebes, M.; Khudoynazarova, M.; Agudo, J.; Tung, N.; et al. Expansion and Activation of CD103(+) Dendritic Cell Progenitors at the Tumor Site Enhances Tumor Responses to Therapeutic PD-L1 and BRAF Inhibition. Immunity 2016, 44, 924–938. [Google Scholar] [CrossRef]
  9. Sánchez-Paulete, A.R.; Cueto, F.J.; Martínez-López, M.; Labiano, S.; Morales-Kastresana, A.; Rodríguez-Ruiz, M.E.; Jure-Kunkel, M.; Azpilikueta, A.; Aznar, M.A.; Quetglas, J.I.; et al. Cancer Immunotherapy with Immunomodulatory Anti-CD137 and Anti-PD-1 Monoclonal Antibodies Requires BATF3-Dependent Dendritic Cells. Cancer Discov. 2016, 6, 71–79. [Google Scholar] [CrossRef] [PubMed]
  10. Spranger, S.; Dai, D.; Horton, B.; Gajewski, T.F. Tumor-Residing Batf3 Dendritic Cells Are Required for Effector T Cell Trafficking and Adoptive T Cell Therapy. Cancer Cell 2017, 31, 711–723.e4. [Google Scholar] [CrossRef]
  11. Lee, J.H.; Choi, S.Y.; Jung, N.C.; Song, J.Y.; Seo, H.G.; Lee, H.S.; Lim, D.S. The Effect of the Tumor Microenvironment and Tumor-Derived Metabolites on Dendritic Cell Function. J. Cancer 2020, 11, 769–775. [Google Scholar] [CrossRef]
  12. Wculek, S.K.; Cueto, F.J.; Mujal, A.M.; Melero, I.; Krummel, M.F.; Sancho, D. Dendritic cells in cancer immunology and immunotherapy. Nat. Rev. Immunol. 2020, 20, 7–24. [Google Scholar] [CrossRef]
  13. Broz, M.L.; Binnewies, M.; Boldajipour, B.; Nelson, A.E.; Pollack, J.L.; Erle, D.J.; Barczak, A.; Rosenblum, M.D.; Daud, A.; Barber, D.L.; et al. Dissecting the tumor myeloid compartment reveals rare activating antigen-presenting cells critical for T cell immunity. Cancer Cell 2014, 26, 638–652. [Google Scholar] [CrossRef]
  14. Roberts, E.W.; Broz, M.L.; Binnewies, M.; Headley, M.B.; Nelson, A.E.; Wolf, D.M.; Kaisho, T.; Bogunovic, D.; Bhardwaj, N.; Krummel, M.F. Critical Role for CD103(+)/CD141(+) Dendritic Cells Bearing CCR7 for Tumor Antigen Trafficking and Priming of T Cell Immunity in Melanoma. Cancer Cell 2016, 30, 324–336. [Google Scholar] [CrossRef]
  15. Mayoux, M.; Roller, A.; Pulko, V.; Sammicheli, S.; Chen, S.; Sum, E.; Jost, C.; Fransen, M.F.; Buser, R.B.; Kowanetz, M.; et al. Dendritic cells dictate responses to PD-L1 blockade cancer immunotherapy. Sci. Transl. Med. 2020, 12, eaav7431. [Google Scholar] [CrossRef]
  16. Small, D.; Levenstein, M.; Kim, E.; Carow, C.; Amin, S.; Rockwell, P.; Witte, L.; Burrow, C.; Ratajczak, M.Z.; Gewirtz, A.M.; et al. STK-1, the human homolog of Flk-2/Flt-3, is selectively expressed in CD34+ human bone marrow cells and is involved in the proliferation of early progenitor/stem cells. Proc. Natl. Acad. Sci. USA 1994, 91, 459–463. [Google Scholar] [CrossRef]
  17. Tsapogas, P.; Swee, L.K.; Nusser, A.; Nuber, N.; Kreuzaler, M.; Capoferri, G.; Rolink, H.; Ceredig, R.; Rolink, A. In vivo evidence for an instructive role of fms-like tyrosine kinase-3 (FLT3) ligand in hematopoietic development. Haematologica 2014, 99, 638–646. [Google Scholar] [CrossRef]
  18. Adolfsson, J.; Borge, O.J.; Bryder, D.; Theilgaard-Mönch, K.; Astrand-Grundström, I.; Sitnicka, E.; Sasaki, Y.; Jacobsen, S.E. Upregulation of Flt3 expression within the bone marrow Lin(-)Sca1(+)c-kit(+) stem cell compartment is accompanied by loss of self-renewal capacity. Immunity 2001, 15, 659–669. [Google Scholar] [CrossRef]
  19. Christensen, J.L.; Weissman, I.L. Flk-2 is a marker in hematopoietic stem cell differentiation: A simple method to isolate long-term stem cells. Proc. Natl. Acad. Sci. USA 2001, 98, 14541–14546. [Google Scholar] [CrossRef]
  20. Higano, C.S.; Vogelzang, N.J.; Sosman, J.A.; Feng, A.; Caron, D.; Small, E.J. Safety and biological activity of repeated doses of recombinant human Flt3 ligand in patients with bone scan-negative hormone-refractory prostate cancer. Clin. Cancer Res. 2004, 10, 1219–1225. [Google Scholar] [CrossRef]
  21. Maraskovsky, E.; Brasel, K.; Teepe, M.; Roux, E.R.; Lyman, S.D.; Shortman, K.; McKenna, H.J. Dramatic increase in the numbers of functionally mature dendritic cells in Flt3 ligand-treated mice: Multiple dendritic cell subpopulations identified. J. Exp. Med. 1996, 184, 1953–1962. [Google Scholar] [CrossRef]
  22. Maraskovsky, E.; Daro, E.; Roux, E.; Teepe, M.; Maliszewski, C.R.; Hoek, J.; Caron, D.; Lebsack, M.E.; McKenna, H.J. In vivo generation of human dendritic cell subsets by Flt3 ligand. Blood 2000, 96, 878–884. [Google Scholar] [CrossRef]
  23. Anandasabapathy, N.; Breton, G.; Hurley, A.; Caskey, M.; Trumpfheller, C.; Sarma, P.; Pring, J.; Pack, M.; Buckley, N.; Matei, I.; et al. Efficacy and safety of CDX-301, recombinant human Flt3L, at expanding dendritic cells and hematopoietic stem cells in healthy human volunteers. Bone Marrow Transpl. 2015, 50, 924–930. [Google Scholar] [CrossRef]
  24. Pulendran, B.; Banchereau, J.; Burkeholder, S.; Kraus, E.; Guinet, E.; Chalouni, C.; Caron, D.; Maliszewski, C.; Davoust, J.; Fay, J.; et al. Flt3-ligand and granulocyte colony-stimulating factor mobilize distinct human dendritic cell subsets in vivo. J. Immunol. 2000, 165, 566–572. [Google Scholar] [CrossRef]
  25. McKenna, H.J. Role of hematopoietic growth factors/flt3 ligand in expansion and regulation of dendritic cells. Curr. Opin. Hematol. 2001, 8, 149–154. [Google Scholar] [CrossRef]
  26. O’Keeffe, M.; Hochrein, H.; Vremec, D.; Pooley, J.; Evans, R.; Woulfe, S.; Shortman, K. Effects of administration of progenipoietin 1, Flt-3 ligand, granulocyte colony-stimulating factor, and pegylated granulocyte-macrophage colony-stimulating factor on dendritic cell subsets in mice. Blood 2002, 99, 2122–2130. [Google Scholar] [CrossRef]
  27. Karsunky, H.; Merad, M.; Cozzio, A.; Weissman, I.L.; Manz, M.G. Flt3 ligand regulates dendritic cell development from Flt3+ lymphoid and myeloid-committed progenitors to Flt3+ dendritic cells in vivo. J. Exp. Med. 2003, 198, 305–313. [Google Scholar] [CrossRef]
  28. Böttcher, J.P.; Reis e Sousa, C. The Role of Type 1 Conventional Dendritic Cells in Cancer Immunity. Trends Cancer 2018, 4, 784–792. [Google Scholar] [CrossRef]
  29. Fu, C.; Jiang, A. Dendritic Cells and CD8 T Cell Immunity in Tumor Microenvironment. Front. Immunol. 2018, 9, 3059. [Google Scholar] [CrossRef]
  30. Balan, S.; Saxena, M.; Bhardwaj, N. Dendritic cell subsets and locations. Int. Rev. Cell Mol. Biol. 2019, 348, 1–68. [Google Scholar] [CrossRef]
  31. Wculek, S.K.; Amores-Iniesta, J.; Conde-Garrosa, R.; Khouili, S.C.; Melero, I.; Sancho, D. Effective cancer immunotherapy by natural mouse conventional type-1 dendritic cells bearing dead tumor antigen. J. ImmunoTherapy Cancer 2019, 7, 100. [Google Scholar] [CrossRef]
  32. Demaria, S.; Ng, B.; Devitt, M.L.; Babb, J.S.; Kawashima, N.; Liebes, L.; Formenti, S.C. Ionizing radiation inhibition of distant untreated tumors (abscopal effect) is immune mediated. Int. J. Radiat. Oncol. Biol. Phys. 2004, 58, 862–870. [Google Scholar] [CrossRef]
  33. Miao, X.; Chen, Y.; Hao, K.; Zheng, M.; Chen, B.; Li, K.; Wang, Y.; Zhang, W.; Zhang, Y.; Mou, X.; et al. CD103(+) Cell Growth Factor Flt3L Enhances the Efficacy of Immune Checkpoint Blockades in Murine Glioblastoma Model. Oncol. Res. 2018, 26, 173–182. [Google Scholar] [CrossRef]
  34. Hammerich, L.; Marron, T.U.; Upadhyay, R.; Svensson-Arvelund, J.; Dhainaut, M.; Hussein, S.; Zhan, Y.; Ostrowski, D.; Yellin, M.; Marsh, H.; et al. Systemic clinical tumor regressions and potentiation of PD1 blockade with in situ vaccination. Nat. Med. 2019, 25, 814–824. [Google Scholar] [CrossRef]
  35. Oba, T.; Long, M.D.; Keler, T.; Marsh, H.C.; Minderman, H.; Abrams, S.I.; Liu, S.; Ito, F. Overcoming primary and acquired resistance to anti-PD-L1 therapy by induction and activation of tumor-residing cDC1s. Nat. Commun. 2020, 11, 5415. [Google Scholar] [CrossRef]
  36. Ohri, N.; Halmos, B.; Garg, M.; Levsky, J.M.; Cheng, H.; Gucalp, R.A.; Bodner, W.R.; Kabarriti, R.; Berkowitz, A.; Yellin, M.J.; et al. FLT3 ligand (CDX-301) and stereotactic radiotherapy for advanced non-small cell lung cancer. J. Clin. Oncol. 2020, 38, 9618. [Google Scholar] [CrossRef]
  37. Lee, Y.S.; O’Brien, L.J.; Walpole, C.M.; Pearson, F.E.; Leal-Rojas, I.M.; Masterman, K.A.; Atkinson, V.; Barbour, A.; Radford, K.J. Human CD141(+) dendritic cells (cDC1) are impaired in patients with advanced melanoma but can be targeted to enhance anti-PD-1 in a humanized mouse model. J. Immunotherapy Cancer 2021, 9, e001963. [Google Scholar] [CrossRef]
  38. Saxena, M.; Sabado, R.L.; La Mar, M.; Mohri, H.; Salazar, A.M.; Dong, H.; Correa Da Rosa, J.; Markowitz, M.; Bhardwaj, N.; Miller, E. Poly-ICLC, a TLR3 Agonist, Induces Transient Innate Immune Responses in Patients with Treated HIV-Infection: A Randomized Double-Blinded Placebo Controlled Trial. Front. Immunol. 2019, 10, 725. [Google Scholar] [CrossRef]
  39. Bhardwaj, N.; Pavlick, A.C.; Ernstoff, M.S.; Hanks, B.A.; Albertini, M.R.; Luke, J.J.; Yellin, M.J.; Keler, T.; Davis, T.A.; Crocker, A.; et al. A Phase II Randomized Study of CDX-1401, a Dendritic Cell Targeting NY-ESO-1 Vaccine, in Patients with Malignant Melanoma Pre-Treated with Recombinant CDX-301, a Recombinant Human Flt3 Ligand. J. Clin. Oncol. 2016, 34, 9589. [Google Scholar] [CrossRef]
  40. Rajakumaraswamy, N.; Dauki, A.; Kuhne, M.R.; Trowe, T.; He, J.; Clarke, C.; Carr, B.; Worth, A.; Vashishtha, A.; Schwabe, C.; et al. GS-3583, a novel FLT3 agonist Fc fusion protein, to expand conventional dendritic cells in healthy volunteers. J. Clin. Oncol. 2021, 39, 2559. [Google Scholar] [CrossRef]
  41. Decalf, J.; Toy, E.; He, D.; Kenkre, R.; Berkley, A.; Wong, D.; Kwong, M.; Kee, Y.S.; Sun, Y.; Myneni, S.; et al. Effectorless Fc-fusion improves FLT3L drug-like properties for cancer immunotherapy combinations. Lancet 2024. manuscript under review. [Google Scholar]
  42. Cao, Y.; Balthasar, J.P.; Jusko, W.J. Second-generation minimal physiologically-based pharmacokinetic model for monoclonal antibodies. J. Pharmacokinet. Pharmacodyn. 2013, 40, 597–607. [Google Scholar] [CrossRef]
  43. Cao, Y.; Jusko, W.J. Incorporating target-mediated drug disposition in a minimal physiologically-based pharmacokinetic model for monoclonal antibodies. J. Pharmacokinet. Pharmacodyn. 2014, 41, 375–387. [Google Scholar] [CrossRef]
  44. Deng, R.; Iyer, S.; Theil, F.P.; Mortensen, D.L.; Fielder, P.J.; Prabhu, S. Projecting human pharmacokinetics of therapeutic antibodies from nonclinical data: What have we learned? MAbs 2011, 3, 61–66. [Google Scholar] [CrossRef]
  45. Zhao, J.; Cao, Y.; Jusko, W.J. Across-Species Scaling of Monoclonal Antibody Pharmacokinetics Using a Minimal PBPK Model. Pharm. Res. 2015, 32, 3269–3281. [Google Scholar] [CrossRef]
  46. Sarin, H. Physiologic upper limits of pore size of different blood capillary types and another perspective on the dual pore theory of microvascular permeability. J. Angiogenesis Res. 2010, 2, 14. [Google Scholar] [CrossRef]
  47. Graddis, T.J.; Brasel, K.; Friend, D.; Srinivasan, S.; Wee, S.; Lyman, S.D.; March, C.J.; McGrew, J.T. Structure-function analysis of FLT3 ligand-FLT3 receptor interactions using a rapid functional screen. J. Biol. Chem. 1998, 273, 17626–17633. [Google Scholar] [CrossRef]
  48. Ribba, B.; Boetsch, C.; Nayak, T.; Grimm, H.P.; Charo, J.; Evers, S.; Klein, C.; Tessier, J.; Charoin, J.E.; Phipps, A.; et al. Prediction of the Optimal Dosing Regimen Using a Mathematical Model of Tumor Uptake for Immunocytokine-Based Cancer Immunotherapy. Clin. Cancer Res. 2018, 24, 3325–3333. [Google Scholar] [CrossRef]
  49. Lu, D.; Yadav, R.; Holder, P.; Chiang, E.; Sanjabi, S.; Poon, V.; Bernett, M.; Varma, R.; Liu, K.; Leung, I.; et al. Complex PK-PD of an engineered IL-15/IL-15Rα–Fc fusion protein in cynomolgus monkeys: QSP modeling of lymphocyte dynamics. Eur. J. Pharm. Sci. 2023, 186, 106450. [Google Scholar] [CrossRef]
  50. Wu, K.C.; Adedeji, A.O.; Zabka, T.S.; Hosseini, I.; Kenkre, R.; Getz, J.A.; Nguyen, T.; Decalf, J.; Bainbridge, T.W.; Chilton, J.A.; et al. Nonclinical pharmacokinetics, pharmacodynamics and safety assessment of a FLT3L-Fc molecule for cancer immunotherapy. Toxicol. Appl. Pharmacol. 2024, 483, 116837. [Google Scholar] [CrossRef]
  51. Leach, M.W.; Rottman, J.B.; Hock, M.B.; Finco, D.; Rojko, J.L.; Beyer, J.C. Immunogenicity/hypersensitivity of biologics. Toxicol. Pathol. 2014, 42, 293–300. [Google Scholar] [CrossRef] [PubMed]
  52. Rojko, J.L.; Evans, M.G.; Price, S.A.; Han, B.; Waine, G.; DeWitte, M.; Haynes, J.; Freimark, B.; Martin, P.; Raymond, J.T.; et al. Formation, clearance, deposition, pathogenicity, and identification of biopharmaceutical-related immune complexes: Review and case studies. Toxicol. Pathol. 2014, 42, 725–764. [Google Scholar] [CrossRef] [PubMed]
  53. van Meer, P.J.K.; Kooijman, M.; Brinks, V.; Gispen-de Wied, C.C.; Silva-Lima, B.; Moors, E.H.M.; Schellekens, H. Immunogenicity of mAbs in non-human primates during nonclinical safety assessment. mAbs 2013, 5, 810–816. [Google Scholar] [CrossRef] [PubMed]
  54. Mager, D.E.; Woo, S.; Jusko, W.J. Scaling pharmacodynamics from in vitro and preclinical animal studies to humans. Drug Metab. Pharmacokinet. 2009, 24, 16–24. [Google Scholar] [CrossRef] [PubMed]
  55. Satyamitra, M.; Cary, L.; Dunn, D.; Holmes-Hampton, G.P.; Thomas, L.J.; Ghosh, S.P. CDX-301: A novel medical countermeasure for hematopoietic acute radiation syndrome in mice. Sci. Rep. 2020, 10, 1757. [Google Scholar] [CrossRef]
  56. Verstraete, K.; Vandriessche, G.; Januar, M.; Elegheert, J.; Shkumatov, A.V.; Desfosses, A.; Craenenbroeck, K.V.; Svergun, D.I.; Gutsche, I.; Vergauwen, B.; et al. Structural insights into the extracellular assembly of the hematopoietic Flt3 signaling complex. Blood 2011, 118, 60–68. [Google Scholar] [CrossRef]
Figure 1. FLT3L-Fc minimal physiologically based PK/PD model and workflow. The minimal physiologically based pharmacokinetic (minPBPK) model and workflow describe the process of describing the in vivo dynamics of cDC1 and cDC2 following the interaction of FLT3L with the FLT3 receptors. (A) FLT3L-Fc and GS-3583 can bind to one or two FLT3 receptors in the central compartment to capture the impact of TMDD; the homodimer binding induces proliferation of different immune cells and leads to the expansion of the total target capacity of FLT3 receptors. (B) A one-compartment model with a subcutaneous compartment and a nonlinear clearance term describes the PK of FLT3L (CDX-301). (C) Free drug in the central compartment is linked to an indirect PD response model via a series of transit compartments to capture the observed expansion of cDC1 and cDC2 cells; dotted lines represent the indirect impact of the drug on cDC1 and cDC2 dynamics. (D) The model-building and dose projection workflow consist of eight steps described in the Section 2.
Figure 1. FLT3L-Fc minimal physiologically based PK/PD model and workflow. The minimal physiologically based pharmacokinetic (minPBPK) model and workflow describe the process of describing the in vivo dynamics of cDC1 and cDC2 following the interaction of FLT3L with the FLT3 receptors. (A) FLT3L-Fc and GS-3583 can bind to one or two FLT3 receptors in the central compartment to capture the impact of TMDD; the homodimer binding induces proliferation of different immune cells and leads to the expansion of the total target capacity of FLT3 receptors. (B) A one-compartment model with a subcutaneous compartment and a nonlinear clearance term describes the PK of FLT3L (CDX-301). (C) Free drug in the central compartment is linked to an indirect PD response model via a series of transit compartments to capture the observed expansion of cDC1 and cDC2 cells; dotted lines represent the indirect impact of the drug on cDC1 and cDC2 dynamics. (D) The model-building and dose projection workflow consist of eight steps described in the Section 2.
Pharmaceutics 16 00660 g001
Figure 2. Time-course data and model fits for FLT3L-Fc PK in cynomolgus monkeys. Individual cynomolgus monkeys (n = 3) were treated with a single dose of 0.1 mg/kg, or two doses of 1 mg/kg, or 10 mg/kg with a three-week interval between doses in preclinical study 1. Individual cynomolgus monkeys (n = 6) were treated with two doses of 1 mg/kg, or 3 mg/kg with a three-week interval between doses in preclinical study 2. (A) Drug concentrations of FLT3L-Fc from preclinical study 1 was used to calibrate the cynomolgus monkey PK parameters in the model. (B) PK from preclinical study 2 served as a validation dataset for the model predictions. Open squares and solid dots represent individual PK samples that were ADA-positive and ADA-negative, respectively. LOQ: limit of quantification.
Figure 2. Time-course data and model fits for FLT3L-Fc PK in cynomolgus monkeys. Individual cynomolgus monkeys (n = 3) were treated with a single dose of 0.1 mg/kg, or two doses of 1 mg/kg, or 10 mg/kg with a three-week interval between doses in preclinical study 1. Individual cynomolgus monkeys (n = 6) were treated with two doses of 1 mg/kg, or 3 mg/kg with a three-week interval between doses in preclinical study 2. (A) Drug concentrations of FLT3L-Fc from preclinical study 1 was used to calibrate the cynomolgus monkey PK parameters in the model. (B) PK from preclinical study 2 served as a validation dataset for the model predictions. Open squares and solid dots represent individual PK samples that were ADA-positive and ADA-negative, respectively. LOQ: limit of quantification.
Pharmaceutics 16 00660 g002
Figure 3. Projected target expression and receptor occupancy in cynomolgus monkeys following FLT3L-Fc treatment. (A) The total estimated target expression in cynomolgus monkeys following a single FLT3L-Fc dose (0.01, 0.03, and 0.1 mg/kg) or two FLT3L-Fc doses (1 and 10 mg/kg) administered three weeks apart is shown. (B) The percent of total FLT3 receptor occupancy (RO), (C) percentage of FLT3 receptors occupied in a single-bound complex, and (D) percentage of FLT3 receptors occupied in a double-bound complex are also shown.
Figure 3. Projected target expression and receptor occupancy in cynomolgus monkeys following FLT3L-Fc treatment. (A) The total estimated target expression in cynomolgus monkeys following a single FLT3L-Fc dose (0.01, 0.03, and 0.1 mg/kg) or two FLT3L-Fc doses (1 and 10 mg/kg) administered three weeks apart is shown. (B) The percent of total FLT3 receptor occupancy (RO), (C) percentage of FLT3 receptors occupied in a single-bound complex, and (D) percentage of FLT3 receptors occupied in a double-bound complex are also shown.
Pharmaceutics 16 00660 g003
Figure 4. Clinical data and model outputs for PK, cDC1, and total DC counts in healthy human volunteers following treatment with FLT3L (CDX-301). PK/PD data from healthy human volunteers treated with FLT3L (CDX-301, Celldex Therapeutics, Inc., Hampton, NJ) was used to calibrate the model to capture the dynamics of cDC1 and total DC counts per mL following different FLT3L exposure levels. (A) Time-course CDX-301 concentration data were measured in subjects following the administration of daily SC injections of CDX-301. Data were digitized from clinical study 1 [23]. (B) Clinical PK of CDX-301 were linked to a PD model and parameters were fitted to digitized human cDC1 time-course data from clinical study 1 [23] and (C) total DC counts digitized from clinical study 2 [22]. (D) The model predictions were tested against total DC counts measured 15 days post FLT3L treatment from clinical study 2 (22). Dots show digitized data from clinical studies 1 and 2.
Figure 4. Clinical data and model outputs for PK, cDC1, and total DC counts in healthy human volunteers following treatment with FLT3L (CDX-301). PK/PD data from healthy human volunteers treated with FLT3L (CDX-301, Celldex Therapeutics, Inc., Hampton, NJ) was used to calibrate the model to capture the dynamics of cDC1 and total DC counts per mL following different FLT3L exposure levels. (A) Time-course CDX-301 concentration data were measured in subjects following the administration of daily SC injections of CDX-301. Data were digitized from clinical study 1 [23]. (B) Clinical PK of CDX-301 were linked to a PD model and parameters were fitted to digitized human cDC1 time-course data from clinical study 1 [23] and (C) total DC counts digitized from clinical study 2 [22]. (D) The model predictions were tested against total DC counts measured 15 days post FLT3L treatment from clinical study 2 (22). Dots show digitized data from clinical studies 1 and 2.
Pharmaceutics 16 00660 g004
Figure 5. Clinical PK data and simulated PK profiles in healthy human volunteers following treatment with GS-3583. The FLT3L-Fc PK parameters were translated from parameters estimated using cynomolgus monkey PK data. FLT3 target-related parameters were estimated from cynomolgus monkeys (dashed line) or refined to improve the fit to human time-course data for the FLT3 agonist Fc fusion protein molecule (GS-3583, Gilead Sciences, Inc., Foster City, CA, solid line). The dots represent the mean concentration-time profiles of GS-3583 following a single IV administration of 225 and 675 μg to healthy human volunteers.
Figure 5. Clinical PK data and simulated PK profiles in healthy human volunteers following treatment with GS-3583. The FLT3L-Fc PK parameters were translated from parameters estimated using cynomolgus monkey PK data. FLT3 target-related parameters were estimated from cynomolgus monkeys (dashed line) or refined to improve the fit to human time-course data for the FLT3 agonist Fc fusion protein molecule (GS-3583, Gilead Sciences, Inc., Foster City, CA, solid line). The dots represent the mean concentration-time profiles of GS-3583 following a single IV administration of 225 and 675 μg to healthy human volunteers.
Pharmaceutics 16 00660 g005
Figure 6. Projected time profiles of clinical PK, cDC1, and total DC cells in healthy volunteers treated with FLT3L-Fc. Simulations were performed to predict clinical FLT3L-Fc PK/PD profiles and compare cDC1 and total DC expansion under both cyno-derived and human-derived scenarios for a q3w dosing regimen. (AC) In the first row, FLT3 target-related parameters were estimated from cynomolgus monkey PK, and (DF) in the second row, FLT3 target-related parameters were estimated using human PK data from the FLT3 agonist Fc fusion protein molecule, GS-3583. Simulations were conducted for (A,D) PK, (B,E) cDC1 cells, and (C,F) total DC cells.
Figure 6. Projected time profiles of clinical PK, cDC1, and total DC cells in healthy volunteers treated with FLT3L-Fc. Simulations were performed to predict clinical FLT3L-Fc PK/PD profiles and compare cDC1 and total DC expansion under both cyno-derived and human-derived scenarios for a q3w dosing regimen. (AC) In the first row, FLT3 target-related parameters were estimated from cynomolgus monkey PK, and (DF) in the second row, FLT3 target-related parameters were estimated using human PK data from the FLT3 agonist Fc fusion protein molecule, GS-3583. Simulations were conducted for (A,D) PK, (B,E) cDC1 cells, and (C,F) total DC cells.
Pharmaceutics 16 00660 g006
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Hosseini, I.; Fleisher, B.; Getz, J.; Decalf, J.; Kwong, M.; Ovacik, M.; Bainbridge, T.W.; Moussion, C.; Rao, G.K.; Gadkar, K.; et al. 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. Pharmaceutics 2024, 16, 660. https://doi.org/10.3390/pharmaceutics16050660

AMA Style

Hosseini I, Fleisher B, Getz J, Decalf J, Kwong M, Ovacik M, Bainbridge TW, Moussion C, Rao GK, Gadkar K, et al. 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. Pharmaceutics. 2024; 16(5):660. https://doi.org/10.3390/pharmaceutics16050660

Chicago/Turabian Style

Hosseini, Iraj, Brett Fleisher, Jennifer Getz, Jérémie Decalf, Mandy Kwong, Meric Ovacik, Travis W. Bainbridge, Christine Moussion, Gautham K. Rao, Kapil Gadkar, and et al. 2024. "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" Pharmaceutics 16, no. 5: 660. https://doi.org/10.3390/pharmaceutics16050660

APA Style

Hosseini, I., Fleisher, B., Getz, J., Decalf, J., Kwong, M., Ovacik, M., Bainbridge, T. W., Moussion, C., Rao, G. K., Gadkar, K., Kamath, A. V., & Ramanujan, S. (2024). 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. Pharmaceutics, 16(5), 660. https://doi.org/10.3390/pharmaceutics16050660

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop