A flexible nonlinear feedback system that captures diverse patterns of adaptation and rebound

Developing a mechanistic translational PK/PD model for a trifunctional NK cell engager to predict the first-in-human dose for acute myeloid leukemia

Natural killer cell engagers (NKCEs), a treatment that stimulates innate immunity, have lately gained attention owing to their favorable safety profile, and their efficacy. Natural killer (NK) cell activation is driven by immune synapse formation between drugs, NK cells, and tumor cells. However, no clear translational modeling approach has been reported for first-in-human (FIH) dose estimation of humanized NKCEs. We developed the first translational mechanistic synapse-driven pharmacokinetic/pharmacodynamic (PK/PD) model for a trifunctional NKp46/CD16a-CD123 (CD123-NKCE) by integrating (i) in vitro target cell cytotoxicity in MOLM-13 tumor cell lines at varying effector-to-tumor cell ratios and incubation intervals; (ii) nonhuman primate PK and profiles of CD123+ cells and NKP46+ NK cells; and (iii) healthy human or patients with acute myeloid leukemia system-specific parameters. To depict direct tumor cell killing by the innate immunity, no transit compartment was included in PK/PD model structures. Model predictions suggested an intrapatient dose escalation of 10/30/100 μg/kg twice weekly to be selected as the starting dose in the FIH trial. However, sensitivity analyses revealed that CD123+ cell growth rate constant and maximal tumor killing rate constant were the key uncertainties to the recommended active dose. This novel translational model structure can be used as the basis to predict clinical PK/PD data for CD123-NKCE, and the translational strategy may serve as a foundation for future advancements of NKCEs.

Population pharmacokinetic/pharmacodynamic models of JNJ-64794964, a toll-like receptor 7 agonist, in healthy adult participants

BACKGROUND

JNJ-4964 is a TLR7 agonist, which, via a type I interferon (IFN)-dependent mechanism, may enhance host immunity suppressed by persistent exposure to hepatitis B antigens in chronic hepatitis B.

METHODS

PK and PD data were pooled from 2 studies involving 90 participants (n = 74 JNJ-4964, dose range 0.2-1.8 mg; n = 16 placebo) in a fasted state. Food effects on PK were studied in 24 participants (1.2 or 1.25 mg). A population PK model and PK/PD models were developed to characterize the effect of JNJ-4964 plasma levels on the time course of IFN-α, IFN-γ-inducible protein 10 (IP-10 or CXCL10), IFN-stimulated gene 15 (ISG15), neopterin and lymphocytes following single and weekly dosing in healthy adults. Covariate effects, circadian rhythms and negative feedback were incorporated in the models.

RESULTS

A 3-compartment linear PK model with transit absorption adequately described JNJ-4964 PK. Bioavailability was 44.2% in fed state relative to fasted conditions. Indirect response models with maximum effect (Emax) stimulation on production rate constant (kin) described IFN-α, IP-10, ISG15 and neopterin, while a precursor-dependent indirect response model with inhibitory effect described the transient lymphocyte reduction. Emax, EC50 and γ (steepness) estimates varied according to PD markers, with EC50 displaying substantial between-subject variability. Female and Asian race exhibited lower EC50, suggesting higher responsiveness.

CONCLUSIONS

PK/PD models well characterized the time course of immune system markers in healthy adults. Our results supported sex and race as covariates on JNJ-4964 responsiveness, as well as circadian rhythms and negative feedback as homeostatic mechanisms that are relevant in TLR7-induced type I IFN responses.

Application of the relationship between pharmacokinetics and pharmacodynamics in drug development and therapeutic equivalence: a PEARRL review

Objectives

The objective of this review was to provide an overview of pharmacokinetic/pharmacodynamic (PK/PD) models, focusing on drug-specific PK/PD models and highlighting their value added in drug development and regulatory decision-making.

Key findings

Many PK/PD models, with varying degrees of complexity and physiological understanding have been developed to evaluate the safety and efficacy of drug products. In special populations (e.g. paediatrics), in cases where there is genetic polymorphism and in other instances where therapeutic outcomes are not well described solely by PK metrics, the implementation of PK/PD models is crucial to assure the desired clinical outcome. Since dissociation between the pharmacokinetic and pharmacodynamic profiles is often observed, it is proposed that physiologically based pharmacokinetic and PK/PD models be given more weight by regulatory authorities when assessing the therapeutic equivalence of drug products.

Summary

Modelling and simulation approaches already play an important role in drug development. While slowly moving away from ‘one-size fits all’ PK methodologies to assess therapeutic outcomes, further work is required to increase confidence in PK/PD models in translatability and prediction of various clinical scenarios to encourage more widespread implementation in regulatory decision-making.

Characterizing the dynamic interaction among gastric emptying, glucose absorption, and glycemic control in nondiabetic obese adults

The effects of altered gastric emptying on glucose absorption and kinetics are not well understood in nondiabetic obese adults. The aim of this work was to develop a physiology-based model that can characterize and compare interactions among gastric emptying, glucose absorption, and glycemic control in nondiabetic obese and lean healthy adults. Dynamic glucose, insulin, and gastric emptying (measured with breath test) data from 12 nondiabetic obese and 12 lean healthy adults were available until 180 min after an oral glucose tolerance test (OGTT) with 10, 25, and 75 g of glucose. A physiology-based model was developed to characterize glucose kinetics applying nonlinear mixed-effects modeling with NONMEM7.3. Glucose kinetics after OGTT was described by a one-compartment model with an effect compartment to describe delayed insulin effects on glucose clearance. After the interactions between individual gastric emptying and glucose absorption profiles were accounted for, the glucose absorption rate was found to be similar in nondiabetic obese and lean controls. Baseline glucose concentration was estimated to be only marginally higher in nondiabetic obese subjects (4.9 vs. 5.2 mmol/l), whereas insulin-dependent glucose clearance in nondiabetic obese subjects was found to be cut in half compared with lean controls (0.052 vs. 0.029 l/min) and the insulin concentration associated with 50% of insulin-dependent glucose elimination rate was approximately twofold higher in nondiabetic obese subjects compared with lean controls (7.1 vs. 15.3 μU/ml). Physiology-based models can characterize and compare the dynamic interaction among gastric emptying, glucose absorption and glycemic control in populations of interest such as lean healthy and nondiabetic obese adults.

Feedback control indirect response models

A general framework is introduced for modeling pharmacodynamic processes that are subject to autoregulation, which combines the indirect response (IDR) model approach with methods from classical feedback control of engineered systems. The canonical IDR models are modified to incorporate linear combinations of feedback control terms related to the time course of the difference (the error signal) between the pharmacodynamic response and its basal value. Following the well-established approach of traditional engineering control theory, the proposed feedback control indirect response models incorporate terms proportional to the error signal itself, the integral of the error signal, the derivative of the error signal or combinations thereof. Simulations are presented to illustrate the types of responses produced by the proposed feedback control indirect response model framework, and to illustrate comparisons with other PK/PD modeling approaches incorporating feedback. In addition, four examples from literature are used to illustrate the implementation and applicability of the proposed feedback control framework. The examples reflect each of the four mechanisms of drug action as modeled by each of the four canonical IDR models and include: selective serotonin reuptake inhibitors and extracellular serotonin; histamine H2-receptor antagonists and gastric acid; growth hormone secretagogues and circulating growth hormone; β2-selective adrenergic agonists and potassium. The proposed feedback control indirect response approach may serve as an exploratory modeling tool and may provide a bridge for development of more mechanistic systems pharmacology models.

Challenges of a mechanistic feedback model describing nicotinic acid-induced changes in non-esterified fatty acids in rats

Previously, we developed a feedback model to describe the tolerance and oscillatory rebound of non-esterified fatty acid (NEFA) plasma concentrations in male Sprague Dawley rats after intravenous infusions of nicotinic acid (NiAc). This study challenges that model, using the following regimens of intravenous and oral NiAc dosing in male Sprague Dawley rats (n = 95) to create different patterns of exposure: (A) 30 min infusion at 0, 1, 5 or 20 μmol kg(-1) body weight; (B) 300 min infusion at 0, 5, 10 or 51 μmol kg(-1); (C) 30 min infusion at 5 μmol kg(-1), followed by a stepwise decrease in rate every 10 min for 180 min; (D) 30 min infusion at 5 μmol kg(-1), followed by a stepwise decrease in rate every 10 min for 180 min and another 30 min infusion at 5 μmol kg(-1) from 210 to 240 min; (E) an oral dose of 0, 24.4, 81.2 or 812 μmol kg(-1). Serial arterial blood samples were taken for measurement of plasma NiAc and NEFA concentrations. The gradual decrease in infusion rate in (C) and (D) were also designed to test the hypothesis that a gradual reduction in NiAc plasma concentration may be expected to reduce or prevent rebound. The absorption of NiAc was described by parallel linear and non-linear processes and the disposition of NiAc by a two-compartment model with endogenous turnover rate and two parallel capacity-limited elimination processes. NEFA (R) turnover, which was driven by the plasma concentration of NiAc via an inhibitory drug-mechanism function acting on NEFA formation, was described by a feedback model with a moderator distributed over a series of transit compartments, where the first compartment (M 1) inhibited the formation of R and the last compartment (M N ) stimulated the loss of R. All processes regulating the plasma NEFA concentration were assumed to be captured by the moderator function. Data were analyzed using non-linear mixed effects modeling (NONMEM). The potency IC 50 of NiAc was 68 nmol L(-1), the fractional turnover rate k out 0.27 L mmol(-1) min(-1), and the turnover rate of moderator k tol 0.023 min(-1). The lower physiological limit of NEFA, which was modeled as a NiAc-independent release (k cap ) of NEFA into plasma, was estimated to 0.023 mmol L(-1) min(-1). The parameter estimates derived in this study were consistent with our previous estimates, suggesting that the model may be used for prediction of the NEFA response time-course following different modes and routes administration of NiAc or NiAc analogues. In order to avoid NiAc-induced NEFA rebound, a slow decline in the NiAc exposure pattern is needed at or below IC (50).

Quantitative analysis of rate and extent of tolerance of biomarkers: application to nicotinic acid-induced changes in non-esterified fatty acids in rats

In this paper we quantitatively evaluate two feedback systems with a focus on rate and extent of tolerance and rebound development. In the two feedback systems, the regulation of turnover of response is governed by one or several moderators. In the basic system, one single moderator inhibits the formation of response. This system has been applied to cortisol secretion and serotonin reuptake inhibition. The basic system has been extended to adequately describe nicotinic acid (NiAc)-induced changes in non-esterified fatty acids (NEFA). In the extended system, the feedback is described by a cascade of moderators where the first inhibits formation of response and the last stimulates loss of response. The objectives of this paper were to analyze these systems from a mathematical/analytical and quantitative point of view and to present simulations with different parameter settings and dosing regimens in order to highlight the intrinsic behaviour of these systems and to present expressions and graphs that are applicable for quantification of rate and extent of tolerance and rebound. The dynamics of the moderators (k(tol)) compared to the dynamics of the response (k(out)), was shown to be important for the behaviour of both systems. For instance, slow dynamics of the moderator compared to the response (k(tol)<<k(out)), resulted in overshoot and pronounced rebound. The extent of tolerance was studied over time at a single constant drug concentration and at steady state for different drug concentrations and was found to be largest at drug concentrations close to IC(50). An upper limit for the response could be identified and included in the quantification of extent of rebound. Especially, for the extended system, the duration of exposure was an important factor affecting size of rebound. The rate of tolerance development was addressed by quantitatively estimating the time to steady state for the two systems, in which the value of k(tol) and the length of the cascade were critical.

Feedback modeling of non-esterified fatty acids in rats after nicotinic acid infusions

A feedback model was developed to describe the tolerance and oscillatory rebound seen in non-esterified fatty acid (NEFA) plasma concentrations following intravenous infusions of nicotinic acid (NiAc) to male Sprague-Dawley rats. NiAc was administered as an intravenous infusion over 30 min (0, 1, 5 or 20 μmol kg⁻¹ of body weight) or over 300 min (0, 5, 10 or 51 μmol kg⁻¹ of body weight), to healthy rats (n = 63), and serial arterial blood samples were taken for measurement of NiAc and NEFA plasma concentrations. Data were analyzed using nonlinear mixed effects modeling (NONMEM). The disposition of NiAc was described by a two-compartment model with endogenous turnover rate and two parallel capacity-limited elimination processes. The plasma concentration of NiAc was driving NEFA (R) turnover via an inhibitory drug-mechanism function acting on the formation of NEFA. The NEFA turnover was described by a feedback model with a moderator distributed over a series of transit compartments, where the first compartment (M₁) inhibited the formation of R and the last compartment (M_N) stimulated the loss of R. All processes regulating plasma NEFA concentrations were assumed to be captured by the moderator function. The potency, IC₅₀, of NiAc was 45 nmol L⁻¹, the fractional turnover rate k_out was 0.41 L mmol⁻¹ min⁻¹ and the turnover rate of moderator k_tol was 0.027 min⁻¹. A lower physiological limit of NEFA was modeled as a NiAc-independent release (k_cap) of NEFA into plasma and was estimated to 0.032 mmol L⁻¹ min⁻¹. This model can be used to provide information about factors that determine the time-course of NEFA response following different modes, rates and routes of administration of NiAc. The proposed model may also serve as a preclinical tool for analyzing and simulating drug-induced changes in plasma NEFA concentrations after treatment with NiAc or NiAc analogues.

Turnover modeling of non-esterified fatty acids in rats after multiple intravenous infusions of nicotinic acid

The objective of this investigation was to use a pharmacokinetic (PK)/pharmacodynamic (PD) approach to describe and evaluate a PK model of nicotinic acid (NiAc) in guinea pigs and a PD feedback model of changes in non-esterified fatty acid (NEFA) concentrations in rats following multiple intravenous infusions of NiAc at different rates and durations of inhouse and literature (NEFA after extravascular NiAc dosing) data. Serial arterial blood samples were taken for evaluation of NiAc exposure in guinea pigs and NEFA in rats. The biophase kinetics of NiAc was assumed to impact on NEFA turnover with feedback incorporated via an inhibitory moderator compartment. The response acted linearly on the production of moderator, which then acted inversely on the turnover rate of response. The potency, expressed as the amount of NiAc in the biophase causing a 50 % inhibitory effect (ID(50)), was 6.5 nmol +/- 31 % and the half-life of response (t(1/2, kout)) 2 min +/- 18 %. The half-life of tolerance (t(1/2, ktol)) was 9 min +/- 27 %. The model can be used to provide information about factors that determine the time course of NEFA response following different rates and routes of administration of NiAc or NiAc analogues.