Emapalumab in primary haemophagocytic lymphohistiocytosis and the pathogenic role of interferon gamma: A pharmacometric model-based approach.

AIM: Primary haemophagocytic lymphohistiocytosis (HLH) is a rare, life-threatening, hyperinflammatory syndrome generally occurring in early childhood. The monoclonal antibody emapalumab binds and neutralises interferon γ (IFNγ). This study aimed to determine an emapalumab dosing regimen when traditional dose-finding approaches are not applicable, using pharmacokinetic-pharmacodynamic analyses to further clarify HLH pathogenesis and confirm IFNγ neutralisation as the relevant therapeutic target in pHLH. METHODS: Initial emapalumab dosing (1 mg/kg) for pHLH patients participating in a pivotal multicentre, open-label, single-arm, phase 2/3 study was based on anticipated IFNγ levels and allometrically scaled pharmacokinetic parameters estimated in healthy volunteers. Emapalumab dosing was adjusted based on estimated IFNγ-mediated clearance and HLH clinical and laboratory criteria. Frequent dosing and emapalumab dose adaptation were used to account for highly variable IFNγ levels and potential target-mediated drug disposition. RESULTS: High inter- and intra-individual variability in IFNγ production (assessed by total IFNγ levels, range: 102 -106  pg/mL) was observed in pHLH patients. Administering emapalumab reduced IFNγ activity, resulting in significant improvements in clinical and laboratory parameters and a reduced risk of adverse events, mainly related to pHLH. Modelled outcomes supported dose titration starting from 1 mg/kg, with possible increases to 3, 6 or 10 mg/kg based on re-evaluation of parameters of disease activity every 3 days. CONCLUSIONS: The variable and unanticipated extremely high IFNγ concentrations in patients with pHLH are reflected in parameters of disease activity. Improved outcomes can be achieved by neutralising IFNγ using frequent emapalumab dosing and dose adaptation guided by clinical and laboratory observations.

Quantification of Radiation Injury on Neutropenia and the Link between Absolute Neutrophil Count Time Course and Overall Survival in Nonhuman Primates Treated with G-CSF.

PURPOSE: To model absolute neutrophil count (ANC) suppression in response to acute radiation (AR) exposure and evaluate ANC time course as a predictor of overall survival (OS) in response to AR exposure with or without treatment with granulocyte colony-stimulating factor in nonhuman primates. METHODS: Source data were obtained from two pivotal studies conducted in rhesus macaques exposed to 750 cGy of whole body irradiation on day 0 that received either placebo, daily filgrastim, or pegfilgrastim (days 1 and 8 after irradiation). Animals were observed for 60 days with ANC measured every 1 to 2 days. The population model of ANC response to AR and the link between observed ANC time course and OS consisted of three submodels characterizing injury due to radiation, granulopoiesis, and a time-to-event model of OS. RESULTS: The ANC response model accurately described the effects of AR exposure on the duration of neutropenia. ANC was a valid surrogate for survival because it explained 76% (95% CI, 41%-97%) and 73.2% (95% CI, 38.7%-99.9%) of the treatment effect for filgrastim and pegfilgrastim, respectively. CONCLUSION: The current model linking radiation injury to neutropenia and ANC time course to OS can be used as a basis for translating these effects to humans.

Emapalumab in Children with Primary Hemophagocytic Lymphohistiocytosis.

BACKGROUND: Primary hemophagocytic lymphohistiocytosis is a rare syndrome characterized by immune dysregulation and hyperinflammation. It typically manifests in infancy and is associated with high mortality. METHODS: We investigated the efficacy and safety of emapalumab (a human anti-interferon-γ antibody), administered with dexamethasone, in an open-label, single-group, phase 2-3 study involving patients who had received conventional therapy before enrollment (previously treated patients) and previously untreated patients who were 18 years of age or younger and had primary hemophagocytic lymphohistiocytosis. The patients could enter a long-term follow-up study until 1 year after allogeneic hematopoietic stem-cell transplantation or until 1 year after the last dose of emapalumab, if transplantation was not performed. The planned 8-week treatment period could be shortened or extended if needed according to the timing of transplantation. The primary efficacy end point was the overall response, which was assessed in the previously treated patients according to objective clinical and laboratory criteria. RESULTS: At the cutoff date of July 20, 2017, a total of 34 patients (27 previously treated patients and 7 previously untreated patients) had received emapalumab; 26 patients completed the study. A total of 63% of the previously treated patients and 65% of the patients who received an emapalumab infusion had a response; these percentages were significantly higher than the prespecified null hypothesis of 40% (P = 0.02 and P = 0.005, respectively). In the previously treated group, 70% of the patients were able to proceed to transplantation, as were 65% of the patients who received emapalumab. At the last observation, 74% of the previously treated patients and 71% of the patients who received emapalumab were alive. Emapalumab was not associated with any organ toxicity. Severe infections developed in 10 patients during emapalumab treatment. Emapalumab was discontinued in 1 patient because of disseminated histoplasmosis. CONCLUSIONS: Emapalumab was an efficacious targeted therapy for patients with primary hemophagocytic lymphohistiocytosis. (Funded by NovImmune and the European Commission; NI-0501-04 and NI-0501-05 ClinicalTrials.gov numbers, NCT01818492 and NCT02069899.).

Clustering and Kernel Density Estimation for Assessment of Measurable Residual Disease by Flow Cytometry.

Standardization, data mining techniques, and comparison to normality are changing the landscape of multiparameter flow cytometry in clinical hematology. On the basis of these principles, a strategy was developed for measurable residual disease (MRD) assessment. Herein, suspicious cell clusters are first identified at diagnosis using a clustering algorithm. Subsequently, automated multidimensional spaces, named "Clouds", are created around these clusters on the basis of density calculations. This step identifies the immunophenotypic pattern of the suspicious cell clusters. Thereafter, using reference samples, the "Abnormality Ratio" (AR) of each Cloud is calculated, and major malignant Clouds are retained, known as "Leukemic Clouds" (L-Clouds). In follow-up samples, MRD is identified when more cells fall into a patient's L-Cloud compared to reference samples (AR concept). This workflow was applied on simulated data and real-life leukemia flow cytometry data. On simulated data, strong patient-dependent positive correlation (R2 = 1) was observed between the AR and spiked-in leukemia cells. On real patient data, AR kinetics was in line with the clinical evolution for five out of six patients. In conclusion, we present a convenient flow cytometry data analysis approach for the follow-up of hematological malignancies. Further evaluation and validation on more patient samples and different flow cytometry panels is required before implementation in clinical practice.

Prediction of Survival Benefit of Filgrastim in Adult and Pediatric Patients With Acute Radiation Syndrome.

Acute exposure to high doses of radiation leads to severe myelosuppression, but few treatments are currently available to treat hematopoietic syndrome of acute radiation syndrome. Granulocyte colony stimulating factors (e.g., filgrastim) stimulate proliferation of neutrophil precursors and enhance mature neutrophil function. Owing to ethical constraints on conducting clinical research in lethally irradiated humans, we developed a model-based strategy to integrate preclinical experience in irradiated nonhuman primates (NHPs) and other clinical myelosuppressive conditions to inform filgrastim dosing to treat hematopoietic syndrome of acute radiation syndrome. Models predicting neutrophil counts and overall survival based on drug exposures were calibrated and scaled from NHPs to adult and pediatric human subjects. Several scenarios were examined investigating variations in filgrastim doses, dose frequency, treatment initiation, and duration, as well as the effect of age and radiation dose rate. Model-based simulations and established safety profiles supported that a subcutaneous filgrastim dose of 10 µg/kg once daily provides a significant survival benefit (50%) over placebo in both adults and children, provided that the treatment is initiated within 1-14 days after radiation exposure and lasts 2-3 weeks. For treatment durations of longer than 3 weeks, filgrastim treatment is not expected to provide significantly greater benefit. This survival benefit is expected to hold for the wide range of radiation doses and dose rates (0.01-1,000 Gy/hours) examined.

Pharmacokinetic-pharmacodynamic modelling of neutrophil response to G-CSF in healthy subjects and patients with chemotherapy-induced neutropenia.

AIM: The objective of the present study was to use pharmacokinetic-pharmacodynamic modelling to characterize the effects of chemotherapy on the granulopoietic system and to predict the absolute neutrophil counts (ANCs) for patients with chemotherapy-induced neutropenia treated with filgrastim and pegfilgrastim. METHODS: Data were extracted from 10 phase I-III studies conducted in 110 healthy adults, and 618 adult and 52 paediatric patients on chemotherapy following administration of filgrastim or pegfilgrastim. The structural model accounted for ANC dynamics and the effects of filgrastim and pegfilgrastim, chemotherapy and corticosteroids. The impact of neutrophils on drug disposition was based on a drug receptor-binding model that assumed quasi-equilibrium and stimulation of the production and maturation of neutrophils upon treatment. The chemotherapy and corticosteroid effects were represented by kinetic-pharmacodynamic-type models, where chemotherapy stimulated elimination of neutrophil precursors at the mitotic stage, and corticosteroids stimulated neutrophil production. RESULTS: The systemic half-lives of filgrastim (2.6 h) and pegfilgrastim (10.1 h) were as expected. The effective half-life of chemotherapy was 9.6 h, with a 2-day killing effect. The rate of receptor elimination from mitotic compartments exhibited extreme interindividual variability (% coefficient of variation >200), suggesting marked differences in sensitivity to chemotherapy effects on ANCs. The stimulatory effects of pegfilgrastim were significantly greater than those of filgrastim. Model qualification confirmed the predictive capability of this model. CONCLUSION: This qualified model simulates the time course of ANC in the absence or presence of chemotherapy and predicts nadir, time to nadir and time of recovery from different grades of neutropenia upon treatment with filgrastim and pegfilgrastim.

The use of the SAEM algorithm in MONOLIX software for estimation of population pharmacokinetic-pharmacodynamic-viral dynamics parameters of maraviroc in asymptomatic HIV subjects.

Using simulated viral load data for a given maraviroc monotherapy study design, the feasibility of different algorithms to perform parameter estimation for a pharmacokinetic-pharmacodynamic-viral dynamics (PKPD-VD) model was assessed. The assessed algorithms are the first-order conditional estimation method with interaction (FOCEI) implemented in NONMEM VI and the SAEM algorithm implemented in MONOLIX version 2.4. Simulated data were also used to test if an effect compartment and/or a lag time could be distinguished to describe an observed delay in onset of viral inhibition using SAEM. The preferred model was then used to describe the observed maraviroc monotherapy plasma concentration and viral load data using SAEM. In this last step, three modelling approaches were compared; (i) sequential PKPD-VD with fixed individual Empirical Bayesian Estimates (EBE) for PK, (ii) sequential PKPD-VD with fixed population PK parameters and including concentrations, and (iii) simultaneous PKPD-VD. Using FOCEI, many convergence problems (56%) were experienced with fitting the sequential PKPD-VD model to the simulated data. For the sequential modelling approach, SAEM (with default settings) took less time to generate population and individual estimates including diagnostics than with FOCEI without diagnostics. For the given maraviroc monotherapy sampling design, it was difficult to separate the viral dynamics system delay from a pharmacokinetic distributional delay or delay due to receptor binding and subsequent cellular signalling. The preferred model included a viral load lag time without inter-individual variability. Parameter estimates from the SAEM analysis of observed data were comparable among the three modelling approaches. For the sequential methods, computation time is approximately 25% less when fixing individual EBE of PK parameters with omission of the concentration data compared with fixed population PK parameters and retention of concentration data in the PD-VD estimation step. Computation times were similar for the sequential method with fixed population PK parameters and the simultaneous PKPD-VD modelling approach. The current analysis demonstrated that the SAEM algorithm in MONOLIX is useful for fitting complex mechanistic models requiring multiple differential equations. The SAEM algorithm allowed simultaneous estimation of PKPD and viral dynamics parameters, as well as investigation of different model sub-components during the model building process. This was not possible with the FOCEI method (NONMEM version VI or below). SAEM provides a more feasible alternative to FOCEI when facing lengthy computation times and convergence problems with complex models.

Population pharmacokinetics meta-analysis of recombinant human erythropoietin in healthy subjects.

OBJECTIVE: The aim of this analysis was to develop a population pharmacokinetic model to describe the pharmacokinetics of recombinant human erythropoietin (rHuEPO) in healthy subjects, after intravenous and subcutaneous administration over a wide dose range, and to examine the influence of demographic characteristics and other covariates on the pharmacokinetic parameters of rHuEPO. METHODS: Erythropoietin serum concentration data were available from 16 studies comprising 49 healthy subjects who received rHuEPO intravenous doses from 10 to 300 IU/kg, 427 healthy subjects who received rHuEPO subcutaneous doses from 1 to 2400 IU/kg, and 57 healthy subjects who received placebo and where endogenous erythropoietin concentrations were measured. Different pharmacokinetic models were fitted to the dataset using nonlinear mixed-effects modeling software (NONMEM, Version V, Level 1). Several patient covariates were tested in order to quantify the effect on rHuEPO pharmacokinetic parameters. Model evaluation was examined using a posterior predictive check. RESULTS: Erythropoietin showed a diurnal baseline variation of +/-20%, described with a dual cosine model. Disposition was described with a two-compartment model with a small volume of distribution (6L) and parallel linear and nonlinear clearance. Total clearance varied between 0.3 and 0.9 L/h over the concentration range studied. A dual absorption model was used to characterise the rHuEPO absorption from the subcutaneous formulation and consisted of a faster pathway described as a sequential zero- and first-order absorption process and a parallel slower pathway characterised as a zero-order process. The bioavailability of subcutaneous rHuEPO increased from 30% at low doses to 71% at the highest dose of 160 kIU and was described using a hyperbolic model. The most important covariate effects were a decrease in the first-order absorption rate constant (k(a)) with increasing age, an increase in subcutaneous bioavailability with increasing baseline haemoglobin, and a decrease in bioavailability with increasing bodyweight. A posterior predictive check showed no systematic deviation of the simulated data from the observed values. CONCLUSION: The population pharmacokinetic model developed is suitable to describe the pharmacokinetic behaviour of rHuEPO after intravenous and subcutaneous administration in healthy subjects, over a wide dose range.

A semimechanistic and mechanistic population PK-PD model for biomarker response to ibandronate, a new bisphosphonate for the treatment of osteoporosis.

AIMS: Ibandronate, a highly potent nitrogen-containing bisphosphonate, is the subject of an ongoing clinical development programme that aims to maximize the potential of simplified, less frequent oral and intravenous (i.v.) administration in osteoporosis. A modelling and simulation project was undertaken to characterize further the clinical pharmacology of ibandronate and identify convenient intermittent oral and i.v. regimens for clinical evaluation. METHODS AND RESULTS: Using selected data from clinical studies involving 174 women with postmenopausal osteoporosis (PMO), a classical multicompartmental pharmacokinetic-pharmacodynamic (PK-PD) model was developed that accurately described the PK of i.v. ibandronate in plasma and urine and urinary excretion of the C-telopeptide of the alpha chain of type I collagen (uCTX), a sensitive biomarker of PD response to ibandronate. To reduce processing times, the classical PK-PD model was simplified using a "kinetics of drug action" or kinetic (K)-PD model (i.e. a dose-response model as opposed to a dose-concentration-response model). The performance of the K-PD model was evaluated by fitting data simulated with the PK-PD model under various dosing regimens. The simplified model produced a virtually indistinguishable fit of the data from that of the PK-PD model. The K-PD model was extended to consider the influence of supplemental therapy (calcium with or without vitamin D) on the PD response and validated by retrospectively simulating the uCTX response in a prior Phase III and Phase II/III study of i.v. ibandronate, given once every 3 months, in 3380 women with PMO. The observed median uCTX responses at the scheduled assessment points in the completed studies were within the distribution of the simulated responses. The K-PD model for i.v. ibandronate was extended further to allow simultaneous fitting of uCTX responses after i.v. and oral administration in 676 postmenopausal women with osteoporosis, and validated by retrospectively simulating the data observed in a Phase I study of oral daily ibandronate in 180 women with PMO. The K-PD model adequately described the uCTX response after oral dosing. CONCLUSIONS: This validated K-PD model is currently being used to evaluate a range of novel intermittent oral and i.v. ibandronate regimens in an ongoing clinical development programme.